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contributions 


The Effect of Echo on the Operation of High-Frequency 
Communication Circuits* 
D. K. BAILEY} 


Summary—Echo on high-frequency communication services is 


defined as the simultaneous reception of signals over both major and 


minor arcs of the great circle connecting a transmitter and receiver. 
Two distinct kinds of echo are recognized according to the illumina- 


tion conditions under which they occur. Echo of the first kind is 


observed when the great-circle path coincides with the twilight zone 


surrounding the earth, whereas echo of the second kind, which can 
occur only on fairly long communication paths, is most severe when 


i ts ee 


the short path is most intensely illuminated. Little can be done to 


obviate echo of the first kind and it is not, like echo of the second 
kind, amenable to prediction by available methods of calculating 


_ sky-wave field intensities. Radio traffic data are cited which cor- 


tf 


- roborate calculations and show both that echo of the second kind is 
more severe at times of maximum solar activity, and is less severe 


on higher frequencies. Conclusions are drawn concerning mode 
of operation and choice of operating frequency to minimize echo 


i 


3 
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- The material in this paper constitutes 


ee. . Te 


- tions Union, t 48 t 
‘It reflects the international character of the Geneva activities. The 


- interference. 


INTRODUCTION 
Rim engineers operating high-frequency com- 


munication circuits over long distances have for 
many years identified a certain type of difficulty 


* Manuscript received by the PGAP, August 12, 1957. Originally 
issued as Natl. Bur. of Standards report (1951), Washington, D. C. 
one of the necessary prelim- 
inary investigations undertaken by the author in connection with 
his duties as Chairman of the Propagation Working Group of the 
Provisional Frequency Board of the International Telecommunica- 
which met at Geneva, Switzerland, from 1948 to 1950. 


_ Propagation Working Group was specifically charged with the duty 


aie in” 


A es 


of establishing on an engineering basis rules for the selection of 


- operating frequencies for high-frequency point-to-point radio com- 
‘munications. ‘ 


Page Communications Engineers, Inc., Washington, D. C. 


Formerly with Natl. Bur. of Standards, Washington, D. C. 


in the operation of their circuits with the simultaneous 
reception of signals over both the major arc (known as 
the long path) and the minor arc (known as the short 
path) of the great circle connecting the transmitter and 
the receiver. The identity of the long-path signals is 
easily established by the determination of their delay 
times with respect to the short-path signals. This re- 
ception condition is known as echo because of the au- 
rally perceptible delays in telegraphy and the hollow 
“large auditorium” sound of telephony. Echo should 
not be confused with multipath distortion of high-speed 
telegraphic communication which is the result of the 
difference in arrival time between the first and last 
significant modes of propagation between transmitter 
and receiver along the same portion of the entire great 
circle. If the very difficult case of antipodal circuits 
(circuits of the order of 18,000 km or longer in length) 
is excluded, it may be said that the simple echo condi- 
tion described above is characterized by delays of the 
order of 10 to 100 msec. Such delays are capable of pre- 
venting any but slow-speed telegraphy and of causing 
annoyance to users of telephony. 

In the course of an investigation of echo and its re- 
lationship to the mode of operation and the choice of 
operating frequencies, it was found to be possible to dis- 
tinguish two quite distinct types of echo, according to 
the conditions under which they occur." 


1 “Second Report of Working Group 1 (Propagation) Committee 
4,” Internatl. Telecommun. Union, Provisional Frequency Board, 
Geneva, Switzerland, Doc. No, 375-E, pp. 4-5; October 6, 1948. 
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EcHo OF THE First KIND 


Echo of the first kind is observed on circuits having 

lengths as short as 5000 km and occurs only at times 
when the great-circle path connecting the transmitter 
and receiver is coincident or nearly coincident with the 
twilight great circle which surrounds the earth. Under 
this particular illumination condition the attenuation 
of radio waves is at times much less than can be ac- 
counted for by use of any of the presently available 
quasi-empirical tools for the computation of skywave 
field intensities.2 The broken-line curve in Fig. 2, con- 
sidered in the light of the discussion of the next section, 
illustrates this failure in a specific manner. In fact nu- 
merous observations are known which have revealed ob- 
servable signals that traveled once and occasionally 
several times entirely around the earth in both direc- 
tions in the twilight zone. For observers located in the 
skip-zone of nearby high-powered transmitters, a special 
case of twilight-zone echo is sometimes observed with 
delays approximating 4 of a second. It is the result of 
the time delay between the scattered signal (often arriv- 
ing from some direction other than the true direction of 
the transmitter), or a weak direct signal, and the signal 
which has passed entirely around the earth in the twi- 
light zone. Hess’ has recently published results of 
studies of the twilight-zone type of echo, and provides 
useful references to early reports. He and the earlier 
workers appear to have observed echo only of the first 
kind. This can probably be ascribed to the fact that 
their observations were made at times of intermediate 
and low solar activity, a circumstance which would ac- 
count for their belief that echo could only occur under 
twilight-zone conditions. 
_ There is at present insufficient understanding of echo 
of the first kind to permit an evaluation of its effect on 
the selection of operating frequencies for high-frequency 
communication. It seems fairly clear however that echo 
of the first kind is independent of transmitter power, and 
that improving the front-to-back ratio of the transmit- 
ting and receiving antennas in view of the complete 
round-the-world possibilities can at best be of limited 
use in obviating it as a source of cochannel interference. 
Echo of the first kind can occur with serious conse- 
quences at all phases of the solar cycle. 


ECHO OF THE SECOND KIND 


Echo of the second kind is observed on circuits having 
short-path lengths greater than 8000 to 10,000 km at 
times when the short path is most intensely illuminated. 
It seems to have its most serious consequences to the 
communications services when the sun is near or at the 
zenith at the midpoint of the short path. For a particu- 


2 “Tonospheric Radio Propagation,” Natl. Bur. of Standards, 
Washington, D. C., Circular 462, pp. 104-150; June 25, 1948, 

3H. A. Hess, “Investigations of high-frequency echoes,” Proc, 
IRE, vol. 36, pp. 981-992; August, 1948. 

4H. A. Hess, “Part 1I—Investigations of high-frequency echoes,” 
Proc. IRE, vol. 37, pp. 986-989; September, 1949, 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


+70 


+60 


ANGULAR ELEVATION 
OF SUN 90° AT 
MIDPOINT OF SHORT PATH 


15 Mc 


+50 


NORMAL PROPAGATION 


+30 


+20 T 


(db) 


+10 


SS 


-20 


ag seitelne Soe pak pm yi rR DR 


= 
' 
' 
! 
' 
nn 
nn 
on 
uw 
age uat ates dal 


Omi x 


| 


cad 


q 
. 


fe) a _—— 
4000 6000 8000 10000 12000 14000 16000 18000 20000 


SHORT-PATH LENGTH (km) 


Fig. 1—Calculated severity of echo of the second kind as produced 
under the most favorable illumination condition. : 


lar operating frequency it is more serious at periods of 
maximum than at periods of minimum solar activity. 

Presently available methods for the computation of 
long-term median sky-wave field intensities? permit the 
prediction of a number of characteristics of echo of the 
second kind. If normal propagation is assumed to be 
possible over both the long and short paths, it is pos- 
sible to calculate, for an omnidirectional transmitting 
antenna, the ratio of the median field intensities of the 
long to the short path at the receiving location. This 
ratio, called L/S, is expressed in Figs. 1 to 3 in decibels. 
It measures the severity of the echo condition. The main 
parameters involved in its calculation are: 


1) Length of short path, 

2) Frequency, . 

3) Geometrical position of the sun relative to the long 
and short paths (uniquely given by two angles, vzz. 
the angular elevation of the sun at the midpoint 
of the short path, and the orientation of the short 
path at its midpoint with respect to the minor arc 
of the great circle which connects the midpoint of 
the short path to the subsolar point), X 

4) Level of solar activity (as measured by the 
smoothed relative sunspot number, abbreviated 
“SS” on the figures), and 

5) Auroral-zone absorption on the long path (taken 
to be zero on all the curves except one on Fig. 3 


for which representative maximum auroral-zone 
effects have been introduced). 
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| Fig. 2—Calculated severity of echo of the second kind for a frequency 
of 20 mc for various illumination conditions. 


Y _ Then, considering a range of short-path lengths from 
8000 to 18,000 km, by varying each of the main param- 
‘eters in the order above in turn, it is seen in Figs. 1 to 


1) Increases in seriousness with length of short path, 

Varies in seriousness as an inverse function of 

frequency, 

Is most serious when the sun is overhead at the 

midpoint of the short path, 

4) Varies in seriousness as a direct function of solar 
activity, and 

5) Decreases in seriousness with increasing auroral 
disturbance to the long path. 


“In addition to the main points above, the following addi- 
ional points are also brought to light by the calcula- 
tions: 

6) There appears to exist a critical short-path length 
below which echo effects vary in seriousness di- 
rectly with frequency. This critical short-path 
length is a function of the geometrical position of 
the sun relative to the long and short paths, but is 
too small for serious echo effects on any frequency 
likely to be of practical use. 

7) The effect of introducing auroral absorption on 
the long path is to depress a particular curve of the 

type shown in Fig. 1 by an amount independent 
of short-path length but dependent upon fre- 
quency. Curves for lower frequencies are more de- 


—e 
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3 for a variety of typical absorption conditions and fre- _ 
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Fig. 3—Calculated severity of echo of the second kind for a frequency 
of 20 mc as produced under the most favorable illumination con- 


dition for representative levels of solar activity, including the 
auroral case. 


pressed. This would have the effect, for example 
in Fig. 1, of increasing the critical short-path 
length as discussed in the preceding point. 


None of these conclusions is valid for very low angles 
of elevation of the sun at the path midpoint, when, 
with suitable path orientation, the necessary condition 
for echo of the first kind is established. 


COMPARISON WITH OBSERVATIONS 


Direct evidence in support of the first five predicted 
characteristics of echo of the second kind, as well as 
examples of echo of the first kind, are available in re- 
ceiving station logs. 

Figs. 4 and 5 show the observed occurrence of echo 
without regard to kind on the circuit from London to 
Buenos Aires as observed at Buenos Aires on the prin- 
cipal day frequency used for this 11,000-km circuit for 
1947, a year of maximum solar activity, and for 1944, a 
year of minimum solar activity. The station being re- 
ceived is GLW operating on a frequency of 19,100 ke. 
The figures show, by the heights of the blocks, the rela- 
tive severity of the echo, if present, which was observed 
at each hour of the day for each month of the year. The 
units are arbitrary. Sunrise and sunset curves are shown 
for each terminal of the path. The occurrence of echo of 
the first kind is thus identifiable during the morning 
twilight most markedly in December and January, and 
during the evening twilight most markedly in June 
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1947 on signals from station GLW, England, 19,100 kc. 
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Fig. 5—Relative severity of echo observed at Buenos Aires during 
1944 on signals from station GLW, England, 19,100 ke. 


and July. At other seasons the twilight zone does not 
coincide very closely with the path at any time of 
day. During 1944 the operating frequency under con- 
sideration was sufficiently high to receive relatively lit- 
tle use at times when echo of the first kind could occur. 

With regard to echo of the second kind Figs. 4 and 5 
reveal that the effects are most severe at maximum 
solar activity, and of no serious consequence for the 
particular case shown at minimum. Furthermore the 
sun is most nearly vertical at the short-path midpoint 
during the months of April and August; the observa- 
tions of Fig. 4 reveal more serious echo of the second 
kind during these months, and show definite minima in 
the solstice months. As would be expected from the fore- 
going and from the actual geographical position of the 
short path, the minimum at the December solstice is 
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Fig. 6—Relative severity of echo observed at Buenos Aires during 
: 1947 on signals from station PGN, Netherlands, 24,465 kc. 
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Fig. 7—Relative severity of echo observed at Buenos Aires during 
1947 on signals from station WOJ, New York, 21,240 ke. 


the deeper. The improvement to be obtained by the use 
in 1947 of a higher frequency is illustrated by Fig. 6 
which presents the echo observations at Buenos Aires 
for a frequency of 24,465 kc as tran$mitted by station 
PGN during its use on the Amsterdam to Buenos Aires 
circuit. The paths London to Buenos Aires and Am- 
sterdam to Buenos Aires are considered sufficiently 
similar to permit direct comparison. It is unlikely that 
the full benefits to be derived from use of a higher fre- 
quency are revealed by the comparison of Figs. 4 and 6 
since the realizable front-to-back gain of the trans- 
mitting Franklin beam array at London is probably sig- 
nificantly greater than that of the rhombic used at 
Amsterdam. Fig. 7 presents the echo observation made 
during 1947 at Buenos Aires of station WQJ used for the 
New York to Buenos Aires circuit and operating on 


+1958 


21,240 ke. This path is about 8500 km in length and 
therefore significantly shorter than the London and 
Amsterdam paths. The long path from New York to 
Buenos Aires passes directly across both polar regions, 
and the observation of only very slight echo disturbance 
is clearly revealed. 


CONCLUSIONS 


Without definite and detailed information on such 
‘Matters as the value of the protection ratio for a par- 
ticular type of telegraph service against cochannel inter- 
ference, the benefits to be derived from the use of di- 
versity reception, and the values of realizable front-to- 
back gains of transmitting and receiving antennas, it is 
not profitable to discuss further the problem of echo of 
‘the second kind as it affects particular communications 
circuits. Nevertheless the following conclusions of a gen- 
eral nature may be expressed with regard to the mode of 
‘operation of high-frequency communications circuits: 


1) The interference effects of echo of either kind ina 
particular situation cannot be reduced by increas- 
ing the power, 

' 2) Increasing the realizable front-to-back gains of 
_ either or both the transmitting and receiving an- 

iS, tennas will reduce the interference effects of echo 

4 of the second kind, 

j 3) Introduction of diversity gain or improvement in 
diversity gain of the receiving antennas will reduce 
the interference effect of echo of the second kind, 
and 

Use of the long path, particularly for larger values 

of short-path length such as 12,000 to 18,000 km, 


directivity, will often give a usable service when 
the short path is unusable with the antennas and 
2 frequency available, because of echo of the second 
4 kind. 


with adequate transmitting and receiving antenna © 
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For circuits having short paths subject to interference 
from echo of the second kind, the following conclusions 
may be tentatively expressed with regard to the selec- 
tion of the operating frequency: 


5) 


6) 


It is important to select the highest possible day 
frequency for use during the minimum phases of 
the solar cycle in order to minimize echo inter- 
ference at these times. It would appear for most 
practical purposes that this frequency should lie 
between approximately 18 and 20 mc. 

An additional and higher day frequency should be 
provided, according to the propagational charac- 
teristics of the path, for use during years of maxi- 
mum solar activity when the day frequency de- 
scribed in the preceding conclusion suffers exces- 
sively from echo interference. It would appear for 
most practical purposes that this frequency should 
lie as far above 22 mc as practicable. For very long 
circuits, for example circuits having short paths 
longer than about 14,000 km, on which the long 
path is nonauroral should receive considerable use, 
this additional frequency is probably unnecessary. 
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Foreground Terrain Effects on Overland UHF Transmissions* 
L. G. TROLESE} anp L. J. ANDERSON} 


Summary—This paper describes an experimental study of the 
influence of the shape of foreground terrain profiles near terminals 
of UHF links on the received field. A gently rounded shape of the 
foreground profile causes a marked diffraction pattern to be super- 
imposed on the normal variation of field strength with height. The 
diffraction geometry shows similarity to knife-edge geometry and 
the rounded terrain feature appears to act geometrically as an equiva- 
lent knife edge. The amplitude of the spatial variations in signal are, 
however, much greater than knife-edge theory predicts. A sizeable 
foreground diffraction enhancement of received field can be realized 
by locating the antenna at the height of the first diffraction maxi- 
mum. Changing refraction due to meteorological variations can 
change both the position in height and intensity of the diffraction 
pattern. 


INTRODUCTION 


ERRAIN features which affect the propagation 
Woo UHF transmissions can be divided into several 
general classes. Large dominant obstacles such as 
mountains and mountain ridges have been studied in 
considerable detail, and many papers have appeared in 
the literature on this subject. The classical knife-edge 
theory has been used with reasonably good success in 
predicting the field to be expected in the shadow region 
behind such large obstacles. t 
Another class of terrain features which has received 
considerable attention is a terrain composed of irregu- 
larities which vary in a random way. For short wave- 
length transmissions the sea surface is usually con- 
sidered to be of this type, and most overland paths have 
irregularities which are certainly random in nature. 
This paper, however, is concerned with terrain features 
which are in neither of the two classes listed above. The 
experimental studies which are described concern the 
influence of the foreground terrain near the terminals of 
UHF links on the received field strength. The shape of 
the foreground terrain was found to have a marked in- 
fluence on the manner in which the field strength varied 
with antenna height. The section of foreground terrain 
under study in this case is neither randomly irregular 
nor is it a dominant obstacle, but takes the form of a 
relatively smooth gently varying departure from a 
spherical earth. The data were all taken on paths in Ari- 
zona over relatively smooth desert terrain. Some of the 
foreground terrain effects on these paths were briefly 
mentioned in previous papers.':? Ament and Katzin,?® 


* Manuscript received by the PGAP, November Pit NSIT 

} Smyth Research Associates, San Diego, Calif, 

*M. D. Rocco and J. B. Smyth, “Diffraction of radio waves 
es the earth,” Proc. IRE, vol. 37, pp. 1195-1203; October, 
1949, 

2J. P. Day and L. G. Trolese, “Propagation of short radio 
waves over desert terrain,” Proc. IRE, vol. 38, pp. 165-175; 
February, 1950. 

’W. S. Ament and M. Katzin, “Signal fluctuations in long-range 
over-water propagation,” IRE Trans. ON COMMUNICATIONS Sys- 
TEMS, vol. CS-4, pp. 118-122; March, 1956. 
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Fig. 1—Profile of Gila Bend to Dateland path. 


in microwave propagation studies in the Pacific, noted 
that foreground terrain features in the form of the mov- 
ing ocean swells caused variations in the field strength 
received on a shore-based receiving terminal. The fore- 
ground terrain was pictured as being “coupled” to the 
receiving antenna. 

In the present paper an attempt is made to relate a 
special shape of foreground terrain to a diffraction pat- 
tern noted in height-gain curves. 


EXPERIMENTAL SITE AND EXPERIMENTAL SETUP 


The data presented in this paper were obtained at an 
experimental station located in the Gila Valley in 
southern Arizona. This experimental station was orig- 
inally set up by the Navy Electronics Laboratory, 
San Diego, and a portion of the data presented was 
taken by the Navy Electronics Laboratory Staff. The 
experimental station is now under the cognizance of the 
Army Electronic Proving Ground, Fort Huachuca, 
Ariz. Three 200-foot towers on which elevator cabs 
were installed for hoisting antennas and associated 
equipment, are located at Gila Bend, Sentinel, and 
Dateland. Fig. 1 shows a profile of the East-West path 
from Gila Bend to Dateland, with the intermediate 
Sentinel tower indicated on the profile. Thus three paths 
could be set up: 1) a 46.3-mile path from Gila Bend to 
Dateland, 2) a 26.7-mile path from Gila Bend to Senti- 
nel, and 3) a 19.5-mile path from Sentinel to Dateland. 
On path 1 the terminals are beyond line of sight for all 
antenna heights available on the 200-foot towers. For 
paths 2 and 3, on the other hand, the antennas are 
within line of sight except for the lower 20 or 30 feet of 
elevator height excursion. Later a fourth link was set up 
by extending the path eastward to Casa Grande, Ariz. 
A mobile transmitter with antenna at fixed height was 
set up at Casa Grande and height-gain measurements 
taken at Sentinel. This 83.9-mile path is dominated by a 
large mountain obstacle, as may be seen from the profile 
of Fig. 13. 

For the tower terminals, antennas were installed on 
the face of the elevator cabs, and the receiving or trans- 
mitting equipment was installed inside the cab. Thus the 
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Fig. 2—Height-gain curves measured at Sentinel tower for 26.7-mile 
Vee from Gila Bend. Gila Bend terminal at fixed height of 190 
eet. 


field strength as a function of height could be recorded 
over a 190-foot height interval. 


EXPERIMENTAL DATA 
Gila Bend-Sentinel Path—26.7 Miles 


Some typical height-gain curves obtained at the 
Sentinel tower with the Gila Bend transmitter antennas 
at a fixed height are shown on Fig. 2 for frequencies of 
170, 520, 1000, and 3300 mc. An examination of these 
curves indicates that a diffraction pattern suggestive 
of knife-edge diffraction is superimposed on the height- 
gain shape normally expected at these frequencies. At 
3300 mc many small amplitude maxima and minima are 
present, with progressively closer spacing as the an- 
tenna height increases. For the other frequencies, the 
‘heights of the maxima and minima are more spread out 
in height as the frequency is lowered. 

For a knife edge in free space the height, h, of succes- 
sive maxima and minima of a knife-edge diffraction pat- 
tern is given by 


bd\n 
= / = (1) 


where 5 is the distance from the probing antenna to the 
knife edge; d, the distance between the probing antenna 
and the source;A, the wavelength; and 7 an integer which 
js odd for maxima and even for minima. If the diffraction 
phenomenon behaves geometrically like knife-edge dif- 
fraction, # should be proportional to the square root of X, 
and also to the square root of 7. On the other hand, if the 
maxima and minima are due to interference between a 
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direct component and one reflected from a smooth re- 


flecting ground, h would vary as the first power of \ and 
according to 


nn 
h=— 
Ay 


where y is the grazing angle of incidence of a plane wave 
on the ground. 

Fig. 3 shows a plot of #? as a function of \. For various 
values of the experimental points lie on straight lines 
passing through the origin in agreement with (1). Fig. 4 
illustrates the dependence of h# on n for 3300-mc data. 
The experimental points lie on a straight line with slope 
z on this log-log plot, again indicating agreement with 
(1). If the maxima and minima were due to interference 
between a direct and a ground reflected component, the 
plot of vs m should lie on a straight line of slope 1 rather 
than 3. Thus the geometrical behavior of the diffraction 
effect is like that of knife-edge diffraction rather than of 
specular reflection. 

Having noted the geometrical similarity of the data to 
knife-edge geometry, one can use (1) to compute the loca- 
tion of the apparent knife edge. However, since the refer- 
ence from which height should be measured is somewhat 
uncertain, it seems better to recast (1) in terms of the 
difference in height between minima or maxima. If this 
is done one obtains, after solving for 0, 


— d(h; — hi)? 
dN(/ 05 — Vmi)® + (lag — hes)? 


(2) 


(3) 


where h; and h; are the heights of any two maxima or 
minima of order 1; and n,. Substituting values from the 
experimental data, an average value of b was found to be 
0.47 mile. That is, the apparent knife edge should be 
located 0.47 mile east of the Sentinel tower. An examina- 
tion of the profile shows the terrain to have a rounded 
drop-off rather than a sharp obstruction like a knife edge 
at this location. 

In order to obtain a clue as to whether the terrain at 
this location was actually responsible for the knife-edge 
type of diffraction pattern, a mobile 3300-mc transmitter 
was set up at several points along the foreground terrain 
and height-gain curves recorded at the Sentinel tower. 
The results are shown in Fig. 5. Curve 4, taken with the 
transmitter 0.32 mile from the Sentinel tower, has max- 
ima and minima which are equally spaced in height; that 
is, h varies directly as m, as in (2), which described the 
geometry of an interference pattern resulting from a di- 
rect and ground reflected component. Curves B, C, and 
D, on the other hand, taken with the transmitter located 
0.68 mile or more from the Sentinel tower, have a dif- 
ferent character. The spacing of maxima and minima 


4 For the plots of both Figs. 3 and 4 a constant 17 feet was added 
to the measured heights in order to obtain better fit to straight lines. 
As long as the same constant is added to all heights, this does not 
seem unreasonable since the reference from which height should be 
measured and the angle of arrival of the wave is not known precisely. 
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Fig. 3—Variation of square of height of maxima and minima of 
Sentinel height gain curves as function of wavelength. n=1 for 
first maximum, 2 for first minimum, etc. 
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Fig. 4—Variation of height of Sentinel 3300-mc maxima and 
minima with the diffraction order number. 


decreases with height more like knife-edge diffraction, 
with # proportional to the square root of n. Clearly 
the terrain between 0.32 and 0.68 mile is responsible for a 
diffraction effect which has geometrical characteristics 
like that of knife-edge diffraction even though the phys- 
ical appearance of the terrain is smooth and rounded 
rather than sharp-edged. 

Another facet of the effects of terrain on radio fields is 
the intensification of the knife-edge type diffraction pat- 
tern when meteorological conditions result in ducting 
phenomena. In this desert area a diurnal meteorological 
cycle is usually evident due to solar heating of the ground 
during day and radiational cooling of the ground at 
night.? This results in a nocturnal ground-based tem- 
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Fig. 5—Height-gain curves measured at Sentinel tower using the 


3300-mc mobile transmitter at four locations along the foreground 
terrain east of Sentinel. 
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Fig. 6—Nighttime 3300-mc Sentinel height-gain curves 
for Gila Bend-Sentinel path. 


perature inversion which reaches its maximum intensity 
during the early morning hours. Occasionally during the 
early morning hours a striking intensification of the dif- 
fraction pattern was noted, as illustrated by the height- 
gain curves of Fig. 6. Although the minima become much 
deeper, the pattern appears to be due to the same dif- 
fraction effect noted in the daytime. This is demon- 
strated by Fig. 4, which shows # plotted as a function of 
m. On this plot the crosses are daytime data and the 
circles, nighttime points. Both sets of data lie on straight 
lines with slope of 4, but the nighttime diffraction pat- 
tern is shifted downward slightly. No doubt this is due 
to the fact that the additional refractive bending of the 
wave changes the angle of arrival slightly, thus moving 
the whole diffraction pattern slightly downward. 

The deepening of the minima are believed to be at- 
tributable to a second wave component which is re- 
flected from the section of terrain (Fig. 1) between 10 and 
17 miles from Gila Bend. This component interferes with 
the direct wave only under conditions of strong surface 
ducting. At heights where destructive interference be- 
tween the two waves reduces the field, the knife-edge 
oscillations become a larger fraction of the total field, 
thus increasing the depth of the knife-edge minima. 
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Fig. 7—Daytime Gila Bend and Dateland 3300-mc height-gain 
curves for Gila Bend-Dateland 46.3-mile path. 


Gila Bend-Daieland Path—46.3 Miles 


= y 


Typical daytime height-gain curves obtained for 3300- 
mc transmissions over the 46.3-mile Gila Bend-Dateland 
path are shown in Fig. 7. The dashed curve is a height- 
gain curve taken at Dateland with the Gila Bend ter- 


- minal at a fixed height of 194 feet. The dotted curve is a 
_ Gila Bend height-gain curve with the Dateland terminal 


at a fixed height of 194 feet. The solid line is a theoretical 
curve computed by use of standard diffraction theory 


_ for a smooth spherical earth. Since the atmosphere is es- 


sentially standard in the daytime in Arizona, both meas- 
ured height-gain curves should agree with the computed 


curve if the path were truly spherical. The Gila Bend 
_ curve is. quite similar to the theoretical curve, but the 
Dateland height-gain shows a bulging maximum near 


100 feet and minimum at 140 feet. During the nighttime, 
when the nocturnal duct forms, the height-gain curve at 
Dateland varies considerably, as may be seen for the 


series of height-gain curves shown in Fig. 8. During the 


night and early morning hours, when the nocturnal duct 
is most intense, the maximum and minimum visible on 
the daytime height-gain curves drop in height and a sec- 
ond maximum and minimum are present in the 190-foot 
height interval. After sunrise the pattern returns to the 
normal daytime condition indicated by the curve for 
1258. This effect is shown more clearly in Fig. 9, which 
is a plot of the height of the minima as a function of time 
of day for a 24-hour cycle. The height of the first mini- 
mum is fairly constant during the day and drops during 
the night, reaching a low value of 100 feet in the early 
morning hours. The second minimum drops in height 


during the night‘in synchronism with the first minimum. 


In this case the pattern is more pees out in height than 


Foreground Terrain Effects on UHF 


333 


150 TIME 2uz, 2248. 2321 0244 


50 40 30 40 30 20 50 40 30 50 ure) 30 


ANTENNA HEIGHT IN FEET 


—= ss — 
30 20 10 40 Ol 20 50 40 30 60 50 40 
FIELD STRENGTH IN DB RELATIVE TO ONE MILLIWATT 


Fig. 8—Dateland 3300-mc height-gain curves for 
Gila Bend-Dateland 46.3-mile path. 
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Fig. 9—Variation of height of Dateland diffraction 
minima with time of day. 


was the case for the 3300-mc pattern at Sentinel. Thus 
if the phenomenon is similar to that at Sentinel, the ap- 
parent knife-edge location must be more distant. With 
the apparent knife edge located at greater distance, 
changes in angle of arrival due to changing refraction 
would shift the diffraction pattern up and down to a 
greater degree than was noted at Sentinel. Since a smaller 
portion of the diffraction pattern appears on the 190- 
foot height interval which is probed, it is more difficult 
to determine if a similarity to a knife-edge diffraction 
pattern exists. However, for the nighttime data, when 
more maxima and minima are present, a plot of h vs n 
can be made, as shown on Fig. 10. Points taken from the 
height-gain curves on Fig. 8, at time 0244, 0445, and 
0805, are shown in Fig. 10. Each set of points lies quite 
close to a line with slope 4 on the log-log plot, indicating 
a dependence of / on n similar to knife-edge diffraction. 

If (3) is again used for values from the Dateland ex- 
perimental data, an apparent knife-edge location 4.7 
miles from the Dateland tower is obtained. In this cal- 
culation d was taken as the distance from the Dateland 
tower to the hill just west of Sentinel (see profile, Fig. 1) 
which is at the horizon. 
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Fig. 10—Variation of height of Dateland 3300-mc maxima 
and minima with the diffraction order number. 


Thus, just as in the East Sentinel foreground effect, the 
diffraction phenomenon is geometrically like knife-edge 
diffraction. However, the amplitude of the signal varia- 
tions with height are much greater than knife-edge theory 
allows. 


Sentinel-Dateland Path—19.5 Miles 


An examination of the profile of Fig. 1 shows that a 
rounded section of terrain is present west of the Sen- 
tinel tower. If the diffraction patterns described previ- 
ously are attributable to the rounded shape of the fore- 
ground terrain, one would expect a similar result from 
the rounded section west of Sentinel. To check this 
point, a 1725-mc link was set up on the Sentinel-Date- 
land path and height-gain curves were obtained at 
Sentinel for different fixed terminal heights at Date- 
land. Fig. 11 shows daytime height-gain curves at 
Sentinel with the Dateland terminal at different fixed 
heights. Most of the curves show two maxima and two 
minima which stay about fixed in height as the Date- 
land terminal height is set at different values. Only for 
the Dateland terminal at the greatest height, 197 feet, 
is this pattern not clearly present. For nighttime data, 
Fig. 12, the situation is similar except that the pattern 
is shifted downward in height slightly and the minima 
are deeper. Thus most of the height-gain curves appear 
to show a diffraction pattern similar to those previously 
described. If the location of the apparent knife edge is 
computed, as was done before using (3), a location 1.7 
miles west of the Sentinel tower is obtained as a median 
value. This location is found on the rounded bulge in a 
position comparable to that in the other two cases. 


Casa Grande-Sentinel Path—83.9 Miles 


A profile of this path is shown on Fig. 13. Propagation 
studies over this path were made for about a one-week 
period to ascertain whether the foreground terrain ef- 
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Fig. 11—Daytime 1725-mcs Sentinel height-gain curves 
for Sentinel-Dateland path. 
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Fig. 12—Nighttime 1725-mc Sentinel height-gain curves 
for Sentinel-Dateland path. 


fects noted at Sentinel would be present on this 83.9- 
mile path as well as on the much shorter Gila Bend- 
Sentinel path. As may be seen from the profile, this 
path is dominated by a large mountain obstacle, the 
Maricopa Mountains. The arrival angle of the wave 
reaching Sentinel should be a little greater than for the 
short path. The Casa Grande antenna was located at a 
fixed height of 10 feet and height-gain measurements 
were made at Sentinel. Over this path the general level 
of the signal seemed to be dominated most of the time 
by diffraction over the mountain obstacle rather than 
by tropospheric scattering. Occasioually rapid fluctua- 
tions characteristic of a scattered field signal were pres- 
ent, but most of the time they were absent. Typical 
height-gain curves are shown in Fig. 14. The maxima 
and minima of a characteristic knife-edge type diffrac- 
tion pattern are clearly visible on the curves. The plot 
of h as a function of m, shown on Fig. 15, again indicates 
that the diffraction pattern obtained is geometrically 
like knife-edge diffraction but the amplitude variations 
are much larger than knife-edge theory predicts. A cal- 
culation of the location of the apparent knife edge would 
place the knife edge 0.39 mile east of the Sentinel tower, 
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Fig. 13—Profile of Casa Grande to Sentinel path. 
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Fig. 14—Sentinel 1800-mc height-gain curves for Casa Grande 
to Sentinel 83.9-mile path. 


| compared to a value of 0.47 mile obtained from the 


previous data on the Gila Bend-Sentinel path. How- 


ever, on the Casa Grande-Sentinel path, the energy ar- 


rives at a little higher angle due to the presence of the 
mountain obstacle, and one might expect the apparent 
knife edge to appear to be closer to the tower. This point 


will be discussed in more detail later. 


General Summary of Experimental Results 


The results from tests on the four paths described 


_ show diffraction effects which can be attributed to the 


shape of the foreground terrain on each of the paths. 


In order to show the similarities in the three sections 


of foreground, profiles plotted to expanded scales are 


; 


indicated on Fig. 16. The Dateland foreground profile is 


at the top of the figure, the foreground west of Sentinel 
in the center, and that east of Sentinel at the bottom. 


The West-Sentinel profile distance scale is expanded 
three times the Dateland scale. The East-Sentinel pro- 


file is plotted on a distance scale expanded tenfold, 
relative to the Dateland plot. The arrows indicate the 
location computed for the apparent knife edge and in 
each case this turns out to be in the same general posi- 
tion on the rounded foreground profiles. It is logical to 
conclude that the three diffraction effects are phenom- 
enologically alike. The similarity may be noted in 
another way by examining the Dateland 3300-mc 
height-gain curve of Fig. 7 and the 170-mc Sentinel 
curve of Fig. 2. The two curves are quite alike in general 
character, indicating that the diffraction effect is simi- 
lar in each case except for geometric scale. 
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Fig. 16—Expanded scale plot of Dateland, Sentinel West, 
and Sentinel East foreground terrain profiles. 


Since the geometry of the diffraction phenomenon 
from the rounded foreground apparently coincides with 
knife-edge geometry, it would be convenient to have an 
empirical method of estimating the location of the ap- 
parent knife edge. If this location is known, the height 
of the first maximum could be computed for a given 
wavelength to insure siting an antenna at optimum 
height. 

In developing a method for predicting the location of 
the knife edge, it was felt desirable to express the dis- 
tance from terminal to knife edge, b, in terms of param- 
eters describing the shape of the terrain bulge and the 
angle of incidence of the energy on the terrain. Fig. 17 
illustrates the concept chosen for this purpose. A cylin- 
der of radius R is fitted to the actual foreground terrain 
profile. Radii are drawn from the center of curvature of 
the cylinder to 1) the terminal location, 2) the knife- 
edge location, and 3) normal to the incoming rays. The 
angles a and @ are defined as shown. Then the angular 
distance from terminal to knife edge }/R=a+6. For 
the terrain features studied thus far, the measured pa- 
rameters 6, R, and @ are indicated in Table ie 
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TABLE I 
: R b b/R (1/R 
Terrain Feature oe : Feet (mr) a(mr) 6(mr) ee Lo 
Dateland 4.7 25,000 5.3 2.5 2.8 0.174 
Sentinel West 0.49 9000 18.5 9.8 8.7 2.02 
Sentinel East 
(to Gila Bend) 0.15 2500 16.7 LD 15.5 6.63 
(to Casa Grande)| 0.15 2100 13u7 —-1.4 15.1 6.63 


From these parameters, 6 has been calculated for each 
feature, together with (1/R—1/a), (where a is the effec- 
tive earth’s radius). The effect of increasing ray inclina- 
tion causing a closer knife-edge location is illustrated in 
the bottom two lines of the table. Although the angular 
distance to the knife edge is reduced 3 mr (16.7 —13.7) 
in the mountain obstacle path, the increased ray inclina- 
tion of 2.6 mr almost completely compensates this effect, 
such that @ remains essentially constant. If one now 
plots 6 vs (1/R—1/a), a straight line relationship is 
found on the log-log plot of Fig. 18. The parameter 
(1/R—1/a), rather than R itself, was chosen to desig- 
nate the shape of the bulge, since this parameter be- 
comes zero as the bulge flattens down to the smooth 
spherical earth. Since @ should also become zero in this 
case, the log-log plot would be expected to be linear over 
a wider range of bulge radii than if R itself is used. 

Using this empirical curve, one can now predict the 
location of a knife edge on any terrain bulge of radius R 
and ray inclination a. For example, if the bulge radius 
R is 1 million feet, (1/R—1/a) would be (1—0.039) 
= 0.96 X10~*/feet. According to Fig. 15, 6 is 6.4 mr. If 
the ray inclination is 2 mr, then b/R=8.4 mr, or 
6 = 8400 feet. It remains to be seen whether this model 
is sufficiently general to cover most foreground terrain 
bulges. Work on similar features is continuing for the 
purpose of checking this point. 

In order to obtain some indication of the practical 
importance of the foreground diffraction effects which 
have been described, the 3300-mc data available on the 
46.3-mile Gila Bend-Dateland path were examined 
statistically. The data for this analysis were obtained 
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Fig. 18—Experimental relationship between (1/R—1/a) and 6. 


from a series of height-gain curves similar to those 
shown on Fig. 7. A height-gain measurement was made 
by recording the field strength as the Gila Bend termi- 
nal was moved in height with the Dateland terminal at 
a fixed height of 194 feet. Immediately thereafter the 
Gila Bend terminal was fixed at 194 feet and the Date- 
land terminal moved in height. Most of the Gila Bend 
height-gain curves agreed with the computed smooth > 
earth curve and exhibited no foreground diffraction 
pattern. The Dateland curves, on the other hand, 
showed a consistent maximum in field strength at about 
70 feet. The histogram of Fig. 19 was obtained by tak- 
ing the ratio of the field strength received on the follow- 
ing two paths: 1) the Gila Bend terminal at 190 feet 
and Dateland at 70 feet; 2) the Gila Bend terminal at > 
70 feet and Dateland at 190 feet. Three histograms are 
shown on Fig. 19, one for daytime data with no duct 
present, a second with nighttime data when a duct was 
present, and a third with all of the data combined to- 
gether. For a spherical earth the histograms should, of 
course, peak up at 0 db. However, the plots show that 
the field on the Gila Bend 190-Dateland 70 path con- 
sistently exceeds the field on the Gila Bend 70-Date- 
land 190 combination. For nonducting conditions the 
bias in favor of the Gila Bend 190-Dateland 70 is most 
pronounced, averaging close to 15 db. 

If one terminal is set at a very low height of 5 feet 
instead of 190 feet, the results shown on Fig. 20 are ob- 
tained for the other terminal at either 50 or 100 feet. 
This shows the effect of the Dateland foreground diffrac- 
tion pattern. By properly siting the antenna on the first 
maximum a diffraction enhancement is obtained due to 
the diffraction effect of the rounded foreground terrain. 

On the 83.9-mile Casa Grande-Sentinel path the 
1800-me data serve to illustrate the effect of the fore- 
ground terrain east of the Sentinel tower for a link hav- 
ing rather low terminal heights. Fig. 21 shows a sum- 
mary of the ratio of the fields for an antenna height at 
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:, Fig. 20—(a) Histogram of ratio in db of field received on Gila Bend 
5 feet-Dateland 100 feet path to that received on Gila Bend 100 
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ceived on Gila Bend 5 feet-Dateland 50 feet path to that received 
on Gila Bend 50 feet-Dateland 5 feet path. 


first maximum as compared to that at the first mini- 
_ mum of the foreground diffraction pattern. This is given 
“in the form of a distribution curve which shows, for 
example, that 50 per cent of the time the signal for an 
antenna height on the maximum (19 feet) is 9 db higher 
than on the minimum (33 feet) and 10 per cent of the 
time it is 17 db higher. Furthermore the fading was more 
severe with the antenna at the height of the minimum 
than of the maximum. This point is illustrated in Fig. 
22, which shows a sample of the signal record for an- 
tenna heights of 6 feet (below first maximum), 21 feet 
(near first maximum) and 31 feet (near first minimum). 
At the minimum the fading range is about 15 db and 
on the maximum less than 5 db. On this particular day 
the received signal appeared to be influenced by 3 dif- 
ferent mechanisms: 1) diffraction over the dominant 
“mountain obstacle, 2) scattering in the lower tropo- 
sphere, and 3) diffraction due to the foreground terrain 
near Sentinel. Most of the time the influence of the 
scattered field component was less evident. 


CONCLUSIONS 


The experimental results which have been summa- 
rized indicate that gently rounded sections of foreground 
‘terrain can cause marked diffraction effects at UHF 
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feet-Dateland 5 feet path. (b) Histogram of ratio in db of field re- 
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Fig. 21—Distribution curve of ratio of field at first maximum to that 
at first minimum. Casa Grande-Sentinel path—1800 mc. 


Fig. 22—Sample field strength record at Sentinel showing time 
fluctuations of 1800-mc signal with the antenna below the first 
maximum (6 feet), on the first maximum (21 feet), and on the 
first minimum (31 feet). 


frequencies. The diffraction phenomena show a striking 
geometrical similarity to knife-edge diffraction. How- 
ever, the amplitude of the spatial variations of field 
strength are considerably greater than for knife-edge 
diffraction. 

Proper siting of antennas overlooking such fore- 
ground terrain features can lead to sizable diffraction 
enhancement of the received signal. 

An empirical method for predicting the foreground 
diffraction geometry from foreground terrain profiles 
having the rounded shape has been devised. If this 
method can be verified experimentally for other cases, 
it should be useful as an aid in proper siting of terminals 
of UHF links. 

Changing refraction due to meteorological variations 
can alter both the position in height and intensity of the 
diffraction pattern. The height variation due to re- 
fraction changes becomes less as the distance of the ap- 
parent knife edge from the terminal becomes smaller. 
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A Rapid Beam-Swinging Experiment 
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in Transhorizon Propagation* 


ALAN T. WATERMAN, JR.t 


Summary—By using a broadside phased array for an antenna, a 
narrow beam can be swung rapidly and in quick succession through 
a limited sector by fast control of the phasing, rather than by move- 
ment of the entire antenna structure. This technique is used at the 
receiving end of a 101-mile beyond-the-horizon transmission path 
in order to probe the portion of the troposphere through which the 
signal is propagated. At the frequency employed of 3.12 kmc, a 0.49- 
degree beam is swung in azimuth through a 4.2-degree sector each 
tenth of a second. 

oA variety of phenomena are observed with this technique which 
have not been directly apparent in slower beam-swinging experi- 

. ments. The beam-broadening effect attributed to atmospheric scat- 
tering is not always evident on any one sector scan. However, the 
change from scan to scan is frequently rapid enough so that a time 
average would show the broadening. At times the scan-to-scan 
changes are systematic and show a continuity indicative of a motion 
of the scattering or reflecting regions; in some cases this motion is 
too rapid to be accounted for by transport of air, thus implying a wave 
motion rippling through the atmosphere. At other times the atmo- 
spheric structure is too fine to be resolved by the beamwidth em- 
ployed, and the time variations are too rapid to reveal a continuity 
from one scan to the next. 


effectively as devices for probing the atmosphere 

in transhorizon propagation experiments.!~7 

The unmanageable nature of the large physical struc- 
ture required for a narrow beam is largely circumvented 
if a phased array is used and a mechanism provided for 
rapid variation of the phase of each radiating element. 
The experiment described here utilizes such an array 
for the receiving antenna. A 3000-mce circularly polar- 
ized signal is transmitted, with sufficiently wide beam- 


IN cectivety a antennas have been used very 


* Manuscript received by the PGAP, June 30, 1958. This paper 
was presented at the URSI meeting, Washington, D. C.; April 26, 
1958. The research reported was made possible through the support 
of the U.S, Army Signal Corps, under Contract DA36(039)SC-73151. 
. A Stanford Electronics Labs., Stanford University, Stanford, 

alif. 
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width (about 5 degrees) to illuminate the pertinent por- 
tion of the atmosphere. The receiving antenna® is an_ 


eight-element broadside array (Fig. 1). Each element 
in the array consists of an offset parabolic reflector and 


helical feed. The helix in each feed is connected to a | 


syncrodrive so that the desired phase shifting can be 


accomplished by mechanical rotation of the helix about 


its axis. The syncros then are driven at multiple rates 
from a centrally located gear train. Thus, the over-all 
rate of beam swinging is controlled by the speed of the 
master motor driving the gear train. 


Fig. 1—Phased array used as the receiving antenna in 
rapid beam-swinging experiment. 


The horizontal beamwidth of the entire antenna de- 
pends on the lateral length of the over-all array. At 
3000 mc, a 30-foot-long array has a beamwidth of 
about } degree. In the array which has been constructed 
each element (consisting of reflector and feed) is mounted 
along a pipe supported on an A-frame structure, so 
that the element can be slid laterally along the pipe, 
thus affording an adjustment of the length of the array 
and consequently of the width of the antenna pattern. 
The interval between successive principal maxima in 
the array pattern depends on the spacing between ad- 
jacent elements in the array or, for an array of fixed 
length, on the number of elements. This angular inter- 
val is also the sector over which the array can scan with- 
out ambiguity. The sector of interest in these trans- 
horizon propagation experiments covers only 4 or 5 de- 
grees in azimuth. Consequently only eight elements are 
needed in the array and these may then be placed ata 
4-foot spacing. Utilization of a parabolic reflector in 
each element does not affect the array factor over this 


*R. E, Miller, G, K. Durfey, and W. H. Huntley, Jr., “A Rapid- 
Scanning Phased Array for Propagation Measurements,” Stanford 
Electronics Labs., Stanford, Calif., Tech. Rep. No. 461-5; July 30 
1958, under Signal Corps Contract’ DA36(039)SC-73151. 
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limited sector but does increase the capture area of the 
over-all antenna. It also establishes the vertical beam- 
width, in this case about 5 degrees. 

For operation in this frequency range, all the com- 
ponent parts are of convenient physical size. The 
helices® consist of 3.4 turns of No. 12 copper wire and 
are 7.8 cm long. The half-power beamwidth for each 
helical feed is 60 degrees, a convenient value for illu- 
minating the offset section of paraboloidal reflector. This 
offset places the feed itself and associated syncrodrive 
out of the path of the reflected radiation and in a con- 
veniently accessible position. 

To provide a continuous azimuthal scan, the helical 
feeds are rotated at multiple rates—or rather, the right 


four helices are rotated at odd multiple rates in one 
sense and the left four at corresponding rates in the 


opposite sense. With this procedure, a rotation rate of 
2100 rpm for the outermost helices corresponds to a 
scan rate in which one principal maximum moves on 
and the next moves into its place in one-tenth second 


_ (t.e., ten sector scans per second). It is to be noted that 


this beam swinging is not a back-and-forth sector scan; 
the principal maxima in the array factor rotate uni- 
formly like the spokes of a wheel and move one after 
another through the sector being scanned. 

At the frequency employed, 3.12 kmc, the horizontal 
beamwidth of the array and the azimuthal sector 
scanned, as measured on a line-of-sight test range, are 
0.49 degree and 4.2 degrees, respectively. 

This equipment has been employed in propagation 
measurements which, though preliminary, are never- 
theless startling enough in their implications to warrant 
reporting promptly. A CW transmitter was located 101 
miles distant from the receiving array, well beyond line 
of sight, and the system response was observed as the 
receiving antenna beam scanned in azimuth through 
the great-circle bearing to the transmitter. 

In the accompanying figures, which are pictures of 


an oscilloscope face, signal level is displayed as a vertical 


deflection and the azimuth in which the beam is point- 
ing is synchronized with the horizontal sweep. Fig. 2 


' shows the response of this antenna system to a low- 


power test transmitter within the line of sight of the 
receiving location. It serves to illustrate the beamwidth 
(0.49 degree) and sidelobe level (better than 10 db 
down) of the array. It also provides an azimuth refer- 
ence for the sector scan. The ten divisions of horizontal 
displacement in the figures correspond to the 4.2-degree 
sector. On this scale, the great-circle bearing to the 101- 
mile distant transmitter (not to the test source in Fig. 
2) is 7.2 divisions from the left-hand edge. 

Fig. 3 is an example of receiving-system response to a 
signal propagated over the 101-mile path. There the 
display is the same, except that a raster containing 
alternate antenna sector-scans is employed. Thus, each 


9W. H. Huntley, Jr., “A Helical Microwave Antenna Feed,” 
Stanford Electronics Labs., Stanford, Calif., Tech. Rep. No. 461-3; 


September 20, 1957, under Signal Corps Contract DA36(039)SC- 


A319), 


Waterman: A Rapid Beam-Swinging Experiment in Transhorizon Propagation 


339 


Fig. 2—Measured antenna pattern of the array 
over a 4,2-degree sector. 
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Fig. 3—Raster display of receiving-system response to 101-mile- 
distant transmitter. Each azimuth scan covers 4.2 degrees (10 
divisions in the figures) in 0.1 second; successive scans displayed 
are 0.2 second apart. Great-circle bearing to transmitter is 7.2 
divisions from left. 


scan consumes one tenth of a second, but there is one 
fifth of a second between successive scans in the raster. 
It will be noticed that the effective antenna pattern 
obtained via the atmospherically scattered signal varies 
rapidly from scan to scan. The beam-broadening effect 
attributed to random scattering! *® is not very marked 
on any one scan, but if a succession of scans is averaged, 
the shift in angle of arrival plus the apparently discrete 
off-path contributions (which show up here looking like 
exaggerated sidelobes) result in considerable beam 
broadening. 

An example in which there appears to be a small 
number of discrete scatterers which change quite rap- 
idly is shown in Fig. 4. A slower beam-swinging experi- 
ment would not reveal this fine structure; only an 
average broadening would be observed. 

The predominant characteristic evident in Fig. 5 is 
the rapid and continuous motion of one discrete scat- 
terer. Such a motion, of course, does not necessarily 
imply a physical translation at this speed (which would 
correspond roughly to 600 or 800 mph), but might be 
caused by a wave motion in a reflecting layer. 

By way of contrast, Fig. 6 shows a constant pattern 
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Fig. 4—Sample of received signal. (See Fig. 3 for interpretation.) 


Fig. 5—Sample of received signal. (See Fig. 3 for interpretation.) 


with very little beam broadening. This constancy, when 
observed, occurred at night, in contrast to the more 
wildly fluctuating patterns seen during daytime. 

A final example has been selected for Fig. 7 in which 
several of the sweeps show an instantaneous beam 
broadening. This type of pattern, which so far has not 
been observed frequently, might be what one would ex- 
pect if there were a large number of small incoherent 
scatterers. 

These figures are samples of observations which have 
been made recently and have been repeated sufficiently 
to establish their reality. A detailed analysis of the 
results has not yet been undertaken. Nevertheless, 
some pertinent remarks can be made. 

There are times when the antenna beam, though 
broad in elevation and only moderately narrow in azi- 
muth, appears to have sufficient resolution to distin- 
guish discrete scatterers. (The term “scatterer” is used 
here in its broadest sense; it might well refer to an ap- 
propriately oriented portion of a partially reflecting 
layer.) These scatterers, in such cases, are small in num- 
ber (one to four); their resolvability implies a spatial 
separation of at least 500 to 1000 feet. At other times, 
the scatterers are not resolved, and few conclusions can 
be drawn at the moment. 
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Fig. 6—Sample of received signal. (See Fig. 3 for interpretation.) 
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Fig. 7—Sample of received signal. (See Fig. 3 for interpretation.) 


Frequently, particularly in the daytime, the at- 
mosphere changes structure in one fifth of a second to 
such an extent that little continuity can be found in the 
changes. In such cases one cannot distinguish a random 
type of turbulent process from a systematic wave mo- 
tion. At other times, the one-fifth second interval be- 
tween scans is sufficiently short to establish a continuity 
of change from one situation to another. In these cases, 
there are many instances in which the apparent motions 
are more rapid than can reasonably be accounted for in 
terms of tropospheric winds; a wave motion of some 
sort is implied, such as ripples along an interface. 

Much of the beam broadening does not occur instan- 
taneously but shows up in a time average owing to the 
rapid changes in the scattering or reflecting centers. 
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Effect of Mountains with Smooth Crests 
on Wave Propagation* 
I. P. SHKAROFSKY}, H. E. J. NEUGEBAUER}, anv M. P. BACHYNSKIt 


Summary—A new method of solving the problem of diffraction 
of EM waves by the smooth crest of a perfectly reflecting cylindrical 
mountain has been previously reported.! This paper presents the 
results in a form more suitable for practical applications. The theory 
is extended, and good agreement with model experiments is ob- 


_ tained for scattering angles up to 5 degrees. The procedure for 

including the effects of reflections from the ground on either side of 
_ the mountain is also indicated. A few examples illustrate cases en- 
_ countered in practice, and exhibit effects up to 8 db compared with 
_knife-edge diffraction. 


Il. INTRODUCTION 


N a previous paper,! the problem of diffraction of 
EM waves by a mountain with a smooth crest of 
approximately cylindrical shape has been solved 


and has been shown to give good agreement with meas- 


urements for small scattering angles. In that theory, 


- additional halo terms were added to the usual knife-edge 


term to account for the effect of the smooth crest of the 


mountain. The appropriate polarization effects were 
obtained using coefficients previously determined by 
- Fock.? The power behind a mountain agreed well with 
model experiments at K-band frequencies for a wide 
range of distances, radii of curvature, and a narrow 


_ range of angles of diffraction. 


This paper presents numerical tables and graphs 


which are useful for quick calculation. The functions 


that determine the real and imaginary parts of the halo 
terms are tabulated. The real and imaginary parts of 
the halo and knife-edge terms are also plotted in graphs 


that serve to simplify calculations. These real and imag- 
_ inary parts, rather than the absolute values, are required 
for calculating the additional effect of reflecting ground 


. 


on either side of the mountain. This is accomplished 
in the usual manner*® by assuming that the reflected 
waves from the ground originate from virtual images. 
The field is thus propagated in four ways and the sums 
of the real and imaginary parts of these four contribu- 
tions are required to calculate the received power. 
When the ground has no effect on the propagation 


and when diffraction occurs only from a smooth crest of 


a mountain, the real and imaginary parts of the knife- 
edge and halo terms can be combined to express the 


* Manuscript received by the PGAP, February 25, 1958. 

+ Res. Lab., RCA Victor Co., Ltd., Montreal, Can. | : 

1H. E. J. Neugebauer and M. P. Bachynski, “Diffraction by 
smooth cylindrical mountains,” Proc. IRE, vol. 46, pp. 1619-1627; 
September, 1958. Some of the material in this paper was presented 
at the General Assembly of URSI in Boulder, Colo., August /Septem- 

r, 1957. 

3 2V.F ock, “The distribution of currents induced by a plane wave 
on the surface of a conductor,” J. Phys. (USSR), vol. 10, pp. 130- 


: February, 1946. 
‘eas 5 eo Schalline, C. R. Burrows, and E. B. Ferrell, “Ultra-short 


“wave propagation,” Proc. IRE, vol. 21, pp. 427-463, March, 1933. 


power as a function of two universal parameters, X and 
Z. These variables‘ are defined by 


X = Wo(2ha)'?,  Z = Wo(2kd/a)}!? (1) 


where 2y is the angle of diffraction, k=27/\, \=wave- 
length, @ is the radius of curvature at the crest of the 
mountain, and d is the harmonic mean between receiver 
and transmitter distances from the mountain. This will 
be defined explicitly later. 

Comparisons of theory and experiment were given only 
for small diffraction angles up to 1.5 degrees.! It is now 
shown that experimental results agree up to angles of 
5 degrees and in some cases 10 degrees, depending on 
the radius of curvature. Although the largest experi- 
mental value of ka is only 250, there is reason to believe 
that even better agreement is obtainable for larger ka. 
The theory is assumed to hold for all ka. The theoret- 
ical graphs of power vs scattering angle are represented 
in a universal manner by plotting power vs X for vari- 
ous (X/Z). 

Finally, some numerical examples are given showing 
the ratio of received power behind a mountain with a 
smooth crest as compared to a knife-edge mountain. 


II. Basis or THEORY 


The electric field diffracted by a smooth cylindrical 
mountain is given by (2) to (5), which are derived from 
(23) and (27) of the previous paper.’ The field can be 
written as the sum of the real and imaginary components 
(Kye, Kim) of the knife-edge contribution and the com- 
ponents (H,e, Him) of the halo contribution. These four 
components are functions of the two universal param- 
eters X and Z. 


E=[1/2(ditd2)| exp [—ja(di tds) (1+Y0?/2)] 
% | [Kr(Z)+IK im(Z)|+ (X/Z) [Hre(X) +7 im(X)]} 
where 


Ke(Z) Sie 41K im(Z) 
= {1— (1+ p[C(Z) — 9S(Z)]} exp (jrZ?/2) (3) 


C(Z) and S(Z) are the Fresnel integrals 


(2) 


C(Z) = feos (ré?/2)dé, S(Z) -{ sin (wé?/2)dé. (3a) 


4In the knife-edge diffraction theory of Fresnel-Kirchhoff, the 
Fresnel cosine and sine integrals have as a limit the parameter 
y=2/kd2/ sin Wo where d, is the distance from the knife edge. If 
the transmitter to the knife-edge distance approaches infinity, d in 
(1) becomes equal to 2d. For this limit and for small yo, the param- 
eter v can be identified with the symbol Z used here. 
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Fig. 1—Geometry for a mountain with a smooth crest. 


Thus 
K,(Z) = [1 — C(Z) — S(Z)] cos (4Z?/2) 


— [S(Z) — C(Z) | sin (Z?/2) (3b) 
Kim(Z) = [1 — C(Z) — S(Z)] sin (4Z?/2) 
+ [S(Z) — C(Z)] cos (#Z?/2) (3c) 


where by (1) 


Z = Wo(2kd/x)!2, = X = Wo(2ka)/8, 


Also 
HX) + 7H im(X) = (29/?/m)e™/4(2ka)?/8p 


{ ib “WOv + ho) |"? exp [—jha(Spop? + 29%) dy 


+ (9/2) { ~ yexp [~2jk0 (Tu4/2 + 2yr)]av! (4a) 


or 


H,«(X) + jHin(X) 
- (28/n)o4 if “(€M)"?(cos P — j sin P) 
-exp [—5/3X#2/4 — #8] dé 
+ (0/28!) f "e(cos Nj sin N) 
-exp [—7./3X'#?/8 — elas} ,_ (4b) 


where 


p is the polarization parameter, 


p = 0.7 for vertical polarization, 


p = — 1.0 for horizontal polarization. (5a) 
M? = (3 + X/3)? + x2, 
P= n/12 — p/2 + 5X#/4 (5b) 
tany = X/(3§ + Xv/3) 
X’ = x28, 
N = w/12 + 7X'€/8. (5c) 


Fig. 1 defines the remaining symbols in the equations. 
The center plane NM through the crest of the mountain 
bisects the angle 7;NT7»2 where T,N and T>N are the 
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pepe teal part, K;., vs parameter Z =yo(2kd/) 1/2. For com- 


RD RO acs em ne pope RR erat 


parison, Px=|K|*/4(di+-d2)* represents total power behind a_ 


knife-edge mountain. 


LO 1.2 


—s» 7 


0 02 04-06 O8 


Fig. 3—Knife-edge contribution to the diffraction field behind a 
mountain: imaginary part, Kim, vs parameter Z =yo(2kd/r) 1/2, 


tangents to the mountain from the transmitter 7, and 
receiver 7». The vertical plane is SM. ° 


Let 


di, dz=11Wi, T,W2=respective transmitter and re- 
ceiver distances measured to the plane VM. 
D,, D2=T1Vi, T2V2=distances from T, and T: to the 
plane SM, 
d= 2did,/(di+d2) is the harmonic mean between 
distances d; and dy, 
MN, H,=SVi1, SV2=distances below top of mountain 
for 7; and TJ». 
2fo=diffraction angle=180°~— / Ti, NT». 
Yo= ZNT{Wi= ZNT2W,. 
2) =W1+WV, where VW, = Z NT,V; and 
W2= ZL NT2Vs. 


Also 


That is, VW, and W, are the angles subtended at T; and 


I>, respectively, by the mountain. For small angles, 


Wield, Vo =~ Hoy Ds, 2f0 ~ Hi/Di + Ho/ Dz. 
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Fig. 4—(a) Knife-edge contribution to the diffraction field be- 

“hind a mountain: curves represent (tZKy.), (—1ZKim), and 
[(rZ)?|K|?/2] against Z and are more accurate than those of 
Figs. 2 and 3 for large Z. (b) Knife-edge contribution to the 
diffraction field behind a mountain: same functions as in Fig. 4(a) 
for large values of Z. The functions tend towards unity. 


A short note of explanation is appropriate. The trans- 
formation from (4a) to (4b) is accomplished by a change 
in variable, £=wei"/*(2ka)"/®, After the transformation, 


it becomes evident that the halo term is a function of X 


only. The independence of the halo components dd ss; 
Him, on distance is obtained by approximating the F 


functions! by the integer 2. This is accurate if 


2poka/kd < 0.8 (6) 


_ which is justified for common values of kd. 


III. VALUES FOR THE KNIFE-EDGE 
AND HALO COMPONENTS 


The functions Ky. and Kim are plotted in Figs. 2-4, 
and H,, and Him are plotted in Figs. 5 and 6. Figs. 2 and 
3 give more accurate values for the knife-edge compo- 
nents for Z<0.7 and Fig. 4(a) and 4(b) for Z>0.7. Also, 


the halo functions are tabulated in Table I, and the 
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Fig. 5—Halo contribution to the diffraction field behind a mountain: 
real part (—H,e)\| for-horizontal and (H,.)1 for vertical polariza- 
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Fig. 6—Halo contribution to the diffraction field behind a mountain: 
imaginary part (Him)|| for horizontal and (—Him) 1. for vertical 
polarization, vs parameter X =wo(2ka) 8. 


knife-edge functions can be calculated by (3b) and (3c) 
with the help of Fresnel integral tables. None of these 
curves exhibit oscillations, most of them are monotonic. 
Although the individual components (1-S-C) and (S-C) 
oscillate, multiplication by the phase factor exp (jrZ?/2) 
removes the oscillations. 

The ratio (X/Z) occurs in (2). It depends only on 
(a/X) and (d/X), but not explicitly on Wo. Fig. 7 is a plot 
of constant values of (X/Z) ranging from 0.05 to 2 vs 
values of (a/X) and (d/d). 

The ranges of the parameters found in practice are 
2yo=0---0.2 radian, ba =0° *- 108 kd = 6X10? o>: 
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TABLE I 
REAL AND IMAGINARY PARTS OF HALO COMPONENTS FOR 
HorizonTat (||) AND VERTICAL (1) POLARIZATION 
x (Aye) || (Him) ll (Hye) eb (ign) a 
0 —0.554 0.148 0.534 —0.143 
0.02 —0.544 0.150 0.521 —0.146 
0.04 —0.539 0.153 0.512 —0.149 
0.08 —0.518 OFTS5 0.487 —0.152 
0.12 —0.505 0.158 0.468 —0.155 
0.20 —0.479 0.161 0.432 —0.157 
0.40 —0.427 0.159 0.370 —0.153 
0.80 —0.364 0.147 0.297 —0.135 
1.20 —0.328 ORASS 0.260 —0.119 
1.60 —0.306 0.128 0.237 —0.108 
2.00 —0.293 0.123 0.221 —0.098 
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Fig. 7—Values for the parameter (X/Z) = (2ka)1/3/(2kd /r)1/2 
vs (a/d) and (d/d). 


2X108, X=0---3,Z2=0---12and (X/Z)=0-.-- 
The graphs cover most of these ranges. 


De 


IV. EFFECT OF AN ISOLATED MOUNTAIN 
A. The Power of the Field 


Eqs. (2), (3b), (3c), and (4b) yield the real and imag- 
inary parts of the electric field. The intensity and phase 
of the field behind the mountain can be computed from 
these parts. When the mountain is the only obstacle 
between the stations, the effect of the ground between 
the stations and the mountain can be disregarded. This 
occurs when either the earth’s surface is too rough to 
support reflections or when the antennas’ directivities 
effectively discriminate against the earth-reflected path. 
The phase then is of no concern and an expression for 
the power is sufficient. From (2), the power diffracted 
by the mountain is° 


Po = [1/4(d: + d2)?]{ [Kre + (X/Z) Hye]? 
+ [Kim + (X/Z)H im|?} 
= [1/4(d1 + d2)?]{ | K |? + (X/2Z)2| A? 
+ 2(X/Z)|Kre Hre + KinHim}} (7) 


5 The basic transmission loss, L;, is related to Py by Ly =P o(d/4r)?. 
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where | 
| K|? = eae + ee | H|? = Aye oe HT gare 


B. Procedure for Calculation 


The procedure for the calculation of the power behind | 
a mountain given k, a, Di, Dz and Hi, Hz, when ground | 
reflections are unimportant, is as follows. | 
1) Calculate di, d:, and Yo. From geometry, these | 
parameters are related to the angles Wi and W2 sub-— 
tended by the mountain at the transmitter and receiver. | 


dy = dD, cos [(¥s —. V,)/2] 
+ (a — Hy) sin [(W2 — W1)/2] ~D, 
ds = Deo cos [(W. = W,)/2| 


— (a — H2) sin [(W2 — ¥,)/2] ~ De 
2vo = Wi + Wo > A/D, + Ae/ D2 
2) Calculate the mean distance d=2did2/(di+d:). 
3) Calculate three parameters in the following 
order.® 
X = po(2ka)?,  (X/Z) = (2ka)'/8/(2kd/m)+/?, 
Z= X/(X/Z) 
or use Fig. 7 for (X/Z). 
4) Calculate K,e and Kim for the Z value from (3b) 
and (3c), or read off their values from Figs. 2 and 3 or 4. 
5) Read off H,. and Him for the X value from Figs. 


5 and 6 or use Table I. 
6) Perform the operations in (7) to obtain the power. 


C. Numerical Examples 


Several examples are included in Table II. In all 
cases, the condition 2~ a/d <0.8 is satisfied. The exam- © 
ple in which a=10 miles, d=20 miles, 2¥)=2°, and — 


\=2 feet is calculated below in detail. 


It is found that ka=1.658X10°, kd=3.317X105, 
X = 1.208, (X/Z) =0.1505, and Z=8.03. All these quan- 
tities are calculated. (X/Z) may also be read off Fig. 7. 
From Fig. 4, rZK,.=1.005, 7ZK,, = —0.995. Hence . 
Kye= 0.0399, Kim= —0.0395, |K|?=0.003145. From | 
Figs. 5 and 6, H,.=(—0.328)), (0.259)1; Him = (0.1348)), 
(—0.1189), where | denotes horizontal and L denotes 
vertical polarization. 


Kye + (X/Z)Hre = (—0.00943)|, (0.0788) 1 
Kim + (X/Z)Him = (—0.01917)|, (—0.0574).. 


From (7), 4(di+d2)?Po = (0.000456)), (0.00950)... 

If the mountain were replaced by a knife edge, the 
power would be 4(d\+d,)?P,= | K| ?= (0.003145. Hence | 
the difference in db from knife edge introduced by the | 
crest of the mountain is (—8.38)), (4.80). db. Often this 
calculation can be avoided by the use of Fig. 8 as will 
be explained in a following section. 


ee eis Sequence avoids an indeterminancy in (X/Z) when 
X=Z=0 fora diffraction angle wo equal to zero. 
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TABLE II 


NUMERICAL EXAMPLES FOR THE EFFECT OF AN 
IsOLATED MOUNTAIN 


_ Example (a) (b) (c) (d) (e) (f) 

a@ miles 1 1 1 10 

d miles 50 20 50 2 |i0 | 3 

2Wo degrees 6 2 6 2 0.5 0.5 

d feet 2 2 20 2 D 1. 

ka/108 16.6 16.6 1.66 |166. 166. 1663. 

x 1.682 0.561 0.781 1.208) 0.302} 0.479 

x/Z 0.0442} 0.0698} 0.0651} 0.151) 0.213] 0.338 

Zz 38.0 Sam t)0 S-030 1.42 | 1.42 

ePo/Px)\ db |\—6.27 |=—5.10 |—6.83 |—8.38 |—3.22 |—5.19 
385 3.80 AS02241 03-00 


(Po/Px)1.db| 5.59 


0 =U NN. Go Ro Oo AN 


f 


2002 04 06 08 10 12 14 16 18 20 
—> X 


_ Fig. 8—Power ratio in db behind an isolated mountain compared toa 
| knife edge vs parameter X for various X/Z. The top set of curves 
are for vertical (1) polarization and the bottom set for horizontal 
(||) polarization. The asymptotic curves are for Z— % or X/Z—0. 


D. Limiting Cases 

he Eq. (7) has an asymptotic behavior for large Z, say 
Z>2.5. Large values of Z occur often as shown in the 
examples. In this limit, it can be seen from Fig. 4(b) 
‘that the knife-edge components approach 


(—Kim) — (Kre) > (1/nZ). (8) 
The expression for the power can be simplified to 

Py = {1/[2(di + de) wZ]?} [1 + 7X Hee)? 

| + (-1 + rXHim)?| 

> = 2{1/[2(d: + de)wZ]?} [1 + (w?X?/2) | a|? 
ete = He) |, ©) 
The dependence of the power on X and Z separate as 
factors. Denote the square bracket in (9) by 

‘ oS XH.) *b (1 + mCi ny. (10) 


The function f(X) gives the dependence on DG-n\emn 
‘tabulated in Table III for both polarizations. It is equal 


TABLE III 
POWER Function f(X) vs X (10) 


x [F(X] f(X)4 
0 1.000 1.000 
0.02 0.957 1.043 
0.04 0.916 1.085 
0.08 0.840 1.169 
Qr12 0.770 1.252 
0.20 0.649 1.412 
0.40 0.427 1.783 
0.80 0.202 2.423 
1.20 0.148 3.008 
1.60 0.210 3.587 
2.00 0.381 4.151 

TABLE IV 
Power FUNCTION g(X/Z) vs (X/Z) (12) 
X/Z [g(X/Z)]} [g(X/Z)] 1 

0 1.000 1.000 
OR 0.893 thea) 
0.2 0.792 1.226 
0.4 0.609 1.476 
0.6 0.454 1751 
0.8 0.324 2.050 
LD 0.221 Deots 
1.5 0.0778 3.289 
2.0 0.0992 4.362 


to the diffracted power relative to knife edge. Thus, for 
example (d) in Table II, the relative power can be 
found directly to be (—8.3) and (4.8). db for X =1.2. 

Another limiting case of (7) occurs for a zero diffrac- 
tion angle, that is for grazing angles. 


X/Z = (Qka)!8/(2kd/x)'?, X = Z= Kim = 0, 
Kre = 1 
Po = [1/4(d1 + do)2|[1 + As,u(X/Z) + Bi,W(X/Z)?] (11) 


where A1,|/=1.068, —1.108; B1i,!=0.305, 0.329. Eq. 
(11) is the same as (38) in the previous paper,’ and has 
been shown to agree well with experiment. 

Let 


g(X/Z) = 1+ As(X/Z) + B1,w(X/Z)?. (12) 
The function g(X/Z) is tabulated in Table IV. 


E. Graphs for the Power and Comparison with Expert- 
ment 


In general, the power in (7) multiplied by 4X (di +d)? 
is a function only of X and (X/Z), or X and Z. The 
ratio of the power relative to that obtained for a knife 
edge is plotted in Fig. 8 vs X for various (X/Z) values. 
The manner in which the power approaches the asymp- 
totes for small (X/Z) or large Z is evident. These curves 
in Fig. 8 can replace steps 4 to 6 in the procedure under 
Section IV-B if a high accuracy is not required. 

The total power relative to [1/4(d1+d2)?] is given in 
Fig. 9 vs X for various (X/Z). This graph is obtained 
by adding the contribution in db for the knife edge to 
the values in Fig. 8. These curves in Fig. 9 are effectively 
a plot of power vs scattering angle w)«X. When the 


346 


—> POWER x 4(4,+d,)° INdb 


12 14 


— xX 


as) 0.2 04 06.08 1.0 16 18 20 22 24 


Fig. 9—Power multiplied by 4(di+-d,)? in db behind an isolated 
mountain vs parameter X for various X/Z. 


radius of curvature or (X/Z) is large, as is usually the 
case in rolling country with relatively small differences 
in altitude, the theory can be used only for small dif- 
fraction angles, because of the approximation given by 
(6) in the analysis. This is not a serious limitation be- 
cause in rolling country large angles are rarely en- 
countered. Conversely, in an area with high mountains, 
larger scattering angles may occur. Usually at the same 
time the radius of curvature is smaller, and the theory 
is applicable to this combination of small radius and 
large angle of diffraction. 

The curves in Fig. 9 are redrawn on a larger scale in 
Fig. 10 and are compared with experimental measure- 
ments of the power vs scattering angles up to 19 degrees 
for both polarizations. The experimental apparatus has 
already been described.! The experimental range of ka 
is 0 to 250, of X is 0 to 0.8, and of (X/Z) is 0.1 to 0.51. 
In all cases, the initial slopes agree well with theory 
up to 5 degrees and in some cases 10 degrees depending 
on radius of curvature. 

For larger scattering angles the experimental curves 
continue as straight lines and the theoretical curves 
slope upward. 


V. EFFECT oF GROUND REFLECTIONS 
A. Amplitude and Phase of Field 


In cases when reflections occur from the ground be- 
tween stations and the mountain, more elaborate cal- 
culations are required corresponding to the well known 
four-ray theory.? Then Table I and Figs. 2 to 6 for the 
real and imaginary parts have to be used. The expres- 
sions for the power in Section IV are inadequate. 

The method for including reflections is as follows. The 
flat ground is eliminated by assuming that the reflected 
rays from the ground come from virtual images 71’ and 
T;', shown in Fig. 11, for the transmitter and receiver. 
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Fig. 10—(a) Comparison of theory and experiment for the power be- 
hind an isolated, perfectly reflecting cylindrical mountain up to 
diffraction angles of 19 degrees (vertical polarization). (The il- 
lumination factors for the transmitter-receiver antennas have 
been taken into account.) The arrows at the experimental curves 
indicate the value of X for which 2~0=5 degrees. (b) Similar to 
Fig. 10(a) (horizontal polarization). 


The four paths that the field is propagated, namely 
T1CT2, TyGiCT2, T1CG2T>, and 11GiCG2T>2 are replaced 
by four other paths, namely T,CT», iy CTs, 1 1Cl aan 
Ty'CT;'. The total field may then be calculated by com- 
bining the four components, taking into account their 


phases and amplitudes and the reflection properties of 
the ground. 


“a 
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Fig. 11—Geometry for a mountain with a reflecting ground. 


Let ai, a2, a3 and ay be the angles subtended by the 
-mountain with the horizon at 71, Ti’, T2, Ts’, respec- 
tively; Hi and H» be the heights below the top of the 
mountain for 71, T2; Di; and Dz be the distances from 
_ Ii and T>; to a vertical line through the peak of the 
- mountain; 6, and 6 be the heights above the ground for 
-T; and Ty» or below the ground. for Ty’ and 7%’; 2, 
22, 2¥3 and 2% be the scattering angles for the four 
_ paths mentioned above, respectively. 

Then the following relations are obtained: 


(WW = ar tas 22 = a2 + a3 
| hs =artaxg Ws = a2 + a4 = 22+ Hs — 1) (13a) 


a= sin— {a 
| [De +(e mye 


— sin7! ‘ Beh oe \ 
[Dit (@— HT" 


For ae, a3, and as, respectively, use H,+26;, Ae, 
- Hz+26, and Di, D2, D2 instead of H; and D, in (13b). 
_ For each of the four paths, there are different values for 
dh and ds. Since d; and d; occur in the phases, they must 
be calculated accurately and cannot be approximated 


(13b) 


f by D; and D2. 
= Let 
B:1 = (a3 — a1)/2 Bs = (as — a2) /2 

4 Bs — (a4 ae ay) /2 Ba 5 (a4 = a2) /2. (13c) 
_ Then the distances are given by 
dy; = D, cos 6B; +(e — Hi) sinB; fort = 1,3 

= D, cos; + (a — Hi — 261) sin B; fort = 2,4 

= D.cosB; — (a — H2) sinB; fora<2= 1,2 


dx; 
| D.cos 8; — (a — Hz — 282) sin 8; for i = 3,4 (13d) 


I 


_ The average distance can be calculated by 


d; = 2dyidei/ (dis + dei). (13e) 


For each path, values can be calculated from (13) for 
Wo, d, di, and dz, and then for X and Z. The total field is 


- the sum of these four rays. 
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Let R: and R, be the magnitudes of the reflection co- 
efficients of the grounds on the transmitter and receiver 
sides of the mountain. Let ¢1 and ¢» be the respective 
phase advances at reflection. Then the total field £ is 
given by 


E = BE, + RyEse-#) + RoEge-#? + RyRoEse-i +92), (14a) 


For a perfectly smooth reflecting ground, the equations 
reduce to 


R, — Rs ae ily. 
=0 for 1 polarization 
i = 2 


=f for | polarization, 
and 


B= by te Lig ieee te (14b) 


where the plus sign must be taken for 1 and the minus 
sign for 1 polarization. 

In (14), the expression for £;, (¢=1, 2, 3, 4,) is given 
by ; 


2(D: + Ds) E:~ exp [—jk(dii + dai) (1 + ¥2/2)] 
x { Kre(Zi) “hi GK im(Zs) 
+ (Xi/Z,)[Hre( Xs) + jHin(Xi)]}. (14c) 


The expression for the electric field already includes 
the amplitude and phase lag introduced by diffraction 
and the phase changes due to different path lengths d; 
and d». Since these expressions are combined in (14b), 
the phase factor k(di+d2)(1-+0?/2) must be included 
in the analysis in addition to the knife edge and halo 
terms. 


B. Procedure for Calculation 


The procedure for calculating the power behind a 
mountain, given k, a, D1, D2, Hi, He, 61, and 62 and in- 
cluding reflections from a perfectly reflecting ground is 
as follows: 


1) Calculate di:, dx, and pW; for the four rays from 
(13): 

2) Calculate the mean distances d; by (13e). 

3) Calculate the parameters Xi, (X/Z); and Z; for 
each ray. 

4) Read off (K;.); and (Kim); from Figs. 2 to 4. 

5) Read off (H,.); and (Hin): from Figs. 5 and 6, for 
each ray. 

6) Calculate for each ray the real and imaginary parts 
of the electric field given by the right hand side of 
(14c), including the extra phase introduced by the 
exponential factor. 

7) Sum the real parts and the imaginary parts of the 
four electric fields according to (14b). 

8) Square the sum of the real parts and the sum of 
thé imaginary parts, and add the squares. The 
result gives the value for 16(D:+D2)*Po. 
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VI. CoNcCLUSION 


The theory for diffraction by a mountain with a 
smooth cylindrical crest is shown to agree with experi- 
mental results up to larger angles of diffraction. The 
smaller the radius of curvature, the larger the angle. 
This situation usually is encountered in practical cases. 

When the stations are far enough away the integrals 
previously derived! depend only on one parameter. 
Hence the results of numerical calculations can be pre- 
sented by simple graphs and tables that simplify the 


Pattern of an Antenna on 
J. R. WAIT} anp 


Summary—Extensive numerical results are presented for the 
radiation fields of electric and magnetic type antennas mounted on 
smooth curved surfaces of finite conductivity. The model chosen is a 
circular cylinder whose surface impedance is specified. A residue 
series representation is employed for the portion of Space deep in 
the shadow while a geometrical-optical representation is used in the 
“fit” region. In the penumbra, the fields are expressed in terms of 
the ‘‘Fock functions.” The results are also applicable to other 
smoothly varying curved surfaces such as spheres, parabolic cylinders, 
and paraboloids. As an application, the E-plane patterns are com- 
puted for a small loop antenna on a spherical earth for both sea and 
land illustrating the so-called cut-back effect. 


INTRODUCTION 
yee the fact that the formal solution of the 


electromagnetic diffraction by a homogeneous 

sphere was obtained a half century ago, certain 
numerical aspects of the problem have yet to be com- 
pletely resolved. The difficulty usually originates in 
the extremely poor convergence of the rigorous har- 
monic series solution when the circumference of the 
sphere, ka, expressed in wavelengths, is large. Employ- 
ing function-theoretic means Watson! demonstrated 
in 1918, for a vertical electric dipole source, that this 
harmonic series could be converted to a more rapidly 
convergent series in certain regions of space. This repre- 
sentation has become known as the residue series as 
each term corresponds to a residue of a complex pole. 
This particular method has been exploited and refined 


* Manuscript received by the PGAP, April 27, 1958. 

rs rs Radio Propagation Eng. Diy., Natl. Bur. of Standards, Boulder, 
olo. 

*G. N. Watson, “The diffraction of radio waves by the earth,” 
Proc. Roy. Soc. London, vol. A95, pp. 83-99, October, 1918; pp. 546— 
563, July, 1919. 

2B. Vvedensky, “The diffractive propagation of radfo waves.” 
Tech. Phys., U.S.S.R., vol. 2, pp. 624-639, June, 1935: ; 


vol. 3, pp. 
915-925, November, 1936; vol. 4, pp. 579-591, August, 1937, 
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application of the theory to practical problems. 

For some typical cases, numerical results are given 
and it is shown that the diffraction field may be up to_ 
3 or 4 db greater than behind a knife-edge mountain of © 
equal height, when the field is vertically polarized, and _ 


up to 7 or 8 db smaller when it is horizontally polarized. 


This theory is one step nearer to explaining the ex- — 


perimental evidence found in field practice. Irregular 


shape of the mountain, vegetation or roughness of the © 


ground, conductivity, etc., are other determining fac- 
tors which have yet to be taken into account. 


a Curved Lossy Surface* 
A. M. CONDAf 


in the late 1930’s by Vvedensky? in Russia, Millington | 
in England, and Van der Pol and Bremmer‘ in Holland © 


The latter two authors in particular developed a highly 
accurate method to compute the values of the complex 
poles occurring in the residue series. Essentially, this is 
based on a representation for spherical wave functions 
by Hankel functions of order one third, and for this 
reason it has become known as the Hankel approxima- 
tion. It is more precise, although more involved, than 
the so-called tangent approximations of Vvedensky 
and Millington which were based on a second order 


representation of the spherical wave functions due to | 


Debye. Van der Pol and Bremmer have also demon- 
strated in an elegant fashion that a geometrical-optical 
representation for the fields can be derived as a special 
case which is valid when the source and observer are 
well within “line of sight.” In 1939 Gray® extended the 
Van der Pol-Bremmer treatment to a vertical magnetic 
dipole source. 

In 1941, Norton® presented a systematic method for 
calculating the ground wave field for both vertical and 
horizontal polarization. Using a number of graphs, the 


*G. Millington, “The diffraction of wireless waves round the 
earth,” (a summary of the diffraction analysis, with a comparison 
between the various methods), Phil. Mag., ser. 7, vol. 27, pp. 517- 
542; May, 1939. 

*B. Van der Pol and H. Bremmer, “The diffraction of electro- 
magnetic waves from an electrical point source round a finitely con- 
ducting sphere, with applications to radiotelegraphy and the theory 
of the rainbow; the propagation of radio waves over a finitely con- 
ducting spherical earth,” Phil. Mag., ser. 7, vol. 24, pp. 141-176, 
July, 1937; vol. 24, pp. 825-864, suppl., November, 1947; vol. 25, 
pp. 817-834, suppl., June, 1938; vol. 27, pp. 261-275, March, 1939, 

°M. C. Gray, “Diffraction and refraction of a horizontally 
polarized electromagnetic wave over a spherical earth,” Phil. Mag., 
ser, 7, vol. 27, pp. 421-436; April, 1939. 

°K. A. Norton, “The calculation of ground-wave field intensity 
over a finitely conducting spherical earth,” Proc. IRE, vol. 29, pp. 
623-639; December, 1941. 
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- method was in a form suitable for use by radio engineers. 


In 1949, Bremmer,’ in a monograph, presented a com- 
prehensive outline of the theory of ground wave propa- 


gation, listed a number of formulas for field computa- 
tions, and gave certain illustrative examples. 


For most applications of ground wave propagation, 


the above mentioned methods of Norton and Bremmer 


are satisfactory. The exception is when a line joining 
the receiving antenna and the transmitting antenna just 
grazes the earth. Norton has proposed that this case be 


‘treated by interpolating between geometrical optics 
and the field computed in the shadow from the first term 

of the residue series. Often this is quite adequate. In 
_ this instance, Bremmer recommends that the complete 


residue series representation should be employed. Un- 


fortunately, it becomes very poorly convergent near 
the boundary of light and shadow, particularly if the 
antenna heights are large. The convergence of the resi- 
_ due series is also poor at short distances at low and 
medium radio frequencies even when the transmitter 


and receiver are both on or near the ground. In the lat- 
ter case, an alternative series has been developed’:® 


_ which contains inverse powers of ka and the leading term 
is the flat earth formula of Sommerfeld. 


A more direct approach to the evaluation of the field 


near the light shadow boundary is to return to the orig- 
_ inal complex integral representation which is actually 
the first stage in the Watson transformation. This is the 
_ line of attack adopted by Fock and Leontovich®:!° in the 
_U.S.S.R. By making approximations to the spherical 
wave functions which were equivalent to the “Hankel 
approximation” of Van der Pol and Bremmer," the 


integral for the field was put in a form suitable for 
numerical integration. Fock!? gave numerical results 


only for perfect conductivity. Fock has shown that the 


- same integral arises in the diffraction of a plane wave 


by a paraboloid and in any generally convex surface if 


the radius of curvature is large and smoothly varying. 


It is the purpose of this paper to further exploit Fock’s 


' method to the computation of the field of a source on a 


aes 


curved surface whose conductivity is finite. Particular 
attention is paid to the case where the receiving antenna 


is at a large distance from the surface and the transmit- 


10, pp. 13-24; January, 1946. 
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Incident Plane Wave 


Fig. 1—Plane wave incident on a circular cylinder. 


ting antenna is on or near the surface. By utilizing rec- 
iprocity, the results will also apply when the location 
of the receiving and transmitting antennas are inter- 
changed. For purposes of presentation, a so-called cut- 
back factor is introduced which is essentially the E- 
plane radiation pattern of a small loop antenna located 
on the curved surface. The term “cut-back” has been 
used to describe the rather pronounced reduction of the 
field at low angles. 


. FORMULATION 


The model chosen as a basis for calculation is an in- 
finitely long circular cylinder of radius a as illustrated 
in Fig. 1 along with the appropriate cylindrical co- 
ordinate system. (The surface of the cylinder is defined 
by p=a and the gz axis is directed into the paper.) The 
incident plane wave which is incident normally on the 
cylinder has a z component of the magnetic field only. 
It is given by 


Hine — Hoe cos ¢ | (1) 
The boundary condition on the surface of the cylinder is 
Es = — ZHs)\,_0 


where Z is by definition the surface impedance. For a 
homogeneous cylinder of conductivity ¢, dielectric 
constant €, and permeability wu, Z can be well approxi- 


mated by 
1Ma) 
ZL= /— =H 
o + 1weE 


where 7 is the intrinsic impedance of the medium. The 
boundary condition in this form is usually ascribed to 
Leontovich!® who showed that for plane wave incidence 
it is valid if |n/no| 1 where no2£1207 is the intrinsic 
impedance of free space. More recently, Monteath,™ 
and Wait and Surtees!® have shown a more general re- 


(2) 


13 M. A. Leontovich, “Approximate boundary conditions for the 
electromagnetic field on the surface of a good conductor,” Bull. Acad. 
Sci., U .S.S.R., serie physique, vol. 9, p. 16; January, 1944. (In Rus- 

ian. ; 
a a G. D. Monteath, “Application of the compensation theorem to 
certain radiation and propagation problems,” Proc. IEE, vol. 98, 
t. 4, pp. 23-30; January, 1951. 
: & qk Wait i W. J. Surtees, “Impedance of a top-loaded an- 
tenna of arbitrary length over a circular grounded screen,” J. Appl. 
Phys., vol. 25, pp. 553-555; May, 1954. 
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striction in that the tangential fields should vary slowly 
in a distance equal to | y-| where 


y = [iww(o + iwe) }4/? 


is the propagation constant of the conductor. Actually, 
a somewhat improved formula for Z for ground wave 
propagation is given by 


Salina eh inh 


for vertical polarization.® 


FORMAL SOLUTION 


To obtain the field scattered from the cylinder, the 
incident field is first written in the form 


H,° = Hy >) ene™*!?In(kp) cos nd (3) 
n=0 

using a well-known addition theorem in Bessel func- 
tions. In the above €)=1, €,=2(n0) and J,(kp) is the 
Bessel function of the first kind of integral order . Since 
the secondary field H,** is to be a solution of the wave 
equation and is to give rise to outgoing waves at infinity 
it is written in the form 


H,8** = >> A,H,®(Rp) cos nd (4) 
n=0 
where H,,(kp) is the Hankel function of the second 
kind. The secondary field, H,°, is now chosen to satisfy 
the approximate boundary condition which can be re- 
written 


— zn| (5) 


noting that H,=H,'+H,. It readily follows that 


dis GJ; 
A, oom Hoeqe'"* t? t) Bi : «) (6) 
A," (x) + GH, © (x) 


where x=ka and G=—iZ/n. On the surface of the 
cylinder, the field can be written in the form 


Fh Eh at AF (x, ) (7) 
where 


€ne"™!2 cos nd 
TIX nad Hy (x) + GH, 2) (x) 


The voltage induced in a small loop placed at p=a is 
proportional to i eee: thus F(x, ¢) can be regarded as 
the radiation pattern of the loop in the principal or E 
plane. 


CoMPLEX INTEGRAL REPRESENTATION 


Unfortunately, (8) for the pattern function F(x, ) 
is very cumbersome for purposes of calculation unless 
x is reasonably small. For example, something of the 
order of 2x terms are required to secure 5 per, céntac- 
curacy. When x is large, it is desirable to represent 
F(x, @) as a complex integral which can be either evalu- 
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ated directly, or by deforming the cont6éur around the 
complex poles of the integrand. Such procedures have’ 
been used by Watson, Van der Pol and Bremmer, and 
Fock for treating the problem of ground wave propaga-_ 
tion from a vertical dipole on the surface of the earth. 
The application of these same techniques to cylinder | 
problems (which for large curvatures are essentially the i 
same as the sphere problems), has been carried out by _ 
Franz,'*—18 Sensiper,!® Baillin,?? and Wait.*»” In most of 
this work, with the exception of Fock’s, the complex | 
integral was evaluated only by a saddle point method | 
yielding useful expressions for the field in the “lit” re- | 
gion of space (2.e., | | <a/2). On the other hand, the | 
residue series representation is only useful deep in the 
shadow (z.e., |¢| >a/2). In 1946, Fock evaluated the 
complex integral for the perfectly conducting cylinder 
by numerical means and thus bridged the gap between 
the geometric-optical domain and the residue series 
solution deep in the shadow. The extension to imper- 
fectly conducting cylinders is carried out in what fol- 
lows. These results enable the pattern F(x, ¢) to be — 
evaluated for the whole domain of ¢. 
The complex integral representation of the function 
F for large cylinders can be written (see Appendix) as 


FS e-tetg(X) 9) 


where 


ext 


nes f di 
Se aol tg TE 


(10) 


where [., the integration contour, runs from © e-27/3 to 
0, then out along the real axis to ©. The other symbols 
are defined as follows: 


W1(t) is the Airy integral defined in the appendix, 


ka \ 18 ka\'8 T 
(OED 
2 2 3) 
ka\1/8 
plane (=) Z/n0. 


‘* W. Franz, “The green’s functions of cylinders and spheres,” 
Z. Naturforsch., vol. 9A, pp. 705-716; September, 1954. 

_ “W. Franz and R. Galle, “Semiasymptotic series for the diffrac- 
tion of a plane wave by a cylinder,” Z. Naturforsch., vol. 10A, pp. 
374-378; May, 1955. 

8 W. Franz and P. Beckman, “Creeping waves for objects of 
finite conductivity,” IRE Trans. oN ANTENNAS AND PROPAGATION, 
vol. AP-4, pp. 203-208; July, 1956. 

*”S. Sensiper, “Cylindrical radio waves,” IRE TRANS. on AN- 
TENNAS AND PRopaGArTIon, vol. AP-5, pp. 56-70; January, 1957, 

UA AL, Bailin and R. J. Spellmire, “Convergent representations 
for the radiation fields from slots in large circular cylinders,” IRE 

TRANS. ON ANTENNAS AND PropaGatTion, vol. AP-5, pp. 374-382: 
October, 1957. 

a hl Wait, “Radiation characteristics of axial slots on a con- 
ducting cylinder,” Wireless Engr., vol. 32, pp. 316-323; December, 

1955, and “Pattern of a flush mounted microwave antenna,” NBS 
J. Res., vol. 59, pp. 255-259; October, 1957. 

_ 2 J. R. Wait and J. Kates, “Radiation Patterns of Circumferen- 
tial Slots on Moderately Large Conducting Cylinders,” Monograph 
No. 167 R, IEE, London; February, 1956, republished in Proc. IEE, 
pt. C, pp. 289-296; February, 1956. 
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The above formula for F is only valid if (ka/2)!/*1>>1 
and | ¢| <1. This means that attention, at least for the 
moment, is restricted to large cylinders and to the region 
near the shadow boundary (i.e, the penumbra). In 
terms of the electrical constant € and o the complex 
factor g is conveniently written 


1/2 
q = ADMs_-ir/4 E re i=] 


o 


A = (ka)¥3 (=)". 
\ 20 


iz When the conductivity is sufficiently large (1.e., o>ew) 


q S& ADNS earls, 


where 


(11) 


_ The geometrical-optics approximation for the func- 
tion F can be written 


= [1 a Rle-*« sin @ (12) 


_ where R(¢) is taken equal to the Fresnel reflection co- 
efficient for a plane wave incident on a plane boundary 
_ whose surface impedance is Z. Therefore,” 


_ sind + Z/n0 8 + Z/no (13) 
sin 8 — Z/n. 6 — Z/n0 
and consequently, 
ae 20 
F & e-ikabgikao?/3y —__ (14) 
i= Z/n0 


for small values of 6. This equation can also be found 
- from a saddle point evaluation of the integral g(X) 
~ when X is negative and not near zero, 

| In terms of X and the electrical constants o and ¢, 
_ the preceding equation becomes 


F = e~ikabeix*/3 Tea (15) 


—1/2 
(1 — ax — 22/84 (1 + i=) 
o 


~ remembering that X = (ka/2)'6. 

Up to this point the discussion has referred specifi- 
cally to a circular cylinder of infinite length. On intuitive 
grounds one would expect the function F to characterize 
the pattern in the penumbral region of an antenna on 
any smooth surface whose radius of curvature is large. 
The quantity a is then regarded as the principal radius 
of curvature. It appears to be difficult to justify the pro- 
ceeding statements on a rigorous basis, for an arbitrary 
smooth surface. As shown by Fock, however, the con- 
cept is valid when applied to a perfectly conducting 
paraboloid and a similar demonstration was made by 
Rice" for the parabolic cylinder. 


23 [t should be noted that if Z is replaced by (?— 102 cos? 0)*?, R 
becomes the Fresnel reflection coefficient for a homogeneous half 


%S, O, Rice, “Diffraction of plane radio waves by a parabolic 
cylinder,” Bell Sys. Tech. J., vol. 33,-pp. 417-502; March, 1954. 
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For applications to antenna radiation over a spherical 
earth, it would be more logical to employ a model of a 
vertical electric dipole located on a homogeneous sphere 
with a specified surface impedance Z. This was the ap- 
proach used essentially in the work of Van der Pol and 
Bremmer, Norton, and Fock. In the present notation, 
the Hertz function (which has only a radial component) 
II of the vertical dipole of unit strength located on the 
earth’s surface is given by 


(16) 


where @ is the angle at the center of the earth subtended 
by the source dipole and the receiver at height /. In the 


above, 
106 Wat — 
ee (/—f griet SAAS Si) (17) 
a JT, Wi'() — qwilt) 
is an attenuation function and 
ka 1/3 _ P 
z=|{ — 7) 18 
ee 08) 
and 
) 1/3 
y= (=) Rh. (19) 
ka 


Now when y is large, the Airy integral function Wi(t—y) 
can be replaced by the first term of its asymptotic ex- 
pansion for the important range of ¢ in the integration. 
Furthermore, if x is large in such a way that x—~/y is 
finite then 


ull? 


pw, p28) 98/2 
eS eek v3/29(X) 


(20) 
when X=x—v+/y is a finite parameter. It is thus shown 
that II expressed as a ratio to the free space Hertz func- 
tion IIp can be written 


He ~F n SOR2T) 
IIo 
where F is identical to the cut-back factor for the cir- 
cular cylinder and X has the same geometrical meaning. 
If the source dipole were at finite height (say fo) then 
the function in (17) should be replaced by*”° 


Lees . 
V =—A/— if em Wit — mo] Wee sy) 
2 M1) Ts 
W(t) — qWolt 
AO TS wie — 9) |e 02 
Wi'(t)—qWMi(t) 
where 


2 1/3 
yo = (—) kho (23) 
ka 
and W,(é) is an Airy integral function. It was shown in 
an earlier communication,® however, that quite gener- 
ally 
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FACTOR 
=| E=field at P'for source at Q 
on spherical surface 
: E-field at Pfor source at Q 
in free space ; 
== art] — 
= 
= a 
- Rateenics 
5° 
c2) 
eee 
x 
: = 
2 —— 
s 
oO 


(24) 


subject to yo being small compared with unity and hy 
small compared with a. 

Summarizing, it can be said that the cut-back factor 
is applicable to a transmitting vertical dipole at height 
hy and a receiving vertical dipole at height h subject to 


kh K a 


and 
kh > (ka), 


When kh becomes comparable to |/Z| it is desirable 
to define a modified cut-back factor F,, as follows 


Z 
Fin = E ale tere bho Ff (25) 


No 


which is valid if kayv<(ka)*/8 and hy<ad. 


DISCUSSION OF THE CuT-BACK FACTOR 


Using the preceding formulas for ew/s=0 or a=0, 
the function | F| described as the cut-back factor, is 
plotted in Fig. 2 as a function of the angular parameter 
X for A ranging from 0 to 10. The geometrical optics 
representation is used for large negative X and the 
residue series representation (see Appendix) for large 
positive X while in the intermediate range, the integral 
g(X) was evaluated numerically. The corresponding 
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phase factor A is plotted in Fig. 3(a) and 3(b). It is de- 
fined as follows 


F = | PF le-tareix*/gia for 9 < 0 (26) 


and 


F = | F| e~*ei4 for 9 > 0, (27) 


This method of normalizing the phase enables one to 
regard A as a phase difference between that of the inci- 
dent field and that of the resultant field in the lit re- 
gion (i.e., <0). On the other hand, A is the phase dif- 
ference between the electrical arc length ka@ and the re- 


y 
ey, 
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Ba aes t 
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THE LOSS FACTOR 
“A= 8(ka)!”? nie as 
| where p= (S92) | 
4 . | & 
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sultant field in the shadow (i.e., 9>0). To facilitate the 


application of these curves to radiation pattern compu- 


tation for a source on a spherical earth, the loss factor 


\" 


A is shown plotted in Fig. 4(a) as a function of the re- 


- ciprocal of the ground conductivity, 1/o, for frequencies 


y 


from 20 to 5000 kc. To allow for normal atmospheric 


- refraction the radius a is replaced by 4/3 times the ac- 
_ tual earth radius. 


At the higher frequencies and when the earth is 


_ poorly conducting, the displacement currents in the 


ground can be appreciable. Following Belkina” the 


~ complex factor g is rewritten in the form 


— ip®!6 
TS ep — 93" we 
4 which is equivalent to (11) if 
| p= [ass 


2% M. G. Belkina,. “Tables to calculate the electromagnetic field 
in the shadow region for various soils,” Soviet Radio Press, Moscow, 
U.S.S.R.; 1949. 
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and 


e/€0 
(6070a)?/® 


€ DNAs 
ee =( ) 
€0 aano 


This seemingly complicated representation for g has the 
advantage that the displacement currents are described 
by the parameter a which is independent of frequency. 

To facilitate the conversion from the electrical con- 
stants o and e, the factor a@ is plotted as a function of 
1/o for various values of the relative dielectric constant 
e/€) again assuming an effective earth radius factor of 
4/3. Numerical values of | F | and A are listed in Tables 
I and II respectively for A ranging from 0 to 10 and 
a from 0 to 0.03. For the smaller A values, the influence 
of a is negligible and, consequently, only the results for 
a=0 need be shown. 

The cut-back factor | FI is shown plotted in Fig. 
5(a) and 5(b) for sea water and average land respec- 


tively. The frequency varies from 20 to 5000 kc. The elec- 


trical constants for sea water are taken as ¢=5 and 
e/€9=80 and for land they are taken as o=0.005 and 
e/é9=15. As an interesting illustration | F| is shown 
plotted in Fig. 5(c) for both land and sea at 100 kc. The 
corresponding geometrical-optical approximation for 
| F | is also shown. At the larger negative values of 0, 
the two sets of curves merge together as they should. 
As @ approaches zero, however, they diverge signifi- 
cantly. It is noted that geometrical-optics would pre- 
dict that | F| tends to zero at the shadow boundary 
(z.e., as 0-0). 


EXTENSION TO H-PARALLEL POLARIZATION 


Up to this point, the incident wave is polarized such 
that the magnetic field is parallel to the axis of the 
cylinder. The other case of interest is when the electric 
field is parallel to the axis of the cylinder. The conver- 
sion of the former to the latter is simply obtained by re- 
placing H, by E, and Z/no by m0/Z. Therefore, if the 
incident field is given by 


Eine = Ey etke cos ¢ — Eve sin 6 (29) 
the resultant tangential field on the cylinder is 
HO ne a OY dg (30) 


where 


ext 


1 
x = dt | e~ika8 31 
Ate feat 1, Wi) — Wi) ile Oy 


ka\33 mo 
r-(8)"S 
2 vhs 


Since in most applications | Z/no|<«1 and ka>>1, the 
quantity | a*| is a very large number and F* can be 
well approximated by 


with 
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TABLE I 
NUMERICAL VALUES OF | F| 
= =0. =0. =03 | A=05 | A=0.7 =1.0 A=1.2 A=1.5 A=2.0 
x are es = : os ab 3 es a a= a= a=0 0.03 a= a=,03 a=0 a=0.03 a=0 a=0.03 
—3.0 | 1.998 | 1.95 1.90 1.85 i715 1.67 1.55 1.54 1.48 1.47 hea eeie seiey IE 
—2.5 | 1.994 | 1.94 1.88 1.82 1.71 1.61 1.47 1.47 1.39 1.38 1.28 Alege Zhe et eTe 1.13 
22950 M1 (0820) A192 1.84 1.78 1.65 1.53 1.38 NeW 1.29 1.28 A toe et 16 eel eae OF 1.01 
Sea er Ilr 1.79 1.71 1.56 1.44 1.25 1.24 1.14 1.13 1.02 | 1.01 | 0.869 | 0.868 
=1.0 | 1.861 | 1.82 1.70 1.61 1.44 1.26 1.13 leo 0.968 0.962 0.850 | 0.840 | 0.707 | 0.706 
—0.5 | 1.682 | 1.59 1.54 1.43 1.27 1.09 0.918 0.917 0.771 0.765 0.669 | 0.664 | 0.539 | 0.537 
0 1.400 | 1.34 7 1.15 1.01 0.901 0.711 0.710 0.579 0.570 0.490 | 0.485 | 0.386 | 0.384 
0.5 | 1.059 | 1.00 0.952 0.904 0.801 0.659 0.500 0.499 0.401 0.339 0.317 | 0.316 | 0.301 | 0.300 
1.0 | 0.738 | 0.722 | 0.692 0.661 0.563 0.446 0.316 0.315 0.251 0.248 0.181 | 0.180 |0:119 | 0.117 
1.5 | 0.488 | 0.477 | 0.466 0.444 0.384 0.290 0.196 0.195 0.152 0.149 0.102 | 0.0997 | 0.0608 | 0.0608 
2.0 | 0.315 | 0.312 | 0.308 0.296 0.249 0.184 0.117 0.116 0.0838 | 0.0829 | 0.0528 | 0.0523 | 0.0282 | 0.0288 
2.5 | 0.202 | 0.202 | 0.202 0.194 0.162 0.114 0.0668 | 0.0664 | 0.0453 | 0.0449 | 0.0266 | 0.0261 | 0.0122 | 0.0127 
3.0 | 0.130 | 0.130 | 0.130 0.127 0.104 0.0706 | 0.0384 | 0.0376 | 0.0246 | 0.0243 | 0.0129 | 0.0127 | 0.00489 | 0.00522 
x A=3 A=5 A=7 =10 
a=10 a=0.03 a= a =0.03 a= a =0.03 a= a=0.01 a=0.02 a =0.03 
Se eOmlLAeO2 1.01 0.768 0.743 0.588 0.586 0.444 0.442 0.450 0.455 The third significant figure is doubtful 
—2'5 | 0.919 | 0.916 | 0.651 0.647 0.512 0.510 0.379 0.382 0.388 0.394 in some cases. 
—2.0 | 0.804 | 0.797 | 0.558 0.555 0.432 0.431 0.315 0.317 0.322 0.326 
—1.5 | 0.664 | 0.662 | 0.459 0.457 0.344 0.344 8.246 0.248 0.253 0.256 
—1.0 | 0.520 | 0.518 | 0.355 0.354 0.255 0.257 0.181 0.183 0.185 0.187 
—0.5 | 0.685 | 0.684 | 0.249 0.249 0.175 0.182 0.122 0.124 0.126 0.127 
0 0.267 | 0.266 | 0.158 0.159 0.110 0.115 0.0742 | 0.0760 | 0.0778 | 0.0791 
0.5 | 0.139 | 0.142 | 0.0779 | 0.0788 | 0.0535 | 0.0560 | 0.0368 | 0.0375 | 0.0386 | 0.0399 
1.0 | 0.0646 | 0.0661 | 0.0339 | 0.0355 | 0.0229 | 0.0246 | 0.0158 | 0.0162 | 0.0166 | 0.0176 
1.5 | 0.0306 | 0.0311 | 0.0149 | 0.0156 | 0.0103 | 0.0111 | 0.00639 | 0.00668 | 0.00695 | 0.00714 
2.0 | 0.0129 | 0.0132 | 0.00603 | 0.00617 | 0.00389 | 0.00426 | 0.00257 | 0.00263 | 0.00275 | 0.00288 
2.5 | 0.00528 | 0.00547 | 0.00247 | 0.00261 | 0.00204 | 0.00216 | 0.00105 | 0.00108 | 0.00112 | 0.00115 
3.0 | 0.00219 | 0.00225 | 0.000898 | 0.000949 | 0.000689 | 0.000693 | 0.000361 | 0.000369 | 0.000385 | 0.000407 
TABLE II 
NUMERICAL VALUES OF A (IN DEGREES) 
(VALUES ARE NEGATIVE ExcEPT WHERE INDICATED BY +) 
ce eR ee ne AT 
x | 4=0 |4=0.1| A =0.2| A =0.3| A =0.5| A =0.7 A=1.0 A=1.2 A=1.5 A=2.0 
a=0 | a=0 | e=0 | a=0 | a=0 | a=0 | e=0 0.03 | a=0 0.03 2=0 0.01 0.02 0.03 a=0 0.01 0.02 0.03 
3-0 Oh 0.52 0.941 252} > 3.619 6.14] 8:4 | 10.9] 90.7 42:43,1 | 13.0 |) 15.6) 49.3.) 45.4) 14.8) dace. 48.0) oars 17.6 
Sos ah 0.85) (0.91) 2-4 329°] 7.0.) 9.4 |) 12.6) 12.4 | 14:4 |. 14.9] “7.244 g6'01 4674 46.4 p94) aoes 19.4 19.3 
=2.0 | 41.5] 0.85) 2.7] 4.4] 8.0] 10.3| 14:4] 14:1] 15:2 | 14:9| 1910] 18:8] 18:6] 18:3] 21:9] 21:7] 2175 21.2 
Siete 2.2) 0,74). 3.01" 8.3 | 9:4.) (42°95) 17.0.5 46.6] 19,2| 18.8) 21.7 | 34.8 | oiean S101 gee oat 33.5 23:0 
Bon oesa7 |) 0063) 4.0 BA > 4252). 160001! 2063.) 19.851) 23.2). 22.8 | 95.8 | 25.5\ 1995-4 los ols 90g. ge eee 26.9 
Seo ers9\) fF a 6.5) 10,3.) 18.1 23.20) 382-1 275 | 33.01 3256 | 38.7 40 a8 ale Sake Seiad > ee wae eens ce 40.0 
0 0 O30 H) S13. 119.8 3125 | 44.81 50.9 11.4917 -.| 9'60.8-| 50.3 | 68.3 | 47 l6-_66.9 ) 66.0% eo eats 67.2 66.1 
0.5 | 11.4] 20.7] 32.0| 38.8] 53.5] 67.3] 83.0 | 79.0] 92.3 | 90:8 | 100.0 | 101.1 | 102.2 | 103.4 | 106.2 | 104.1 | 102.7 101.3 
1.0 | 26.6] 40.3] 53.9 | 61.3] 81.2} 98.4 | 119.0 | 116.0 | 128.4 | 124.4 | 137.6 | 136.0 | 134.8 | 133.1 | 144.1 | 141.8 | 14010 | 1380 
1.5 | 42.7 | 59.5] 76.9| 86.3 | 112.4 | 132.8 | 158.8 | 155.5 | 170.0 | 166.8 | 178.2 | 178.2 | 177.0 | 175.5 | 186.5 | 184.6 | 182.7 180.9 
2.0 | 58.0 | 77.9 | 99.0 | 110.7 | 141.1 | 167.4 | 198.6 | 195.0 | 212.0 | 210.0 | 221.5 | 220.3 | 219.1 | 217.9 | 227 : : 
-1| 225.5 | 223.9 | 299.4 
2.5 | 72.9 | 96.5 | 116.0 | 134.8 | 170.5 | 202.0 | 236.8 | 234.8 | 251.0 | 249.2 | 262.1 | 263.9 | 262/5 | 261.1 | 269.8 | 268.2 | 260.6 | 265.1 
3-0 | 87.6 | 114.9 | 142.7 | 158.7 | 199.9 | 236.3 | 276.9 | 274.5 | 293.2 | 290.5 | 309.2 | 307.6 | 306.0 | 304.4 | 313.5 | 311.6 | 309.4 | 307.5 
A =3.0 A=5.0 A=7.0 = 
x a=0 0.01 0.02 0.03 a2=0 0.01 0.02 0.03 a=0 0.01 0.02 0.03 2=0 ob O02 0.03 
pe On 28.41 22.8 | 22.2) (20.4) 20.4 | 27:61 "2611! 24,51 33.0 | 30.9] 98.8 | 96,9 | 35.4 9 41.9) ee peas 5 
Soh 2552) 24,4 | 23.6), 22.9) > 30.7 | :2914,| 27.5 | 25.8 1 33.6 | 34.5 | ocd}. oF 3 : ; ; ; ea eae 
Bee eri) 26.3) 25.51 24-7 |" 32.21 30.6 | 2920 |. 27.3 | 34.7 |) 32.6 | sonssec an a | - at al aecaae ane $4 € \ oct Gee ee 
T1-5 | 29.2 | 28.6] 28.0] 26.3 | 33.2] 31.5 | 29.9] 28.2] 36:1 | 34:0] 31.9| 29:8] 3810| 33.2| 20:0] 2416| doubthtin 
1.0 | 37.0) 36.1] 35.2 | 34.2) 37.6 | 36.1] 34.6 | 33:2] 42:2] 40:1| 38:0| 35.9| 43.2| 39:1| 33:8] 28:9 arth ati gat > 
—0.5 | 54.2) 52.5) 50.8) 49.1 | 50.8 | 48.9] 47.0] 45.0] 55:3: 53:1 | 51.0] 48.9| 56.0] 51.2| 45:7| 42.6| “> 
0 Bi ez | 79.8 17.8 | 76.0) 75.7) 73.1 | 70.5 | 68.0 | 77.2 1.95.1.) 72.9 | TOR | 77.2) 921 67 | es.e 
0.5 | 118.0 | 116.0 | 114.0 | 113.7 | 107.9 | 105.2 | 102.6 | 99.9 | 106.8 | 104:7 | 102/5 | 100.4 | 104.9 | 99.9 | 04°75 | go's 
1.0 | 146.0 | 143.9 | 142.0 | 140.2 | 143.0 | 140.4 | 137.8 | 135.1 | 140.2 | 136.8 | 133.4 | 130.0 | 138.2 | 134.1 | 1280 | 12270 
1.5 | 187.0 | 184.0 | 181.4 | 180.3 | 181.1 | 178.5 | 175.9 | 173.3 | 178.2 | 174.5 | 170.8 | 167.1 | 176.8 | 171.2 | 165.3 | 150°0 
2.0 | 226.0 | 224.0 | 221.8 |,219.8 | 219.7 | 216.8 | 213.7 | 210.9 | 215.8 | 212.0 | 20811 | 204.4 | 213.6 | 206.5 | 20271 | toy > 
2.5 | 268.1 | 265.5 | 263.0 | 260.0 | 257.6 | 254.3 | 251.0 | 247.8 | 253:9 | 249.8 | 245.7 | 241.6 | 250.8 | 243.9 | 23875 | 239 °5 
_ 3:0 | 306.2 | 304.1 | 302.0 | 300.1 | 296.1 | 292-9 | 289.7 | 286.4 | 290.7 | 296.4 | 29211 | 277.9 | 286.4 | 280.1 | 274.3 | 269.1 
ee 1 where hy is the n lized ic fi 
FY — — U(X) et (32) $ ormalized magnetic field. For purposes 
Vr of presentation, it is useful to introduce a phase factor 
where A;, in analogy to A, as follows 
en iXt ; 
er snl eae , 
u(x) =f —— an (33) he = | U(X) | e-MererTsei for @< 0 (34) 
Vr Jr, W(t) 
and 
The integral U(X) has been evaluated numerically by pes | U(X) | ited gid, ¢ 
Rice.*4 Noting that sy CN eS OTe: (35) 
1 According to geometrical-optics 
Hy, = 7 Dp at p= 2 
Hy = — (sin 0) (Eo/no) eta sin ot p=a (36a) 
. . laws 4 Ay} . 3 ! 
it is seen that = — O( Eo/no)e—#a%eika0?/3! (36b) 


no(ka/2)'!8H 
eke >— ES aT BK OG) = he 
0 


which is equivalent to 


hy = — 2X, (37) 
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Fig. 5—The cut-back factor as a function of angle above or below the horizon. (For great heights 
this angle can be determined adequately by ray tracing methods.) 
_ The function | U(X)| which is proportional to the am- os) ka \'/8 
; : ; & enika = 
plitude of the tangential magnetic field, | HZ,| is shown BS e's) 2 
plotted in Fig. 6(a) as a function of X on a semi-log by 
scale. The corresponding geometrical-optics approxima- 4. gn ika(n—8) g | — 6) (=) | (38) 
2 


tion is also shown. The function | U(X)| along with the 
phase factor A, is plotted in Fig. 6(b) using linear scales. 
In this case, the geometric-optics result is simply a 


straight line and the phase is zero. 


APPLICATION TO MODERATELY LARGE CYLINDERS 


In all the preceding discussion the domain of @ and 
the radius of curvature is large enough that the tan- 
gential magnetic field is described adequately by (9). 
When ka is not excessively large, the formula for F 
should be replaced by (53) in the Appendix. The leading 


term of this expression which is usually adequate, reads 


and is valid for the domain 
r/2=|6|>—A6 where (A6)*<1. 


The second factor in the above equation can be inter- 
preted as a wave which has “crept” around from the 
rear of the cylinder. To illustrate the behavior for 
cylinders of moderate size | F| and A based on (38) are 
plotted in Fig. 7(a) and 7(b) respectively for ka =8 us- 
ing the normalization for phase defined by (16) and 
(17). Calculated points from the rigorous harmonic 
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series solution [#.e., (8)] are shown in Fig. 7(a) and 
7(b). In view of the fact that ka is not large, the agree- 
ment between the two sets of calculations is very reason- 
able. The corresponding curves for ka=21 are shown in 
Fig. 8(a) and 8(b) where the agreement between, the ap- 
proximate and the rigorous theory is much improved. 


CONCLUDING REMARKS 


The numerical results presented herein can be used 
to predict the radiation patterns of electric and mag- 
netic type antennas located on smooth curved surfaces 
of finite conductivity. For example, in certain investi- 
gations of ionospheric wave propagation, it is of im- 
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portance to know the amplitude of the incident wave at 
the lower edge of the ionosphere. The curves of the cut- 
back factor can be used to estimate this quantity. 
Furthermore, the field measured at the receiving an- 
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tenna on the ground can also be estimated from the cut- 
back factor assuming of course, that the amplitude of 
the downcoming wave is known. It is particularly im- 
portant to use these antenna cut-back factors in the in- 
terpretation of long distance ionospheric sky wave prop- 
agation at iow radio frequencies. Considerable error may 
result if the radiation patterns are simply computed on 
the basis of a flat earth. 
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An excellent discussion of the application of cut-back 
factors to the calculation of sky wave field strength has 
been presented recently by Norton.” 


MATHEMATICAL APPENDIX?’ 


The far field pattern in the equatorial plane or E 
plane located on the surface of circular conducting cyl- 
inders was shown to be described by the function 
F(x, 6) where x=ka and ¢ is the angular or azimuthal 
coordinate. This function is given by 


pO £5 
F(a, ¢) = — 
ON pe oe. HME GHae 


eint!2 ging 


(39) 


where G= —1i(Z/no), Z is the surface impedance and 
noaz120r. 

The summation is over both positive and negative 
integers. H,,© (x) is a Hankel function of the second kind 
and the prime indicates a derivative with respect to x. 
The first step is to express F(x, @) as a contour integral 
in the complex v plane. The contours chosen are the lines 
Cz and C; which are located just above and just below 
the real axis of vy as shown in Fig. 9(a). Some considera- 
tion shows?’ that 


1 2 
F(x, o) = (-)= 


e vt /2eve 
; f ; dy. (40) 
cc, [Hy ®! (x) + GH, (x)] sin vr 
This integral has poles at yv=0, +1, +2, 43,--- as 


indicated in Fig. 9(a) as equispaced points on the real 
axis of v. It can readily be verified that the sum of resi- 
dues of the integrand at these poles leads back to the 
series expression. 

Replacing vy by —y in the integral over Cy, it is seen 
that 


dv ei (o—4 12) 
‘p> sin yr [H,®! (x) + GH, (x)] 
d en (o—-4/2) 
3 f : (41) 
c, sin vm [H_,®’ (x) + GH, (a)] 
and since 
H_, (x) = eH, (x) 
the integral over C2 is equal to 
iv (84 /2—6) 
: (42) 


dy 
fe sin yr H,’(x) + GH, (x) 


%K, A. Norton, “Low and medium frequency ionospheric 
propagation,” sec. 4.2 in “Transmission Loss in Radiowave Propaga- 
tion,” to be published. See also A. D. Watt et al., “Observations on 
some LF propagation paths in Arctic areas," to be published. 

27 MKS units and a time factor exp (zwt) are employed. 

28 J. R. Wait, “Blectromagnetic Radiation from Cylindrical Struc- 


tures,” to be published as a monograph. 
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Fig. 9—Contours in complex plane. 
Consequently, 


it eur e—¥ (¢—1/2) + e7 Br/2—9) 
P(asi6) = ap | ow (43) 
mé&Jo,sin vel H,'(x) + GH,® (x) 


For this contour, Im.v <0 and therefore, it is permissible 
to write 


eur 


—— = 2j De e- 2am (44) 
sin vr Sar 
which leads to 
Vip 1 3a 
TX m=0 D; 2 
where 
e~ iv (B+2rm) 
rao) = f 
4 auths!)(a) ORO eee 
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with f being either 6 —7/2 or 30/2 —o. 
The so-called third order approximation for the 
Hankel functions is now introduced. This enables one 


to write?’ 


1 a acd _ 
H,@(«) = + (=) W(t) . (47) 
T x 
and . 
i fit ee 
H,@" (x2) = — (=) Wi (t) (48) 
ra x 


where Wi(t) and W,’'(t) are Airy integrals defined by”® 


1 - 
W(t) = — fh ET Ri (49) 
eile IT, 


and 


1 3 
W(t) = —{ sdset— ti? (50) 
T r, 


in the notation of Fock where the quantity ¢ is defined | 


by 
Mba ices 


In the foregoing text a function W2(t) was also em- 
ployed. It has the same form as Wi(#) but the contour 
I, is replaced by I’; shown in Fig. 9(b) and 9(c), respec- 
tively. The function T,,(8) can now be written 


ar ae 

2 ee ee (51) 
where 

Core yp oes (52) 

§ a Wty = STN : 

Vrdr, W(t) — qW (0) 
with 
gq = — i(x/2)!8Z/n9 
and 
x \W3 
NER raed (8 = 2mm)( =) "4 
Therefore, 
00 a x 1/8 
F(x, ¢) ~ > e— tka ((p—m /2)+20m) g E — — + 2am (=) | 
m=0 2 - 
Pa NLS 
ae Se ee de mL | E — d+ 21m (=) | (53) 
m=0 2 : 

; ie The Airy integral a related to the Hankel function of order 
3 DY 


Wilt) = caning (—1)¥2 Ay), [ (5) (9s | : 
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which is valid for =| | =x/2—Ad where (Ad)?<1. 
It is also readily shown that if m(x/2)!/*>>1, only the 
term for m=0 need be retained. 

For the shadow region (X >0) it is desirable to cal- 
culate the integral for g(X) by the evaluation of the resi- 
dues at the poles of the integrand. The poles #, of the 
integrand are determined from 


Wy' (ts) = qgWil(ts) = 0. (54) 


Expanding the denominator in a series around the pole 
t=t,, leads to 


Wii) — qWi) & [Wi''(ts) — gW'(te) 1 — te) 


when terms containing (t—#,)?, (t—é,)? have been 


omitted. But since 

W1'(ts) — isWilts) = 0 
and 

Wi'(ts) — qWilts) = 0 
_ one may write 


(55) 


Wy (te) — qWi'(ts) = (ts — gq’) Wis) 
and consequently 
Wi'(t) Ry qW x(t) = (ts a gq?) Wi(te) (t = i) ae (56) 


higher order terms in (t—t,). 

For X>0, C can be deformed to a path enclosing the 
poles which are contained (as will be shown) in the 
fourth quadrant of the t plane. g(X) is thus equal to the 
sum of the residues of the poles at ¢=#,. Thus 

en iXte 


rey (57) 
: T(t — @)Walh) 

The numerical evaluation of the poles ¢,, being a solu- 
tion of 


Wi'(t) — qWilt) = 0 


- has been discussed in detail by Fock. Actually, how- 
ever, the task had been completed by Van der Pol and 
Bremmer a decade earlier. (The equation for the roots 
i, arises also in the problem of a dipole located on the 
surface of an imperfectly conducting sphere.) Van der 
Pol and Bremmer’s equation is 


Heys ([)(=2r.)*2] 


6 in/3 __ ee (—2r,)71/? (58) 
Sys ((4)(— 20)" 
where 
- (no/Z) 
6=]-121 ; 
(ka)*!* 


Ho and Hin™ are Hankel functions of the second 
kind of order 2/3 and 1/3. But, since 
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Wil) = mi y/’ z (=a (=) (—o| (59) 


it readily follows that 


eee EME O. 


(60) 


Although the notation of Van der Pol and Bremmer 
should have a historical priority, the Russian notation 
will be retained in the remainder of this Section. When 
q=0, the roots are determined simply from 


Wy (i) = 0 
[or from H»/3 [2(—2r,)3/2] = 0]. 
These particular roots designated ¢,° are well known.* 
The first five are given by 
£:°)==| 1,01879es" t4° 
to9 == 3:24820e—*"* [3° = 
ts? = 4.82010¢e—**/8 


I 


6.16331¢~ 7/8 
(CEN CAN ert 


A series expansion for f, in terms of ascending powers 
of g can be obtained readily. For example, since 


d ae Wilt | 
digi ais te ) 
dt dt 
Wi'() — — Wil) + WY) —=0 (61) 
dq dq 
and 


W(t) = —- tW(t), Wy (t) = qWi(t) 
it is seen that the differential equation for ¢, is 


iyi ¢ ae 
The resulting expansion is 
te = t,° + ae q = 1 Ge a (—.+ Ve (63) 
te 2A)? Shoe NOR 


+ higher terms in q. 

Other expansions for ¢, and a very extensive tabula- 
tion is given in a monograph by Belkina.” The results 
given here should be adequate for most conditions in 
practice. There is some difficulty however, when X be- 
comes small; this is near the shadow boundary. Here 
one must resort to a numerical integration of the inte- 
gral for g(X) since the residues converge poorly for small 
positive X and not at all for negative X. 
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Nonresonant Slotted Arrays* 
ANDRE DION} 


Summary—The distribution of slot conductance of nonresonant 
arrays is obtained by considering the array as a continuous line 
source. Distributions of conductance per unit length for three Taylor 
aperture distributions are thus obtained. However, the discreteness 
of the array is retained for a discussion of second-order beams and 
for the development of a method leading to their suppression. The 
performance of an experimental array is described. 


INTRODUCTION 
BT abo DESIGN of linear arrays having lengths of 


about ten wavelengths or more is generally facili- 

tated by considering the array as a continuous 
line radiator. Also, provided the array is not too long, 
it is permissible to neglect the effects of small losses oc- 
curring inside the feeding arrangement, which further 
simplifies the design. However, because the array is 
made of discrete elements, it is possible that more than 
a single principal maximum will appear in the radiated 
field ;! therefore it is necessary to retain the discrete ele- 
ments representation of the array at least for a discus- 
sion on secondary maxima. Slotted arrays of lengths 
conforming to the above description are considered here. 
Design data for these arrays are presented in conven- 
ient forms and an experimental array based on them 
agrees well with expectation. 


SECOND-ORDER BEAMS OF LINEAR ARRAYS 


The radiation pattern of linear arrays of discrete non- 
interacting elements may be expressed analytically as 
the product of two factors. One, the element factor, is 
the radiation pattern of a single element and generally 
has little directivity. The other, the space factor, is the 
radiation pattern of a similar array of point sources and 
is the factor controlling the directivity of the antenna. 
Referring to Fig. 1, where a linear array of NV identical 
radiators with spacing d is considered, the product of 
these two factors is? 


(G0) D Ca exp [in (s — = a cos s)] 


n=0 


with F(¢, @) as the element factor and the summation 
as the space factor. In this expression, C, is the excita- 
tion coefficient of the uth radiator, 6 is the phase dif- 
ference between the exciting field of adjacent radiators, 


* Manuscript received by the PGAP, February 9, 1958; revised 
manuscript received, June 23, 1958. 

+ Canadian Armament Res. and Dev. Establ., Valcartier, Quebec, 

an. 

+ D. W. Fry and F. K. Goward, “Aerials for Centimeter Wave- 
ee Cambridge Univ. Press, Cambridge, Eng., pp. 101-107; 

2 See, for example, S. Silver, “Microwave Antenna Theory and 
Design,” M.I.T. Rad. Lab. Ser., McGraw-Hill Book Co., Inc., New 
York, N. Y., vol. 12, ch. 9; 1949, 
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i 
and ) is the free-space wavelength. The factor (6 — (27/A)d_ 
cos ¢) in (1) is readily seen to be the phase difference ) 
between the far-field signals of two adjacent elements. 
The space factor is maximum when this difference cor- 


responds to an integral number of wavelengths, 7.e., 


2a 
MRR ie op a (2) 


where m=0, +1, +2, - - - and @,, is the direction of the | 
corresponding field maximum. Rearranging the terms, | 
this direction is obtained from 


Phere “(+ e m). (3) 


There may be more than one value of m yielding real 
values of ¢m, the amplitude of each of the resulting 
beams in the radiation pattern being proportional to 
the element factor. In general, the unwanted beams, 
called second-order beams, are appreciably reduced by 
this factor. However, in some designs the second-order 
beams may still be many times larger than the specified 
sidelobe level and therefore a method of eliminating or 
reducing them is required. Gruenberg? has proposed 
some methods of suppressing second-order beams, but 
in the following analysis the simpler method of limiting 
slot spacing is used and is shown to be practicable. 


Fig. 1—Geometry of array of N identical elements. 


In a slotted array antenna it is generally required to 
have the beam directed in some specified direction with 
respect to the array axis. This requirement is met bya 
proper choice of 6, d, and m in (3). The direction of the 
secondary beams nearest to the principal one is then 
obtained by replacing m by m+1 and m—1 in the same 
equation. Then 


$ nN E 1 r 
COS Oia = ee = eae 
eae) re (m+ | COS dm ‘ (4) 


3 H. Gruenberg, “Second-order beams of slotted waveguide ar- 
rays,” Can. J. Phys., vol. 31, pp. 55-69; January, 1953. 
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and 


ATS r 
COS $m—1 = +|+ Set 1)| = COS dm + 3 (5) 
where $m 41 and m1 are the directions of the second- 
order beams situated one on the positive side, the other 
on the negative side of the principal beam, respectively. 
It is seen that the modulus of (4) and (5) may be larger 
than unity provided )/d is large enough, and then the 
second-order beams are suppressed as they are in the 
imaginary region. For arrays made of elements radiating 
zero energy in the direction of the array axis, second- 
order beams can be avoided by choosing the slot spacing 
so that the second-order beam closest to the principal 
beam will be directed along this axis. The spacing do 
_ corresponding to this condition is 


r 


er nee 
; 1+ | cos om | 


(6) 


- and is seen to be dependent on the direction of the 
principal beam. This dependence is illustrated in Fig. 2. 
Limiting the slot spacing to this chosen value will be a 
practical method of suppressing second-order beams, 
_ provided that the principal beam direction as given by 
(3) is practically obtainable by a proper choice of 6. 
This is generally the case, as the following discussion 
shows. 


do/n 
© 


o 
° 
AS) 


END-FIRE 


SPACING IN WAVELENTGHS do 


| 
) 
0 cme SO CHO 50 *eCOme 70 “80° ~90 


180 170 160 150 140 -130 120 110 100 90 


DIRECTION OF BEAM er (DEGREES ) 


Fig. 2—Spacing between elements of a linear array 
vs principal beam direction. 


In a slotted array fed from one end (negative z axis 
in the present analysis) and terminated by a matched 
load at the other end, the phase difference 6 between 
the exciting field of adjacent slots is 


2nd 


(7) 


g 


where X, is the guide wavelength. 

Subjecting d to the relation expressed by (6), the re- 
- quired value of 6 may then be obtained by a proper 
choice of the guide wavelength. Substitution of (6) and 
- (7) in (3) yields the following expressions for \/Ag: 


A/Xq = COS bm + m(1 + COS $m) (8) 
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for 0 degrees<$n<90 degrees, i.e., for a beam tilted 
toward the load end of the array (forward-looking beam) 
and 


A/Xg = COS bm + m(1 — cos bm) (9) 


for 90 degrees <dm<180 degrees, i.e., for a beam tilted 
toward the generator end of the array (backward-look- 
ing beam). 

The relations expressed by (8) and (9) are plotted in 
Fig. 3. These plots give the relation that must exist 
between \/A,, m, and ¢@m in order that the nearest sec- 
ond order beam be directed along the 0-degree axis for 
a beam tilted toward the generator end, or along the 
180-degree axis for a beam tilted toward the load end 
of the array. 


BEAM TILTED TOWARD BBAM TILTED TOWARD 
LOAD GENERATOR 
1 


i 


0) 
Go © 15= SON esos 75, OR S05 
BROADSIDE 


Fig. 3—The relation between X/d,, the beam order m and 
the direction $m of the beam. 


It is noted that for m=0 the method is useful for 
forward-looking beams only. For m>0 it can be used 
for all directions of the beam. The value m=1 has a 
peculiarity in that for a backward-looking beam the 
required ratio \/X, is equal to unity for all values of $n. 
In effect for these particular values of m and )/), the 
second-order beam is directed along the 0-degree axis 
independently of \/d. This can be seen by substituting 
in (5) the value of 6 from (7) to give 

\ : 
COS dm-1 = — — — (m — 1) = 1. 
Xo 

Practical values for m are 0 and 1. It is not advisable 
to use values of m much larger than 1 as loading of the 
guide with a high dielectric constant material would be 
required to produce the extra phase shift in the line 
and such material presents losses that are generally 
objectionable. 

For a forward-looking beam, the choice between m =0 
‘and m=1 depends on practical considerations. For in- 
stance, for a near-broadside beam, taking m=0 would 
be impractical as it corresponds to large values of guide 
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wavelength where sensitivity to frequency shift is high. 
The choice of m=1 is then logical. 

The required value of \/d, is obtained by a proper 
choice of the guide dimensions and of the dielectric 
constant of the material filling the guide. For the par- 
ticular case of rectangular waveguides operated in the 
dominant mode, the ratio of the wavelength in the un- 
bounded medium of dielectric constant € to the wave- 
length in the guide is normally near 0.7 and therefore 
the ratio of free-space wavelength to guide wavelength 
is 


N/dy = 0.7+/¢, 


a relation from which the dielectric constant is deter- 
mined. The availability of low-loss dielectric materials 
of adjustable dielectric constant makes the above 
method of eliminating second-order beams practicable, 
at least for arrays of lengths such that dielectric losses 
are not prohibitive. 


(10) 


AMPLITUDE DISTRIBUTION 


The shaping of the radiation pattern of a linear array 
is generally done by exciting the elements of the array 
in a distribution such as Taylor or Dolph-Tchebycheff. 
A Taylor‘ distribution, which implies that the array of 
discrete elements is considered as a line source, was 
- selected for this investigation. This approximation is 
well justified in arrays made of about twenty elements 
or more. Eq. (1) for the radiation pattern of linear ar- 
rays is easily converted into the form studied by Taylor. 
Consider the array of Fig. 1 to be replaced by a line 
source of length 2/ and having a similar but continuous 
phase and amplitude distribution. The amplitude of the 
field produced at a distant point in a direction ¢, 0, by 
a distribution f(x) of the amplitude along such a line 
source, where x is the distance from the center of the 
line source, is proportional to 


u j24rx 
F(@, 8) [ Sta) exp = eics ae | dx (11) 


where ¢m is the direction of the beam maximum ob- 
tained from (2) and F(q, 9) is the element factor. 
Substituting 


al 
“= a (cos dm — cos ¢) (12) 
and 
1x 
C= ie (13) 


yields for the space factor an expression proportional to 
f soemap 


_4T. T. Taylor, “Design of line-source antennas for narrow beam- 
widths and low sidelobes,” IRE Trans. oN ANTENNAS AND PRopa- 
GATION, vol. AP-3, pp. 16-28; January, 1955. 
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which is the form studied by Taylor. The Taylor am- 
plitude distributions of line sources have been computed 
for space factors of 20, 25, and 30 db sidelobe ratio and 
are: 


(p)215 = 1+ 0.260 cos p + 0.043 cos 2p . 
— 0.049 cos 3p + 0.026 cos 4p (14a) 
f(p)2san = 1+ 0.443 cos p — 0.011 cos 2p : 
— 0.013 cos3p+0.010cos4p  (14b) 
f(p)s0an = 1 + 0.581 cos p — 0.030 cos 2p : 


+ 0.003 cos 3p + 0.002 cos 49. 


The parameter % of Taylor’s paper, which determines 
the boundary of the region of uniform sidelobes, was 
taken as 5 for reasons given in his paper. This parameter 
also determines the term at which the Fourier series — 
representing the amplitude distribution f(p) terminates. 


BEAMWIDTH OF NONRESONANT ARRAYS 


(14c) 


The beamwidth of a line source antenna of optimum _ 


design and with broadside radiation (resonant array) 
was shown by Taylor to be a function of the sidelobe 
ratio and of the parameter #7. The beamwidth of non- 
broadside arrays, however, also depends on the squint 
angle of the beam. As the beam direction departs from 
the normal to the array, its beamwidth increases. The 
amount of broadening may be obtained from (12). The 
space factor [", f(p)e*”"dp is maximum for u equal to 
zero or 6=¢m and symmetrical about the u~=0 axis. 
Let uo be the value of u for which the space factor is 
halved and £ the corresponding beamwidth; then 


21 
uy = =e | co om — COS (+n == 


e2 | 8 
Uy = ‘ co On Cos (on = =)| 


and for small values of the beamwidth the above ex- 
pressions yields 


and 


2 ALO v Bo 


~ Tsin dm 7. (15) 


sin dm 


where (8 is the beamwidth of a similar broadside array. 
The ratio of beamwidth 6 of nonresonant arrays to 

the beamwidth 6» of equivalent broadside arrays as a 

function of the squint angle, ¢n, is shown in Fig. 4. 


CONDUCTANCE OF ELEMENTS OF Non- 
RESONANT ARRAYS 


The computation of the conductance of the slots of a 
waveguide array to produce a given aperture distribu- 
tion is a tedious task when performed for each slot in 
succession. Considerable effort is saved if the array is 
considered again as a continuous radiator and the con- 
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Fig. 4—Beamwidth of nonresonant arrays as a function 
of the direction of the beam. 


ductance per unit length calculated. This is permissible 


for nonresonant arrays made up of elements having 


small or no mutual interaction (e.g., longitudinal shunt 


slots) and comprising a sufficient number of slots. The 
conductance of each slot is then small and each one 
produces only negligible reflection. When such an array 


is terminated in a matched load, it can be considered 
_ to be matched throughout its length; under this condi- 


tion, the conductance may be found from the relation® 
P(x) 
12) o(x) 


g(x) = (16) 


where P(x) is the power radiated per unit length, 


P(x) is the power in the guide, and g(x) is the normal- 
ized conductance per unit length at a distance x from 


- the center of the array. 


The power radiated per unit length is obtained from 
the desired amplitude distribution f(x) and may be 
written 


P(x) = K[f(«)? (17) 


where K is a constant. The power in the guide at a dis- 
tance « from the center of the array may be immediately 


_ written 


Po(x) =1—- if Ks) ae (18) 


where unit input power to the array is assumed and the 
guide is considered to be lossless. In order that the slot 
conductance be small, it is necessary that the power 
radiated by a given slot be small when compared to 


' the power in the guide at this slot position. This im- 


poses the condition that some power be absorbed in a 
load terminating the array. If 7 is the fraction of power 
absorbed by this load, the following relation may be 


written 


(19) 


1-r=f KL se) ae 
=} 


from which the proportionality constant K is evaluated. 
Substituting this value in (17) and (18) and using the 


6 Fry and Goward, op. cif., p. 113. 
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resulting expressions with 


in (16) yields the following expression for the normalized 
conductance per unit length times half the array length: 


wl fle)? 
1 a P 
ad leer - f" rores 


This result, it is recalled, applies only to cases where the 
guide is filled with nonlossy material and where wall 
losses are assumed to be negligible. Substitution in (20) 
of the expressions for f(p) given by (14) yields, after in- 
tegration, expressions for the normalized conductance 
per unit length of slotted waveguide arrays of optimum 
design with sidelobe ratios of 20, 25, or 30 db. The ex- 
pressions thus obtained for the conductance per unit 
length multiplied by half the array length are plotted in 
Figs. 5-7 for various percentages of power absorbed in 
the terminating load. 

The normalized conductance of a given slot of the 
array is obtained by multiplying the normalized con- 
ductance per unit length at that slot position by the 
spacing do, between two consecutive slots, or 


digray 
slot normalized conductance = /g(p) ave = ist?) (21) 


y N/2 


lg(p) = 


(20) 


where JV is the total number of slots. Therefore, to ob- 
tain the normalized conductance of each slot of the 
array from the graphs of Figs. 5-7, it is only necessary 
to divide the abscissa between —7z and 7 into a number 
of spacings equal to the number of slots minus one and 
to read the corresponding ordinate divided by a number 
equal to half the total number of slots. 


VOLTAGE STANDING WAVE RATIO OF 
TERMINATING LOAD 


It is necessary, for reasons expressed above, to ter- 
minate a nonresonant array into a matched load. When 
this load is not perfectly matched, a reflected wave prop- 
agates in the reverse direction in the waveguide and is 
partly radiated outside, thus interfering with the field 
produced by the direct wave and deteriorating the pat- 
tern. Because of the reversal of the direction of propa- 
gation, the reflected wave produces a radiation pattern 
with a beam pointing in a direction opposite (with re- 
spect to the array normal) to the design beam direction 
but tilted by the same amount. In this direction the 
amount of interference is largest and the total field is 
the phasor addition of the sidelobes generated by the 
direct wave and of the beam generated by the reflected 
wave. An estimate of the maximum VSWR allowable 
may be obtained by calculating the load VSWR for 
which the beam produced by the reflected wave is equal 
in amplitude to the design sidelobe level. The result of 
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VSWR is not too critical. 


EXPERIMENTAL RESULTS 
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A slotted waveguide array was constructed to deter- 
mine the errors introduced by considering the array of © 
discrete elements as a line source. The design param- 
eters chosen for this array are: sidelobe ratio 30 db, 
beamwidth 6 degrees, tilt angle 30 degrees, 1.¢., dm=60 — 
degrees, efficiency 65 per cent, frequency X band. 

The array was made of resonant slots cut longitudi- 
nally in the wide face of a guide, the spacing between 
the slots obtained from (5) or Fig. 2 being 2/3. For this 
value of spacing, the second-order beam is directed 
along the array axis and its amplitude is reduced to zero © 
since the radiation of a longitudinal slot in this direc- 
tion is zero. However, because of the relatively large 
beamwidth of the second-order beam in this direction — 


POWER 10 


: ™% AE +t 2 a "and because the field of a slot is still appreciable at 
Fig. 6—Conductance distribution of elements of a nonresonant angles close to the array end) it 1s necessary to make the 
slotted array designed for 25-db sidelobe ratio. spacing between slots a little less than the value indi- 


cated above. Thus the second-order beam may be 
thrown further into the imaginary region so that no 
part of it greater than the design sidelobe level remains 
in the visible range. The amount by which the spacing 
must be decreased depends on the sidelobe ratio, on 
the beamwidth, and on the radiation pattern of a single 
element, and is generally small in most cases. It is not 
advisable to make the slot spacing much less than the 
value obtained from these considerations because of the 
increased amount of interaction which results when the 
slots are too close. Higher order waves inside the guide 
are excited by the slots. Increased interaction results 
when the higher order waves are not sufficiently attenu- 
ated in the distance between two consecutive slots. 

The tilt angle of the beam determines the ratio of 
free-space wavelength to guide wavelength. From Fig. 
3 it is seen that a value of 3 for this ratio satisfies the 


Fig. 7—Conductance distribution of elements of a nonresonant requirement. For this particular value the required 
slotted array designed for 30-db sidelobe ratio. guide wavelength was easily obtained by a suitable 
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choice of the guide dimensions at the design frequency. 
No dielectric loading was required. 
The design beamwidth was found to be realizable 
with an array of 20 slots. The displacement y of a slot 
from the center of the guide was obtained from Steven- 
son’s formula,®:7 


an aN 
g = 2.09 — cos? (=) sin? (=) 
DN 2Xq a 


where a and 0 are the wide and narrow dimensions of 
the guide and the conductance g is obtained from the 
proper curve of Fig. 7. The distance between two con- 
secutive slots was too small to permit introducing an- 
other slot midway between them and displaced on the 
other side of the center line, and therefore the array was 
single-row. 

The experimental H-plane pattern of the antenna is 
_ presented in Fig. 9 where it is compared with the theo- 
_ retical pattern. This pattern was taken in the longitudi- 
nal plane perpendicular to the guide surface (1.e., 9=0). 
Measurements taken in planes inclined with respect to 
this last plane were also taken and did not reveal the 
presence of second-order beams. Because the radiation 
is omnidirectional in one plane, the accurate measure- 
ment of small sidelobes was awkward due to the diffi- 
culty of eliminating small reflections from the ground 
and surrounding objects. These reflections may be re- 
sponsible for the small discrepancies observed close to 
the beam. At large angles from the beam, agreement is 
not to be expected because the theoretical pattern is 
derived on the basis of a continuous line source, a deri- 
vation which does not reveal the effect of second-order 
beams and therefore cannot predict the behavior at large 
angles. The proper matching of the terminating load is 
indicated by the absence of an increase of radiation at an 


(22) 


6 A. F. Stevenson, “Theory of slots in rectangular waveguides,” 
| J. Appl. Phys., vol. 19, pp. 24-38; January, 1948. a 

7A. A. Oliner, “The impedance properties of narrow radiating 
slots in the broad face of rectangular waveguide,” parts I and II, 
IRE TRANS. ON ANTENNAS AND PROPAGATION, vol. AP-5, pp. 4-20; 
January, 1957. 
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Fig. 9—Principal H-plane pattern of experimental array. 


angle of 120 degrees which is the direction where the 
beam produced by a load reflection should be found. 


CONCLUSION 


Data useful to the design of nonresonant: arrays have 
been developed and presented in convenient form. It 
has also been shown that limiting the slot spacing is an 
appropriate method for suppressing second-order beams. 
In cases where dielectric loading of the guide is needed, 
the method is still practicable for arrays of lengths in 
which dielectric losses are not prohibitive. 
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Gains of Finite-Size Corner-Reflector Antennas* 
H. V. COTTONY} anp A. C. WILSONT 


Summary—An experimental corner reflector was erected at the 
Table Mesa antenna range near Boulder. The aperture angle of this 
antenna was made adjustable to any value between 20 and 180 de- 
grees. The widths and lengths of the reflecting surfaces were each 
adjustable from 0.4 to 5.0 wavelengths. Measurements of gain 
were made for numerous combinations of lengths and widths of re- 
flecting surfaces. These measurements were made with a half-wave 
dipole in the first, second and third maximum positions. The aper- 
ture angle was adjusted to maximize the gain. The principal results 
are presented in the form of contours of constant gain plotted for a 
range of widths and lengths of reflecting surfaces from 0.4 to 5.0 
wavelengths. These graphs should be useful to a designer of corner- 
reflector antennas. 


INTRODUCTION 
\ Y ARIOUS aspects and applications of corner-re- 


flector antennas have been described in numer- 

ous papers. The literature! presents the more 
important features of this type of antenna. While the 
volume of technical literature dealing with this antenna 
is extensive, a systematic and detailed investigation of 
the effect of the lengths and widths of reflecting surfaces 
on the gain realizable from a corner-reflector antenna 
has not been reported. Neither is there any record of 
experimental measurements on corner-reflector antennas 
having values of aperture angle other than those ob- 
tained by dividing 180 degrees by an integer, 7.e., 90, 
60, 45, etc., degrees. 

In 1952 and 1953 measurements of gain of various 
corner-reflector antennas were carried out at the NBS 
model antenna range at Sterling, Va. The purpose of 
these measurements was to aid in selecting the size of 
reflecting surfaces and the value of the aperture angle 
for corner-reflector antennas to be used with some of the 
VHF ionospheric scatter circuits, then being installed 
for the U. S. Air Force. Because the corner-reflector 
antenna continues to be a useful type for this applica- 
tion, additional measurements were carried out in 1956. 
These experiments combine the study of the size of 
reflecting surfaces with the value of aperture angle 
required to maximize the gain. 

The term “gain” as employed here denotes power 
gain as defined by IRE (1948). However, a half-wave 
dipole rather than an isotropic radiator is used as the 
standard. 


* Manuscript received by the PGAP, March 26, 1958. This work 
has been carried out on behalf of and has been supported by the 
U. S, Air Force Communication System, Strategic Communication 
System Engineering Organization. The results represent a part of the 
supporting work of the Natl. Bur. of Standards in connection with 
establishment of the U. S. Air Force transatlantic VHF ionospheric 
scatter system. 

+ Natl. Bur. of Standards, Boulder, Colo. 

+ J. D. Kraus, “The corner-reflector antenna,” Proc, IRE, vol. 
28, pp. 513-519; November, 1940. 

2 E. B. Moullin, “Radio Aerials,” Oxford University Press, Lon- 
don, Eng.; 1949. 

* J. D. Kraus, “Antennas,” McGraw-Hill Book Co., 


Inc., New 
York, N. Y., pp. 324-352; 1950. 
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EXPERIMENTAL PROCEDURE 


The measurements were carried out at 400 mc. The © 
corner reflector consisted of two lattice-type wooden 
frames, each five wavelengths (12.28 feet) wide by five © 
wavelengths long, supporting the reflecting surfaces. © 


The reflecting surfaces were made of overlapping strips 
of sheet aluminum fastened to these frames in such a 
manner that each strip added 0.2 wavelength to the 
length of the reflecting surface. Fig. 1 illustrates the 
terminology. 


DISTANCE OF DIPOLE 
FROM APEX 


Fig. 1—Diagram of a corner-reflector antenna showing the 
nomenclature and the symbols used in the text. 


Fig. 2 is a view of the experimental corner-reflector 
antenna employed in these measurements. The lengths 
of the reflecting surfaces were varied by adding or re- 
moving the strips. The width was varied by trimming 
the lengths of the strips. The two frames were pivoted 
along the same axis and the aperture angle was adjust- 
able from 20 to 180 degrees. 

The driven element was a folded half-wave dipole. 
Its position within the corner, S, measured from the 
apex could be varied from approximately 0.07 to 2.5 
wavelengths. The support for the dipole served also as 
a balun transformer. It was designed to transform the 
impedance of the dipole to approximately 50 ohms. A 
two-stub tuning unit was employed to attain an exact 
match to the 50-ohm line. The gain of the corner-re- 
flector antenna was measured by comparing its response 
with that of an antenna, the gain of which was ac- 
curately known. Corrections for variations in system 
sensitivity with time were made by comparison with a 
nearby reference antenna. The details of the method of 
gain measurement were described in a recent paper. 

The bulk of the measurements consisted of deter- 
minations of gains realizable with the reflecting surfaces 
of given widths and lengths, with the dipole in each of 
the first three positions which gave maximum response 
and with the aperture angle adjusted to maximize the 


rs; iH. V. Cottony, “Method for accurate measurement of antenna 
gain,” presented at URSI Meeting, Washington, D. C.; May 1, 1956. 
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Fig. 2—A view of the experimental corner-reflector antenna. 


gains for each position. Additional measurements were 
carried out to trace the changes in the locations of the 
three dipole “maximum positions” with changes in the 
aperture angle. Some measurements were carried out to 
determine the effect of the aperture angle on the gain, 
with the size of the reflecting surfaces held constant. 


RESULTS 


The principal results of this series of measurements 
are represented by the data on the gains realizable from 
corner reflectors of various lengths and widths. This 
material is presented in a manner which should permit 
a designer to select the size of corner reflector best 
- suited to his application, considering both performance 
and economy in size. The work also led to a better un- 
derstanding of the various “maximum positions” of the 
dipole and of the changes in gain with continuous varia- 
tions of aperture angle. 


The Effect of Aperture Angle on the Optimum Position 
of the Dipole 

As is well known, in order to get maximum gain in the 
forward direction from a corner-reflector antenna hav- 
ing a given angle of aperture, the dipole may be at one 
of several specific distances from the apex of the corner. 
These locations are usually referred to as “first posi- 
tion,” “second position,” “third position,” etc. The first 
position is relatively broad and involves a compromise. 
The directivity actually improves as the dipole is 
brought very close to the apex; however, the radiation 
resistance becomes quite low and the losses reduce the 
gain. If the position of the dipole is displaced from the 
optimum in the direction away from the apex, the gain 
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is reduced first by the broadening and, when the dis- 
placement is further increased, by the splitting of the 
main beam. The second and higher order maxima corre- 
spond to points at which the direct radiation and re- 
flected radiation add most favorably in the forward 
direction. These points are very much more sharply 
defined than is the location of the first position. 

The earlier of these series of measurements were 
carried out using reflecting surfaces 5 wavelengths wide 
by 5 wavelengths long, with the purpose of tracing the 
location of the first, second, and third maximum posi- 
tions with varying aperture angles. However, it was 
found that the dimensions of the reflecting surfaces 
did not materially affect the optimum position of the 
dipole. Fig. 3 presents the results of these tests. The 
first position is found to be rather broad and poorly de- 
fined with the minimum value of 0.10 to 0.15 for an 
angle of 180 degrees. It is a continuous function, the 
distance from the axis increasing as the angle of aper- 
ture is made smaller. 
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APERTURE ANGLE DEGREES 


Fig. 3—Locations of the three dipole positions as functions of aper- 
ture angle. A) Point of greatest gain for second position (¢=65°; 
Genlaili| ne B) Point of greatest gain for third position (¢=90°; 
S=1.50 A). 


The second position is found to be a discontinuous 
function. Its position roughly parallels the first position, 
but for aperture angles between 120 and 100 degrees 
the second position moves out to meet the third posi- 
tion. The second position is found to reappear again for 
an aperture angle of approximately 75 degrees and 
continue to the smallest angle méasured (20 degrees). 

For an aperture angle of 180 degrees, the third posi- 
tion is located at 1.25 wavelengths. The variation in its 
position with changes in aperture angle is somewhat 
irregular, being affected by the second position which 
joins it at a point corresponding to an angle of aperture 
equal to 100 degrees. The third position is found to 
disappear when the aperture angle is decreased below 
70 degrees. It cannot be traced beyond this point be- 
cause the limit of available dipole movement, 0-2.5 
wavelengths, was insufficient to reach it. (At an aper- 
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Fig. 4—Contours of constant gain for corner-reflector antenna with 
dipole in first position for various reflector sizes. The aperture 
angle adjusted for maximum gain in each case. The optimum 
value of aperture angle, in general, diminishes with increase in 
the length of reflecting surfaces. 
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Fig. 5—Contours of aperture angle required to optimize the gain of 
a corner-reflector antenna, with dipole in first position, for various 
reflector sizes. 


ture angle of 60 degrees the third position is theoreti- 
cally predicted to be 2.6 wavelengths away from vertex.) 

With the dipole in the first position, the aperture 
angle required to maximize the gain is a function of 
the length of the reflecting planes, the larger the planes 
the smaller the optimum angle and the greater the 
realized gain. With the dipole in the second position, on 
the other hand, the maximum gain is realized when the 
angle of aperture is adjusted to 65 degrees regardless of 
the size of the reflecting planes, provided these are larger 
than about 1 wavelength in length and width. The maxi- 
mum is well defined and a change of aperture angle in 
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Fig. 6—Contours of constant gain for corner-reflector antenna with 
dipole in second position for various reflector sizes. The aperture 
angle adjusted for maximum gain in each case. For lengths greater 
than 1 wavelength, the optimum value of aperture angle has been 
found to be 65 degrees. 
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Fig. 7—Contours of constant gain for corner-reflector antenna with 
dipole in third position for various reflector sizes. The aperture 
angle adjusted for maximum gain in each case. For lengths greater 
than about 0.5 wavelength, the optimum value of aperture angle 
has been found to be 90 degrees. : 


either direction by 5 degrees results in a reduction in 
gain of the order of 0.3 db. With the dipole in the third 
position there is likewise one optimum value of aperture 
angle for which the gain is maximized. This value is 
found to be 90 degrees and is valid for all sizes of reflect- 
ing planes longer than 0.5 wavelength. 

It should be noted that with moderate-size corner 
reflectors, the relationship between gain and the size 
of the reflecting planes and that between gain and the 
angle of aperture, although continuous, is not smooth, 
but contains numerous small-scale (of the order of 0.1 
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db) real variations superimposed upon the major varia- 
_tion. The results presented in this report have been 
smoothed to remove these small-scale variations. 


The Effect of Dimensions of the Reflecting Surfaces on the 
Gain 

The gains realizable with a corner-reflector antenna 
having various combinations of length and width of re- 
flecting surfaces are presented in Figs. 4, 6, and 7. Fig. 
4 represents results obtained with the dipole in the first 
position. Since the optimum aperture angle in this case 
depends on the length and, to a lesser extent, on the 
_ width of the reflecting surfaces, the optimum values of 

aperture angle are presented in Fig. 5. Figs. 4 and 5 
should be used together, 7.e., to realize the gain of a 
_corner-reflector antenna with reflecting planes having 
length, Z, and width, W, which is given by Fig. 4, the 
‘aperture angle must be adjusted to values given in Fig. 
5 for corresponding length and width. The location of 
the dipole, in its turn, should be selected from Fig. 3 

corresponding to the aperture angle which, in turn, is 
- obtained from Fig. 5. 

The gains realizable using a corner-reflector antenna 
with its dipole in the second position are presented in 
Fig. 6. For all sizes of reflecting surfaces greater than 
approximately 1 wavelength long by 1 wavelength 
wide, the maximum gain was realized with the aperture 
angle optimized at 65 degrees, and with the dipole 1.17 
wavelengths from the apex. For very small corner re- 
-flectors the optimized value of aperture angle was irreg- 
ular and because of the relative unimportance of such 

corners no attempt is made to present the data on the 
exact values of the aperture angle for small corner-re- 

flector antennas with dipole in the second position. 

Fig. 7 presents the gains of the corner-reflector an- 
tenna with the dipole in the third position. As in the 
} preceding case, the aperture angle optimizes at very 
nearly the same value for all combinations of length 
and width greater than about 0.5 wavelength. In the 
case of the dipole in the third position, the optimum 
value was found to be 90 degrees with the dipole posi- 
Y tion 1.5 wavelengths away from the corner. 
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CONCLUSIONS 


To be of practical use, a corner reflector should have 
a width of at least 0.5 wavelength. Considerable im- 
provement in gain is realized if the width is increased 
to a full wavelength. Beyond that, the increase in gain 
with increase in width is much slower, and beyond a 
certain point an additional incremental increase in 
width may actually result in a decrease in gain. The 
phenomenon can be understood by applying the con- 
cept of Fresnel zones on the reflecting planes. 

In the case of the dipole in the first position the rela- 
tionship between the length of the surfaces and the gain 
is simple: The gain increases monotonically with in- 
crease in length. In the cases of the dipole in the second 
and third positions, the same relationship holds true 
generally except for lengths of less than 1.5 wave- 
lengths. For these shorter lengths the relationship is 
irregular. Using the dipole in the first position, the 
optimum aperture angle depends on the size of the cor- 
ner, particularly its length. Using the dipole in the sec- 
ond position, the optimum aperture angle is 65 degrees 
for reflectors of length greater than 1 wavelength. With 
the dipole in the third position the optimum aperture 
angle is 90 degrees for reflectors having length greater 
than approximately 0.6 wavelength. 

If the designer’s goal is maximum gain along the for- 
ward axis, there may be an advantage in employing a 
corner-reflector antenna with the dipole in the second 
position. To realize this advantage the lengths of the re- 
flecting surfaces should be greater than 2 wavelengths. 
The widths of the reflecting surfaces should be ap- 
proximately 2 wavelengths. For smaller corner reflec- 
tors there is an advantage from the standpoint of gain 
in the use of an antenna having its dipole in the first 
position. There do not appear to be any advantages in 
placing the dipole in the third position. 
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Scanning Surface Wave Antennas—Oblique Surface Waves 
Over a Corrugated Conductor* . 


R. W. HOUGARDY} anv R. C. HANSENT 


Summary—The existence of a surface wave which propagates 
across a corrugated metallic surface at an oblique angle with the 
teeth is investigated both experimentally and theoretically in this 
paper. Expressions are derived which give the variation of the wave 
velocity and amplitude with the change of wave direction. Experi- 
mentally measured values of the surface velocity compare favorably 
with the theory. 

The radiation pattern of an experimental antenna is given which 
demonstrates that a low sidelobe, narrow azimuth beam scannable 
to +30 degrees with a cosecant-squared elevation pattern is attain- 
able. A method of feeding this antenna to give a low silhouette, mak- 
ing the corrugated scanner antenna suitable for flush mounted 
applications is illustrated. 


INTRODUCTION 


URFACE wave antennas in their most common 
S form consist of a finite width and length of trap- 
ping surface, usually planar, which may or may not 
be imbedded in a ground plane. Some device, e.g., a 
waveguide horn, dipole array, or pillbox, is used to ex- 
cite a surface wave at one end of the trapping surface 
_ with the energy end firing off the opposite end. 
- Patterns produced by surface wave antennas are fan 
beams, with the pattern in the plane of the surface 
being approximately that of the feed. Gain may be 
maximized by adjusting the phase velocity of the sur- 
face wave along the trapping surface to give a net phase 
shift of 2.921 radians (the Hansen-Woodyard condi- 
tion) ; however, in this region the gain is changing slowly 
so a small deviation from the 2.921 figure will not re- 
duce the gain significantly. The gain is directly propor- 
tional to the length of the surface, and is G~4.83 1/). 
Beamwidth,! in a plane perpendicular to the surface, is 
approximately inversely proportional to square root of 
length, with typical minimum figures being 15 degrees 
for /=4) and 10 degrees for 10A. A finite length of travel- 
ing wave radiator either in free space or attached to a 
finite ground plane has the characteristic end fire pattern 
tilted above the surface plane, due to diffraction of the 
wave at the end of the surface; the maximum beam tilt 
is inversely proportional to square root of length (as is 
the beamwidth) and is given by O4:42%124°/4/Bl. 
Traveling-wave end fire antennas trade the height of 
an equivalent two-dimensional slot array or dish for 


* Manuscript received by the PGAP, August 2, 1957; revised 
manuscript received, June 5, 1958. This research was supported by 
the AF Cambridge Res. Center, Air Res. and Dev. Command, Cam- 
bridge, Mass., under Contract AF19(604)-1317, and is an extract 
from Hughes Aircraft Co., Culver City, Calif., Scientific Rep. 13 of 
the same title. 

{ Microwave Lab., Hughes Aircraft Co., Culver City, Calif. 

*R. S. Elliott, “On the theory of corrugated plane surfaces,” 
IRE Trans. ON ANTENNAS AND PRopaGarTIoNn, vol. AP-2, pp. 71-81; 
April, 1954. 
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length in the beam direction, and while the end fire an- — 
tennas are not capable of producing very-marrow eleva- _ 
tion beams, their low height or silhouette (often less — 
than two wavelengths) adapts them to flush mounting 
on air frames and missile bodies. The most often used — 
trapping surface is a corrugated metallic sheet where © 
the teeth can, if desired, be filled with a low loss ceramic — 
material. 

Two points justify the use of surface, or slow, waves. 
First, the use of a trapping surface allows the elevation | 
pattern to be shaped, with a cosecant square being © 
typical of what can be accomplished by slight modifica- _ 
tion (curving and tapering) of the trapping surface. — 
Second, the maximum end fire gain condition requires a _ 
slow wave with velocity given by v/c=1/(1+2.921/kl); | 
e.g., for a maximum gain antenna of length 10A, the — 
ratio should be 0.955. 


QUALITATIVE DISCUSSION OF THEORY 


Surface wave propagation over a corrugated metallic 
ground plane has been studied quite extensively during 
recent years;'~’ Elliott! has presented a list of earlier © 
papers. These authors have been concerned with the 
problem where the direction of wave propagation across 
the surface is normal to the corrugations. 

Present requirements for flush mounted antennas de- 
mand that it be possible to scan the beam of such an — 
antenna in the plane of the surface, which, in the case of 
a corrugated conductor antenna, means that the wave 
propagation direction will not always be normal to the 
corrugations, but must be along the direction of the scan 
angle. (See Fig. 1.) Questions have been raised about 
whether surface waves can exist when propagated at an 
angle across the corrugations. It was the original pur- 
pose of this work to demonstrate that such modes could 
satisfy the boundary conditions. In the last section the 
propagation of a skew surface-wave across a corrugated 
metallic surface is analyzed. The analysis indicates that 


*L. Brillouin, “Wave guides for slow waves,” J. Appl. Phys., 
vol. 19, pp. 1023-1041; November, 1948. 

*D. K. Reynolds and W. S. Lucke, “Ridge and Corrugated An- 
tenna Studies,” Second Quart. Prog. Rep., Stanford Res. Inst. 
Menlo Park, Calif., Project No. 199; January, 1950. : 

*D. K. Reynolds and W. S. Lucke, “Corrugated end fire an- 
ee Proc. Natl. Electronics Conf., vol. 6, pp. 16-28; September, 

®W. Rotman, “A study of single-surface corru ated gui Z 
Proc. IRE, vol. 39, pp. 952-959; August, 1951. ee 

6A. L. Cullen, “The excitation of plane surface waves,” IEE 
Monograph 93 R, vol. 101, part 4, pp. 225-234; February, 1954. 

R. A. Hurd, “The propagation of an electromagnetic wave along 


an infinite corrugated f. 3 - 
734; December, 1054. Mame Aah a8 
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Fig. 1—Oblique propagation over a corrugated surface. 


_as the direction of propagation is varied from normal to 
the corrugations (@=0 degrees) to the parallel direction 
_ (8@=90 degrees), the free space velocity to surface wave 
velocity ratio decreases to the plane wave value of unity. 
Fig. 2 depicts the propagation constant ratio or velocity 
ratio as a function of angle @ for two depths of teeth and 
for gap widths d/d of 0.10 and the limit value zero. For 
the low values (near unity) required for most practical 
end fire antennas, the velocity ratio changes slowly 
_ with angle as seen in the curves of Fig. 2. The limit 
curves for h/ equaling the cutoff value at@=0 degrees 
are also shown. Recent independent work by Kay has 
considered the efficiency of skew surface-wave excita- 
tion.® 


EXPERIMENTAL RESULTS 


Some experimental verification of the theory was at- 
tempted by making actual measurements of the wave 
velocity over a corrugated surface. The surface wave 
was excited using a linear waveguide slot array which 
operated at X-band frequencies. The corrugated surface 
which supported the surface wave has as its geometry 
d/X=0.0953 and h/A=0.0648 where d and h are the gap 
width and tooth depth, respectively. The wavelength 
of the surface-supported wave was measured for scan 
~ angles of from 0 to 30 degrees (in 5-degree increments). 
Values of c/v calculated from these data are plotted 
- with the corresponding theoretical values in Fig. 3. 
Since the wave velocity ratio does change slowly 
- and continuously for values of c/v near unity as the 
wave propagation direction changes, it was desirable to 
determine what effect the scanning of a practical corru- 
gated surface antenna had on the radiation patterns. 
To this end a corrugated surface was constructed with 
a design very similar to that of an existing search radar 
_ antenna which, however, cannot be scanned. Both an- 
_tennas were constructed to give approximately a cose- 


8 A. F. Kay, “Excitation Efficiency of Surface Waves Over Cor- 
rugated Metal and Doubly Corrugated Metal and in Dielectric 
Slabs on a Ground Plane,” Sci. Rep. No. 5, Tech. Res. Group, 
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Fig. 2—Velocity ratio vs scan angle: limiting cases. 
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Fig. 3—Experimental velocity ratios. 


cant squared elevation pattern. The experimental sur- 
face was fed by a linear waveguide slot array (Fig. 4) 
designed for beamwidth of 9 degrees and 25-db side- 
lobes, and scanned by insertion of spacers in the snake. 
The azimuth patterns for this antenna are shown in 
Fig. 5 for scan angles of 0, 8, 15, and 25 degrees, and it 
is apparent that the main beam remains unchanged as 
the surface is scanned. Sidelobes are below the recorder 
threshold. The secondary beam which appears for larger 
scan angles is a result of slot spacing of the exciting 
source and, of course, is not a consequence of scanning 


over the corrugated surface. 


372 


Fig. 4—Prototype corrugated scanner antenna. 


Of greater interest, however, are the elevation pat- 
terns of this antenna (Fig. 6). The shape of beam is not 
greatly changed as the antenna is scanned, and although 
no exhaustive effort was made to obtain an exact cose- 
cant squared beam, the elevation patterns remain within 
plus or minus 3 db of cosecant squared from 10 to 80 
degrees. 

For the experimental antenna, it was convenient to 
mount the snake feed waveguide slot array in a plane 
perpendicular to the plane of the corrugated surface, as 
shown in Fig. 4. However, the snake feed array or other 
feeding device may be placed parallel to and just under 
the corrugated surface with a small splash plate convey- 
ing the energy up to and over the edge of the corrugated 
surface. Such an antenna could be constructed with a 
low silhouette, as indicated by the laboratory model 
shown in the photograph of Fig. 7. Patterns of this an- 
tenna were taken and found to be substantially the same 
as those of the previous model. 


THEORETICAL ANALYSIS 


The method employed is that used by Hurd? on 
propagation normally across a corrugated surface, and 
includes 1) writing the complete fields both above the 
corrugations and in the gaps, 2) matching the tangential 
field components at the boundaries, and 3) using a cal- 
culus of residues technique to evaluate the resulting 
infinite set of simultaneous equations. To facilitate the 
solving of this problem, an idealized physical structure is 
assumed on which the corrugations or teeth are vanish- 
ingly thin. Another important assumption is that all 
modes but the dominant mode in the slot region are at- 
tenuated to negligible amplitude in traveling down the 
slot to the ground plane and back to the mouth of the 
slot. The conditions for this assumption to be good will 
be given later. The approach outlined above is not cog- 
nizant of finite tooth thickness, which is a second order 
effect, but does account for the tooth spacing which has 
an appreciable effect on the propagation constant. 

The physical structure is shown in Fig. 8. A conduct- 
ing plane of infinite length and width forms the base 
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Fig. 5—Azimuth patterns of corrug. 
angle=0 degrees. (b) Scan angle = 
degrees. (d) Scan angle =25 degrees. 


ated scanner antenna. (a) Scan 
8 degrees. (c) Scan angle =15 


for an infinite number of vanishingly thin, evenly 
spaced parallel metal vanes. The vanes are of height h 
and separated by a distance d. The origin of coordinates 
is at the top of a vane with the g axis lying along the 
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Fig. 6—Elevation patterns of corrugated scanner antenna. (a) @=0 
degrees. (b) 0=8 degrees. (c) 6=15 degrees. (d) 9=25 degrees. 


vane and parallel to the ground plane. Fields produced 
by a source of finite dimensions will not be considered; 
only the surface: wave mode on an infinite structure is 


assumed. 
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Fig. 8—Infinite corrugated plane surface. 


The allowable electromagnetic field above the corru- 
gations is characterized by the absence of an E, com- 
ponent, and thus the field can be expressed as 


E= curle,f — H = — curl curl of (1) 
JOM 
where f is a solution of the scalar wave equation and 
a, denotes a g-directed unit vector. 

To obtain the correct solution of the scalar wave 
equation for the region above the corrugations, the de- 
sired wave structure is examined. The fields are to be 
attenuated for increasing y, since only surface waves are 
to be considered. Due to the nonuniform surface struc- 
ture in x, higher order surface modes will appear in this 
direction; however, along the z axis a field variation 
with but a single propagation constant is valid, since at 
the y=0 plane the fields must match for all z. A time 
dependence of e*** and loss free media throughout is as- 
sumed. 

The solution of the scalar wave equation which satis- 
fies these conditions is 


fn = An exp (—JBn% — jtz — any) (2) 


where Bn, T, and a, are related to the free space propa- 
gation constant by 


kh? = By? + 7? — an’. (3) 


Since Floquet’s Theorem applies to the periodic paral- 
lel-plate configuration in x the propagation constant in 
that direction must take the form 


Bn = Bot 2nn/d (n=0, +1, 42,---) 


where Bo, which defines the fundamental mode, must be 
determined by the corrugation geometry. Thus it is 
seen that for the dominant mode, the direction of propa- 
gation is given (see Fig. 1) by 


6 = tan (r/Bo) (5) 
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where @ is measured from the x axis. The free space to 
surface wave velocity ratio for this mode is therefore 


Vo Bo ey? 
er iE. (5 
Vs P k 
The required solution of the scalar wave equation for 
the region between the corrugations is readily deter- 
mined from the boundary conditions imposed upon the 


field components at the bottom and sides of the gaps. 
Thus 


ip Bin cos x[cos ¥m(y -+ h) |e-#7 (7) 


where 


mr 


4 2 
a(R) pet taet 012-5). &) 


Equating tangential electric and magnetic components 
of the field at the plane y=0, three infinite sets of simul- 
taneous equations result, of which only two sets are in- 
dependent. The orthogonality properties of the trigono- 


metric functions allow these to be reduced to the form? 


d €mY¥mPm sin Ymiht = —JanBnAn 
Fy: ee oa ee Eh ame eae (9a) 
2 he Le exp (—jBod) Sail OMEN eS Ym 

and 
d €mBm COS Ymht Re lope 
mer : eB corey apr te 
2 (4) exp (—j@od) spat! n=—0 An se Ym 


where the Neumann number is 
A? m= 0 
En = 
1, m > 0 


Multiplying (9b) by jym and then adding and subtract- 
ing (9a), there results, after some simplification, 


—d €mVmDBme iim a Bn m 
aan 5 : ee (10) 
CS Wrexp(7eidea el no, Oh Fyn 
and 
d €nY¥mDmeryn" Es nA 
= > , = 2) ——. (11) 
2 G1) exp (—jBod) = 1 n=—«0 An aie IVm 


If d</2, then all modes except m=0 are evanescent. 
Furthermore, since ina practical system h> d, the factor 
exp (—JYmi) in (10) becomes negligible even for m=1. 
Thus, neglecting the left-hand side of that equation 
when m>0, 


od Be-770h Baan 
SO) aa ——————— 12 
exp (—jGod) =) 4! n=—«o An — RV ( ) 


Mes 


where 69m is the Kronecker Delta symbol. 


From the two infinite sets of simultaneous equations 
(11) and (12) must be found the amplitude coefficients 
A, and B, and the wave number Bo, the latter yielding 
directly the propagation constant v. 

These equations may be solved using a calculus of 
residues technique? which enables one to write the — 
right-hand side of (11) and (12) as the sum of the resi- | 
dues of a function f(w). This function f(w) is defined as i 


1) having simple poles at w=a,(n=0, See ses), i 
2) having simple zeroes at w=j7Yn(m=1, 2, 3,---), : 
3) tending to zero as | w — oo except at w=Qn. : 


By choosing as the path of integration contours of in- ; 
creasingly large radius which avoid the poles w=a, and | 
enclose w=7Yo | 


ae AA ne 3 eS + Somf(jro) = 0 (13) | 


LrjJeoW-fi¥n nb-a & — YY. 


where R(an) is the residue of f(w) at w=a,. Comparing — 
this equation with (12), the coefficients are obtained in 
terms of the function: | 


0d Boe 179% 


f(jyo) = — ay ar (14) 
and 
R(an) 
A, = . 15 
3 (15) 


With the same definition for f(w) 


1 f(w)dw be R(an) 
ae SS = he ee — (— Ym) — 0. 16 
29] c Ww ia 2) Qn + IVm I : ( ) 
When this equation is compared with (11), the same 
formula results for 4,, but also 


ds ém¥mBmetr™* 


t(—d¥m) = 9 (— 1) ™¢—ibod | : 


(17) 


Now if f(w) were known, the mode amplitudes could be 
found from (15) and (17). 

It is essential that the function has the prescribed 
poles and zeroes, and must satisfy (14). In addition to 
conditions 1), 2), and 3), the function f(w) must also 
satisfy the edge criterion which requires that the field 
at the edge of each corrugation be singular in the allow- 
able fashion," as r—/?, Thus it can be shown! that the 
unique expression for f(w) is 


*R. C. Hansen, “Single slab arbitrary polarization surface wave 
structure,” IRE TRANs. ON MICROWAVE THEORY AND TECHNIQUES, 
vol. MTT-5, pp. 115-120; April, 1957. 

°C. J. Bouwkamp, “Diffraction theory,” Reps. Prog. Phys., 
vol. 17, pp. 35-100; May, 1954. 
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—-yedBo exp (—joh) exp 
f(w) = 


~ where 


eet Gai erg exp( =) (19) 


p=1 pr 


Taw) = TI (wo — a5)(@ — a) (*) exp or (20) 


2pr 


p=1 


In these infinite products the exponential terms are 


_ necessary to produce convergence. 


The combination of (14) and (17) for m=0 yields the 
formula 


fio) = — f(—Jro) exp (— 2jyoh) (21) 


from which the phase constant for the dominant mode 


- may be evaluated. Since (21) must be satisfied in argu- 
_ ment as well as modulus, 


RS Geng! ee ee 
cos-! — + > {sin — — sin-! — — sin! 
Dev apat pn Bp 


vod In 2 
LAPlhtted 


Yo 
| 6» } 
= yoh. (22) 


T 


The inverse sines in the summations may be approxi- 
mated by their arguments, and so after some simplifica- 
tion the following expression is obtained which relates 


_ the propagation constant v to the scan angle 6 (since 


y?=k?+8,2—~.”) and to the physical dimensions d and h. 


d d 
Yoh = cost + *lv(1 7 =) 


0 Tv 


Bod 2 
+ x (cot = —) + 2(¢ + In 2]. (23) 


-y denotes the digamma function,” and ¢ is Euler’s 


~ McGraw-Hill Book Co., Inc., 


constant. 


eee 


5 
Bo 


If the gap width d approaches zero, (23) reduces to 
cos yok = Y0/Bo 


which is just that obtained by assuming only a single 
(lowest) mode in each region, and then equating these 
fields. 


2 A. Erdelyi, ef al., “Higher Transcendental Functions, Vol. I,” 
New York, N. Y., p. 15; 1953. 
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(Jvo — w)(Fvo — ao) IIa (w)TM2(770) 
(18) 


Differentiation of (22) with respect to 8» shows that 
the maximum allowable value of Bod is 7. When this 
value is substituted into (23), the cut-off value of the 
tooth depth h, is obtained: 


a d 
— — —|n 2 — Order (y0"d?). 
270 T 


he = (24) 


Thus /, is a function of the wave propagation direction, 
For propagation normal to the corrugations, the cut-off 
depth is exactly that given by Hurd,’ but as the propa- 
gation angle moves from the normal position, h, in- 
creases if d remains fixed. 

For propagation across a surface which has constant 
values for the tooth depth and spacing, there is no cut- 
off angle for surface wave propagation. This is illus- 
trated in Fig. 2, where v/k vs @ is plotted for d/A=0.1 
and for d/A=0. The upper curves for h/A=h-/d (at 
@=0 degrees) gives the maximum values for v/k which 
can be obtained with the given tooth and gap dimen- 
sions. Any h/X less than this limit value gives a corre- 
spondingly lower wave velocity ratio as indicated by 
the curves for 4/A=0.159. 

The normalized attenuation factor ao/k is plotted 
against the scan angle in Fig. 9. As is to be expected 
from the variation of v/k with 6, the attenuation of the 
“trapped” wave also decreases as the direction of propa- 
gation moves from the normal (@=0 degrees). Thus as 
the free space to surface wave velocity ratio decreases, 
a smaller percentage of the power in the surface wave 
is concentrated close to the corrugated surface. 

The mode amplitudes can be found from (15) and 
(18) and an approximate solution of these equations 
gives 


Tv 


d(an — Bo) 
a 21 | 


! and and In 2 
yod (an + a0) (On - 1) (Qn = a1) (2 = =) exp (-“>*) 


D(a, “1 
iran | Br | (an + jy)T E + mee | 1 E 


(25) 


Two ratios of particular interest have been computed 
and plotted (Fig. 10) for a specific corrugation geom- 
etry. These are B:/Bo| A1/Aol, the ratio of first higher 
order mode to dominant mode above the corrugations, 
and |4 ache the ratio of the dominant mode above to 
the dominant mode in the corrugations. At 6=90 de- 
grees the corrugations have no effect on the dominant 
mode which becomes a plane wave propagating with 
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Fig. 9—Transverse attenuation constant vs scan angle. 


the velocity of light in the z direction. The higher order 
modes, however, are highly attenuated in y, and their 
amplitudes are small compared to the amplitudes of the 
dominant mode. 


CONCLUSIONS 


Theory has predicted that a surface wave can propa- 
gate across a corrugated surface at an oblique angle, 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


1.00 


d/X = 0.10 
h/d = 0.159 


0.60 


6, DEGREES 


Fig. 10—Amplitude ratios vs scan angle. 


and that the wave velocity increases smoothly to the 
free space value as the wave direction approaches that 
parallel to the corrugations. Experiments have shown 
that a corrugated surface antenna may be scanned 
through an appreciable angle, at least +30 degrees, 
with the surface maintaining control over the elevation 
pattern. Thus an additional degree of flexibility has 
been added to the corrugated surface antenna. 
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Measurements of the Bandwidth of Radio Waves Propagated 
by the Troposphere Beyond the Horizon* 
J. H. CHISHOLM}, L. P. RAINVILLE}, J. F. ROCHE}, anv H. G. ROOT} 


ceived well beyond the horizon are propagated 
by a tropospheric mechanism which produces in 
effect a multiplicity of equivalent plane wave compo- 
nents arriving simultaneously at the receiving antenna 
with random phase relationships and angular directions 
_ which vary rapidly in time. The description of the com- 
munications capacity of this propagation mechanism 
has been referred to as “bandwidth” which has an ex- 
_ plicit meaning to communications engineers for circuit 
_ applications based upon parameters of resistance, capac- 
‘itance, and inductance which do not vary in time. If 
‘the troposphere could be “frozen,” an instantaneous 
transmission response could be obtained by transmit- 
ting an infinitely narrow impulse and using an infinitely 
wide bandwidth receiver or tuning a CW transmitter 
over the spectrum of frequencies of interest and receiv- 
ing the signal with a narrow band receiver. In order to 
approximate the latter technique in the case of the actual 
troposphere, a series of consecutive frequency sweeps 
over a finite band of frequencies with a period which is 
small in terms of the fading periods should show the 
near instantaneous successive transfer function H(w) 
of the mechanism. 
We are presenting a brief summary of experiments 
designed to investigate the transmission characteristics 
as a function of time and frequency using a frequency 


ean is general agreement that radio waves re- 


* Manuscript received by the PGAP, May 29, 1958. The research 
reported in this document was supported jointly by the Army, Navy, 
and Air Force under contract with the Massachusetts Institute of 
Technology. : 

t Lincoln Lab., Mass. Inst. Tech., Lexington, Mass. 


modulated 2290-mc klystron transmitter and a wide- 
band traveling-wave tube receiver over the 188-mile 
path between Round Hill, Mass. and Crawfords Hill, 
N. J. The transmitter frequency was swept over a 20- 
mc band at a rate of 100 cps and the results were re- 
corded with a 16-mm movie camera at a rate of 24 
frames per second. A sample sequence of motion pictures, 
shown in Fig. 1, summarizes the results of these experi- 
ments. The amplitude appears to vary in a systematic 
manner. The variation of amplitude as a function of 
time and frequency is shown simultaneously in the form 
of a three-dimensional model, Fig. 2, constructed from 
1500 frames of the aforementioned motion pictures and 
indicates variation of amplitude for a one-minute period 
over the range 2283-2297 mc. 

Analyses of these data in both the time and the fre- 
quency domains were performed by comparing the 
spacing of the maxima in time and frequency obtained 
from experimental data with the spacing predicted 
with a given scattering model.! The results of the time 
analysis indicate a mean separation of 2.2 seconds com- 
pared to a separation of 0.5 second predicted theoreti- 
cally. Conversely, the result of analysis of data in the 
frequency domain indicated a separation of maxima of 
3 mc which was greater by a factor of about 2 than the 
separation of 1.7 mc predicted by theory. A more rigor- 
ous analysis in the frequency domain consisted of deter- 
mining the correlation of amplitude as a function of 
spacing. The results of this analysis indicated a correla- 

1S, O. Rice, “Statistical fluctuations of radio field strength far 


beyond the horizon,” Proc. IRE, vol. 41, pp. 274-281; February, 
1953. 
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Fig. 1. 


tion coefficient of 0.4 at 2 mc separation compared with 
a 0.02 correlation coefficient predicted theoretically. 

Since it is not realistic to represent bandwidth as a 
simple figure when applied to tropospheric scatter 
transmissions, a statistical form is used. The maximum 
variations of amplitude within frequency increments of 
1.0 to 10 mc were determined and the distributions com- 
piled (Fig. 3). 

The results of the analyses in the time and frequency 
domains are consistent and both suggest that the multi- 
path delays are less than predicted by the scattering 
model used by Rice in his theoretical paper. 
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Fig. 3—Distribution of maximum variation of signal. 


Remarks on the Fading of Scattered Radio Waves* 


RICHARD A. 


an expression was derived for the envelope fading 
rate of radio waves scattered by dielectric turbu- 
lence, under the tacit assumption that the scattering 
volume V is effectively at infinity as seen by the trans- 


‘le a previous paper,' hereinafter designated as A, 


* Manuscript received by the PGAP, June 19, 1958. The research 
reported in this paper has been sponsored by the Electronics Res. 
Directorate of the AF Cambridge Res. Ctr., Air Res. and Dev. 
Command, under Contract AF 19(604)3495, 

} Inst. Math. Sci., New York University, New York 3, N. Y. 

* R. A. Silverman, “Fading of radio waves scattered by dielectric 
turbulence,” J. Appl. Phys., vol. 28, pp. 506-511; April, 1957; errata, 
ibid., p. 922; August, 1957. With a few exceptions, we adhere to the 
notation and terminology of this reference, 
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mitter and receiver, 7.e., that V is accurately character- 
ized by a single scattering angle. Since this condition is 
rarely met in practice, the fading rate formula derived 
in A has only limited applicability. The purpose of this 
note is to extend the treatment of A by including the 
additional contributions to fading which result when 
there is appreciable scattering angle variation within ¥; 
as well as when there is appreciable shear of the mean 
wind within V. 

As in A, we begin by dividing V into a large number 
N of cubical subvolumes Vi, - ++ Vy, each of side d. 
We denote the center of each V, by P,, and we associate 


a aanaiiell dthdiierntned, aitacin itil le tdci 


mR tn hide 
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with every P, its own scattering angle 0, and scattering 
wave vector’ K,. In addition to assumptions 1-4 and 6 
of A, we make the following two assumptions: 1) The 
scattering wave vectors associated with the different 
points of any subcube V, differ from each other by so 
little that we can regard every point of V, as having the 
same typical scattering wave vector, which we choose to 
be K,,; 2) The quantity d is such that d< Zp, the macro- 
eddy dimension, and K,d>>27 for 1<n<WN. Then, if 
the incident radiation has angular frequency w, by a 
standard approximation’ the contribution to the total 
scattered current from the dielectric turbulence in V, 
is given by the real part of 


T,(t) = Bk? exp (—iwt) exp (in) (Sn8n'/dndn’) 


‘ 


. f exp [i(Kn-tn)|Ae(rn, t)drn, (1) 
where the position vector r, is referred to P, as origin, 
and Ae(r, ¢) is the random departure from its mean of 
the relative dielectric constant at the point (r, ¢). Here 
d, is the distance from the transmitting antenna T to 
P,, and g, its gain in the direction 7P,; d,’ and g,’ are 
the corresponding quantities for the receiving antenna. 
As we shall see immediately, the values of the constants 
B and y, do not concern us here. 

We now assume that* 


(Im(t)In*(t!)) ~ 0, if m#n and |t—?t| <7, (2) 
where Tp is of the order of the correlation time of the 


received scattered current. To justify (2), we note that 
assumption 2 of the preceding paragraph implies that 


the “scattering eddies” are small compared to the cube 


size d. Thus, if as usually assumed, the fine structure of 
dielectric turbulence is such that the small eddies are 
substantially uncorrelated when separated by distances 
large compared to their size, and if in time 7» (the great- 
est time separation of interest) the dielectric noise is 
convected through a distance small compared to d, then 
(2) should be a good approximation.® Relying on the 
same kind of heuristics, we can establish the formula’ 


| EL exp [i(Kn- tn) ]Ae(rn, l)drn 
Va 
f exp [—i(Kn- rn’) |Ae(rn’, Ndr.) 


n 


~ ((A6)*)La8 f “Role, 7) exp [i(Ks* 0) ldo 


—2 


= (2m)*((Ae)?)Lo°Pn(Kn; 7), 


2 Defined as & times the difference between the unit vector from 
the transmitter to P, and a unit vector from P, to the receiver, 
where & is the wave number of the incident radiation. ; : 

3 See, e.g., R. A. Silverman, “On Radio Scattering by Dielectric 
Turbulence,” New York University, Inst. Math. Sci., Div. of Elec- 
tromagnetic Res., New York, N. Y., Res. Rep. EM-98; September, 
1956. 

4 The angular brackets denote the ensemble average, and the 
asterisk denotes the complex conjugate. ; : 

5 These conditions seem to be met in practice, at least at the higher 
of the carrier frequencies used. 
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where R,(0, 7) is the space-time correlation function of 
the dielectric noise as measured in V, (the variables 9 
and 7 are defined by ep=r—r’, r=t—-?’), and the quan- 
tity ®,(Kn, 7) is the semitransform used in A. Finally, 
using (1)—(3), we find that the correlation function of 
the total scattered current J(t) = > Ae I,,(¢) is given by 


C(r) = (II*(“))/(| TO |?) 


N N 
2 En exp (—ian) (Ky, 2) / Y Ena(Kyy 0), (4 
n=1 


I 


HT 


n=1 


where En = (gngn//dndn’)*. 
To take into account the effect of convection, we re- 
call (7) of A, to the effect that 


&,(k, r) = exp [i(k-vn)7]®,(k, ah (5) 


where ®,(k, 7) is the semitransform of the dielectric 
noise as measured in a coordinate system which in each 
V, moves with the velocity v, that remains after veloc- 
ity fluctuations on the scale of the smaller eddies (in- 
cluding the scattering eddies) have been smoothed out. 
(It is to be expected that this smoothing will have little 
effect on the univariate velocity distribution, which 
is approximately Gaussian.®) In A this moving coor- 
dinate system was called the “Lagrangian system,” but 
since the small eddy components of the velocity field 
have been averaged out, the term is somewhat of a mis- 
nomer. The possibility of separating the self-motion of 
the smaller scattering eddies from their motion due to 
convection by the larger eddies of size d, at least in an 
approximate way, is essential to our approach, and 
seems intuitively reasonable; however, the legitimacy 
of this separation, which is also implicit in the work of 
Kolmogorov,’ has recently been challenged.* We now 
substitute (5) in (4) and obtain 


N 
C(t) = 5G exp (—iwr) exp [i(Kn-Vn)t]¢(Kn, 7), (6) 
n=l 
where $(k, 7) =®,(k, 7)/®,(k, 0) is the normalized 
semitransform in the moving system. The weights Gn 
are defined by 


N 
Gn = En®(Kn, 0) / 21 En®s(Kn, 0) 


n=1 


and sum to unity. 

Since the wind v, has a turbulent component, the 
quantity C(r) is still random. Thus, if we center v, about 
its mean value by writing v.=v,+vn’, it remains to 


6 G. K. Batchelor, “The Theory of Homogeneous Turbulence,” 
Cambridge University Press, Cambridge, Eng., ch. 8; 1953. The effect 


‘of the small eddies on the bivariate velocity distribution is crucial, 


however (ibid.). net 

7 A_N. Kolmogorov, “The local structure of turbulence in incom- 
pressible viscous fluid for very large Reynolds number,” C. R. 
(Doklady) Acad. Sci. URSS, vol. 30, pp. 301-305; 1941. ; 

8 R.H. Kraichnan, “Irreversible statistical mechanics of incom- 
pressible hydromagnetic turbulence,” Phys. Rev., vol. 109, pp. 1407— 
1422: March, 1958; “The Structure of Isotropic Turbulence at Very 
High Reynolds Numbers,” New York University, Inst. Math. Sci., 
Div. of Electromagnetic Research, New York, N. Y., Res. Rep. 
MH-9; February, 1958. 
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average C(r) with respect to the random variables 
v,’, an operation we denote by an overbar. The possi- 
bility of doing this averaging independently of the 
averaging we have already done over realizations of 
the dielectric noise seen in the moving coordinate sys- 
tem rests on the reasonable assumption that the veloci- 
ties of eddies of size d are substantially independent of 
the dielectric fluctuations of the smaller scattering 
eddies. Assuming that the turbulent velocity field is 
homogeneous, so that the random variables v,’ are 
identically distributed, we obtain 


C(r) = DG, exp (—iwr) 


n=1 


“exp [i(Kn Vn) 16(Kn, T)W( Kn, T)\. (7) 


The quantity W(k, 7) is defined by 


Wk, ») =exp [eve] = fi exp lite-vyrlov)av', 8) 
where p(v’) is the probability density of the common 
distribution of the v,’; p(k, 7), like ¢(4, 7), is a correla- 
tion function in 7 for every k. As already remarked, 
p(v’) is effectively the univariate probability density of 
the turbulent velocity field, which is approximately 


- Gaussian. Hence the correlation function p(k, 7) is 


real, and is in fact itself approximately a Gaussian 
function. 

The calculation of the envelope fading rate from (7) 
proceeds in standard fashion. First, we note that since 
I(t) contains only negative frequencies, the envelope of 
its real part is just | Z(2)| . Next, we note that since J(t) 
is approximately a Gaussian process, | aa (t)| is approxi- 
mately a Rayleigh process. It follows that the rate at 
which | (2) | crosses the level o\/<|JI|?> is given by® 


Z(c) = 2\/O/ne exp (—o”), o > 0. 


®V. I. Bunimovich, “Amplitude peaks of random noise,” Zh. 
Tekh. Fiz., vol. 21, pp. 625-636; June, 1951. 
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The fading parameter Q is defined by 
Q = — RO - SOP = [-@/4)|C@)| 1-0 9) 


where R(r) and S(r) are the real and imaginary parts 
of C(r), and the dot denotes differentiation with respect — 
to tr. We see from (9) that any unimodular complex ~ 
factor of C(r) can be replaced by unity without chang- — 
ing the envelope fading rate. Thus, we are at liberty to © 
set equal to unity the RF factor exp (—iwr) appearing — 
in (7). 

The remainder of the calculation is straightforward — 
and consists in calculating the fading parameter Q from 
(7)-(9). The result is 


N 


i—1 n=1 


-| LetK. | + 0 GKny 


n=1 


(10) 

n=1 
It follows from (10) that the mean wind has an effect -_ 
on the envelope fading rate only when different parts 
of the scattering volume are characterized by appreci- 
ably different scattering wave vectors or when there is 
appreciable shear in the mean wind within V. To see 
this, we note that the second and third terms of (10) 
cancel each other when we set K,=K»y and v,=vo, 
1<n<WN. Thus these terms represent the effects of scat- 
tering angle variation and wind shear omitted in A, 
where it was assumed (rather unrealistically) that the 
effect of variation of K and v within V is unimportant. 
The first and last terms of (10), which were retained in 
A, represent the fading due to the random self-motion 
of the scattering eddies as seen in a coordinate system 
which effectively moves with the local wind velocity and 
the fading due to the convection of the scattering eddies 
by the mean wind and the macro-eddies. Finally, we 
note that for some purposes it may be convenient to 
approximate the sums in (10) by integrals over V, with 
terms which are the continuous analogs (defined in the 
obvious way) of G,, K, and Vn. 


Phantom Radar Targets at Millimeter 
Radio Wavelengths* 
C. W. TOLBERT}, A. W. STRAITON}, anv C. 0. BRITT 


Summary—This paper describes the techniques and the meas-_ 


ured radar return from ‘‘phantom targets” using 8.6 and 4.3- 
millimeter radars. The radar returns are compared to the measured 
back-scattering cross section of water drops, insects, steam and 


_ * Manuscript received by the PGAP, February 21, 1958; re- 
vised manuscript received, June 16, 1958. 
} Elec. Eng. Res. Lab., University of Texas, Austin, Texas. 


other materials at 8.6-mm wavelength as measured by a CW radar 
at this wavelength. 

Within the limits of the conditions used in the laboratory, it was 
impossible to produce returns from synthesized refractive index 
gradients of sufficient magnitude to account for those noted on the 
radar. It is concluded, therefore, that for the millimeter wavelengths 
and the short ranges considered, the observed phantom returns were 
due to solid or liquid particles in the atmosphere. 
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INTRODUCTION 


ANY observers have noted radar returns from 
phantom targets at frequencies throughout the 
radio spectrum. These observers have accumu- 

lated a large amount of data on the characteristics of 

the return from such targets and attempts have been 
made to simulate conditions hypothesized as respon- 
sible for the radar return.! 

A recent and most complete resume of the data per- 
taining to the characteristics of these phantom radar 
returns over a wide range of frequencies has been made by 
Plank.? In this summation and in other reports, a large 
number of the returns have been attributed to refrac- 
tive index gradients. The strength of the refractive index 
gradients and the size and shape of the atmospheric in- 
homogeneities which would produce these returns, 
which are called “angels,” were considered by Plank. 
As a general rule, the measured returns are from one to 

two orders of magnitude greater than have been cal- 
culated for the angel models at short centimeter and 
‘millimeter wavelengths. The index gradients used in 
the angel models are based on maximum refractive in- 
dex measurements as reported by Crain,? Cunningham* 
and others, and are assumed to exist over a sufficiently 
large area to give the effect of an infinite plane. 

In order to obtain data on the simultaneous occur- 
rence of phantom returns on 2-mm wavelengths, 8.6- 
‘mm and 4.3-mm radars were operated simultaneously. 
The sensitivity of the radars used was such that most 

of the observations were restricted to the first 1000 
meters above the antenna. In addition, return from a 
rocket exhaust was measured. 

To complement the radar return measurements, 
‘back-scattering cross sections of objects that might be 
airborne were measured with CW techniques.’ The ob- 
jects measured included insects, grass, water droplets, 
soap bubbles filled with helium, and water vapor jets or 
columns of steam and the fog produced with the at- 
mosphere by cooling with dry ice and liquid nitrogen. 


~ INSTRUMENTATION AND MEASURING TECHNIQUES 
OF THE RADAR OBSERVATIONS 


The 8.6-mm and the 4.3-mm radars are shown in Fig. 
4. Bistatic antennas were used for transmitting and re- 
ceiving at each frequency. The effective beamwidth 
was 1 degree for the 8.6-mm radar and 0.5 degree for 


1A. B. Crawford, “Radar reflections in the lower atmosphere,” 
Proc. IRE, vol. 37, pp. 404-405; April, 1949. 


2V. G. Plank, “A Meteorological Study of Radar Angels,” Geo- ~ 


physics Res. Paper No. 52, Air Force Cambridge Res. Ctr., Cam- 
bridge, Mass.; July, 1956. a 

8 J. R. Gerhardt, C. M. Crain, and H. Chapman, “Microwave re- 

- fractive index fluctuations associated with convective activity in the 

atmosphere,” Bull. Amer. Meteorological Soc., vol. 37, pp. 251-262; 


June, 1956. , 
4R. M. Cunningham, V. G. Plank, and C. M. Campen, “Cloud 
_refractive index studies,” Geophysics Res. Paper No. 51, AFCRC- 
-TR-56-210, ASTIA No. AD110259, Air Force Cambridge Res. Ctr., 
‘Cambridge, Mass.; October, 1956. - : 
i 5 C. W. Tolbert, et al., “Back-scattering cross sections of water 
droplets, rain and foliage at 4.3 millimeter radio wavelengths, 


EERL Rep. No. 91, Univ. of Texas; April 30, 1957. 
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Fig. 1—View of millimeter radar. 


the 4.3-mm radar. These antennas were pointed at a 
common volume. The sensitivity of the two radars was 
adjusted to produce the same deflection on an “A 
scope” per unit back-scatter cross section at each wave- 
length when the antennas were pointed at a point 
source. This permitted matching the sensitivity of the 
two systems in the direction of the main beam of the 
antenna patterns, but did not provide matching in other 
directions. 

Antennas with the same effective beamwidth would 
have been very desirable, but were not considered essen- 
tial since the emphasis here is on the presence of phan- 
tom targets at the two wavelengths and on the meas- 
ured back-scattering cross section of possible airborne 
reflectors and simulated refractive index gradients. 

The pulse length of the radars was 0.15 microsecond 
and the recovery from the “main-bang” was effective 
within 250 feet of the antennas. The two radars were 
adjusted to “fire” alternately, with the return from the 
4.3-mm system producing a positive return on the A 
scope and the 8.6-mm system producing a negative re- 
turn on the A scope. The face of the A scope was photo- 
graphed with a motion picture camera at 24 frames per 
second. 

Measurements also were made of the back-scattering 
cross section of the exhaust from a simulated rocket 
motor burning a 6 to 5 ratio of zinc and sulfur as fuel. 
The stem or straight high velocity portion of the rocket 
exhaust filled the antenna pattern between the half- 
power points, and the plume or flared portion more than 
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Fig. 2—Samples of radar recordings on “A scope.” (a) Photograph 


of A scope showing scales. (b) Negligible return on both wave- 
lengths. (c) Single return from 750 feet noted on both wave- 
lengths. (d) Returns from two elevations with 4.3-mm return 
suppressed. (e) Returns from two elevations with one 4.3-mm 
return suppressed. (f) Returns from two elevations with 8.6-mm 
return suppressed. (g) Multiple returns seen alternately on each 
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(k) (1) (m) 


wavelength. (h) Returns from two elevations with 8.6-mm re- 
turn suppressed on one. (i) Multiple return noted on both 
wavelengths. (j) Multiple return from several elevations with 
some suppression on 4.3 mm. (k) Volume return centered at 500 
feet with joint return to 800 feet. (1) Similar to (k) with 4.3-mm 


return suppressed. (m) Volume return following inverse square 
law. 
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Corner Reflector (4.6cm.) 


Rain (10 mm/hr.) 


Polyfoam Sphere (|.22 meter dia) 


Rubber Balloon (| meter dia.) 


Rocket Exhaust (Stem) 


Rocket Exhaust (Plume) 
Johnson grass (green,0.40 cm.) 
Phantom Target Moth (dead | day) 
8.6 mm. (max.). a alive 
Phantom Target ater droplet (3mm.dia.) 
4.3mm. (max.) 


Johnson grass, dry(0.40.cm.) 


6| Minimum _ detectable House fly ,(dead 2 days) 
level for Radar , 
(Range 152 meters) 


mm. Water droplet 


Mininilendeteetabie Bubble, Helium filled (10cm. dia) 


level for solid model. 
(Range .045 meters) 


BACK SCATTERING CROSS SECTION (square meters) 


Fig. 3—Measured back-scatter characteristics at 
8.6-mm wavelengths. 


filled the pattern between the first nulls. The escaping 
gases and particles were semiopaque to light in both the 
stem and the plume during the period when radar re- 
turns were detected. 


RESULTS OF THE RADAR MEASUREMENTS 


The simultaneous vertical soundings of the 8.6 and 
4.3-mm radars were made on nine days in January, 
1958. The days represented a variety of meteorological 
situations ranging from completely clear to heavily 
overcast. On no occasion did the radar fail to pick up a 
phantom return within a matter of a few seconds. How- 
ever, there was _considerable range in the frequency 
with which the returns occurred. Samples illustrative 
of some of the conditions observed on the radars are 
shown in Fig. 2(b) through 2(m). The horizontal scale 


shows altitude in increments of 250 feet. The returns” 


measured are characteristic of those reported by other 
investigators with durations of up to several seconds, 
flutter characteristics of a coherent target, and occur- 
rences at more than one range simultaneously. As would 
be expected from the antenna patterns, more 8.6-mm 
returns were noted than 4.3-mm returns. 


In all, 3000 feet of film were taken, which will be 
analyzed in detail to obtain the statistical properties of 
the soundings. The analyses will be compared with 
the meteorological conditions under which they were 
made. This analysis is beyond the scope of this present 
work and will be presented in a later report. From 
these soundings, however, we may select the maximum 
return for comparison with the cross sections measured 
in the model range studies described in the next section. 
The scattering cross sections were calculated for the 
phantom returns and are included in Fig. 3. 

In a similar manner, the scattering cross sections for 
the rocket exhausts were measured for the column fall- 
ing within the antenna pattern, both for the stem and 
the plume and these results are also shown in Fig. 3. 

In previous tests, this Laboratory has had occasion 
to measure the cross section of a 4.6-cm corner reflector, 
rain® as a function of rate, individual water drops as a 
function of size, and polyfoam spheres and rubber bal- 
loons which have been inflated with air. For the sake of 
comparison with the results of the present work, sam- 
ples of these previous measurements have been added 
to Fig. 3. 


INSTRUMENTATION AND MEASURING TECHNIQUES 
OF THE MODEL STUDIES 


A continuous wave radar operating at 8.6-mm wave- 
length was set up for laboratory measurements. The 
system had a minimum point back-scattering cross sec- 
tion detecting level of 10-7 square meters at a distance 
of 0.45 meter from the antenna. The 0.45-meter dis- 
tance was in the Fraunhofer region for the antenna used. 

Measurements were made of the back-scattering cross 
section of water droplets, insects, and grass. The ob- 
jects were allowed to fall through and/or were sus- 
pended in the antenna pattern. The water droplets were 
produced with hypodermic needles. The size of the 
drops was controlled by varying the water pressure and 
needle size.® 

A water droplet fog condensed from the air with dry 
ice and liquid nitrogen was allowed to pour down across 
the volume illuminated by the antenna. Soap bubbles 
were filled with helium, steam, and water vapor, and al- 
lowed to rise through the illuminated volume in front 
of the antenna. The distance from the antenna in these 
measurements was such that the minimum detectable 
reflectivity was 10-!° square meters. The reflected sig- 
nals were recorded on a Brush recorder. The scattering 
cross sections calculated for each of the measurable 
cases are plotted in Fig. 3. 


SUMMARY OF CROSS-SECTION MEASUREMENTS 


The fog droplets produced by cooling of the air and 
the soap bubbles filled with water vapor and steam gave 
no measurable signals and therefore had a power re- 
flection coefficient of less than 107. 
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Single drops of water, 2 mm in diameter, produced a 
small but measurable return of 2.5 X 1077 square meters, 
while a drop 3 mm in size produced a return of 7.5 X107° 
square meters. It is noted that the 3-mm drop gives ap- 
proximately the same return as a medium-sized house 
fly. In Fig. 4 this similarity has been used to compare 
the return measured for the water drop or house fly to 
returns calculated from assumed refractive index gradi- 
ents in the atmosphere. 

The steam presented the case of the sharpest artificial 
gas boundary with a difference across the surface of ap- 
proximately 150 NV units. It is difficult to imagine any 
index changes in the atmosphere except for solid or 
liquid particles with gradient as great as this. Since no 
return at the millimeter wavelength was measurable 
at the steam boundary, it is felt that it would be very 
difficult to get a return on millimeter radar from refrac- 
tive gradients. 

The various targets represented in Fig. 3 provide a 
comparison between reflectors ranging over six decades 
in scattering cross section, and thus provide a broad 
scale which helps to orient the cross section of the vari- 
ous materials considered. 


BACKSCATTER AS A FUNCTION OF RANGE 


The potential return from a single millimeter water 
drop or a common house fly is plotted in Fig. 4 asa 
function of the range from 100 meters to 2000 meters. 
Two scales are used for the return indication; the one 
on the left is the ratio of the power transmitted to the 
power received, normalized to 100 meters, and the one 
on the right is the power reflectivity of a plane discon- 
tinuity which would produce the same return at the 
same range. This second scale also may be expressed in 
terms of the refractive index discontinuity which would 
be required to produce this return. The ordinates cor- 
responding to AW of 6.3 and 20 are indicated. It is noted 
that a smaller index change would be required at a 
greater range since the assumption of an infinite plane 
reflector causes the power to decrease only as the square 
of the range, whereas the point reflector would have 
the power decrease as the fourth power of the range. 

The use of this scale is not intended to imply that a 
large plane reflecting surface in the atmosphere is real- 
istic. The extent of the flat plate for the reflections 
shown would need to increase with altitude and the 
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8.6-mm wavelengths. 


deviations from flatness would need to be small com- 
pared to a wavelength. It would seem very unlikely 
that these conditions would be met in the atmosphere 
for millimeter wavelengths at least. The plane reflector 
does present, however, an upper limit to the power that 
might be returned from a single index gradient. For 
the short ranges, below approximately 1000 meters, 
single water drops or insects would seem to present 
much better targets than refractive index gradients. 
For much greater ranges and longer wavelengths, index 
gradients would seem to play a more significant role. 
Since many of the meteorological conditions under 
which longer wavelength returns have been noted are 
conducive to the formation of water drops, the possibil- 
ity of this type of return should not be overlooked. 
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Telephone Remote Control Circuit for 
an Antenna Test Site* 
L. YOUNG} anp G. M. WARD} 


Summary—A circuit is described which permits the remote con- 
trol of several quantities, such as transmitter frequency, direction of 
polarization and RF match, at an unattended antenna tower. It’re- 
quires only a single telephone line, which already exists for com- 
munication purposes at many test sites. Each quantity is controlled 

by a separate motor. The motors are controlled separately from the 
control end. This is achieved by adapting and installing a dialing 
circuit which allows, first, any desired motor to be ‘‘dialed’’; next, 
to be controlled; and, finally, the line cleared so that another motor 
can be called. 

The parts for the dialing circuit were purchased at a total cost of 
about $150. Detailed circuit diagrams are given showing how the 


circuit is constructed. 

qt 

A tween the transmitting and receiving antennas 
is so great as to make the remote control of the 
system parameters (frequency, polarization, matching, 
etc.) at one of the locations desirable. No suitable re- 
mote control equipment seems to be available commer- 
‘cially, and none has been described for this specific ap- 
plication. A system built for an antenna test site at the 
Westinghouse Electronics Division is described in this 
communication. 

At most sites, a telephone line is available between the 
two antenna towers for communication purposes. The 
usual procedure is to instruct a man making the adjust- 
ments at the remote antenna. By adapting and install- 
ing a dialing circuit connected to several relays con- 
trolling motors, the remote tower can be left unat- 
tended and the several system parameters can be ad- 

_justed independently over the one telephone line. 


INTRODUCTION 


many antenna test sites, the separation be 


FREQUENCY REMOTE CONTROL 


The most commonly required adjustment at the re- 
‘mote site has been changing the frequency of the signal 
generator. The first remote control circuit, therefore, 
-merely used the two wires of the telephone line to drive 
a slow motor which was suitably geared down to tune 
“the remote signal generator. The frequency was meas- 

ured at the receiving antenna with the superheterodyne 
receiver by bracketing the transmitter frequency with 
the two local oscillator frequencies which peaked the 
signal. It was arranged so that the remotely controlled 
motor could be driven in either direction by pressing one 
of two buttons. 
The success of this arrangement led to more ambitious 
plans. It was desired to adjust by remote control not 


* Manuscript received by the PGAP, January 27, 1958; revised 


manuscript received, June 2, 1958. : 
+ Westinghouse Electric Corp., Baltimore, Md. 


only the frequency, but also the polarization, and a 
double stub tuner; the possibility of making other ad- 
justments was also contemplated. A means had to be 
devised which allowed any one of several motors to be 
controlled remotely, still using only the one telephone 
line. This was achieved by adapting and installing a 
dialing circuit which allowed any desired motor to be 
first “dialed,” and then controlled, and which finally 
cleared the line until another motor was called. Each 
motor makes one of the desired adjustments. 

The parts for the dialing circuit can all be purchased 
(total cost from surplus suppliers amounted to about 
$150), and can readily be constructed. 


DIALING AND CONTROL CIRCUIT 


The circuits at the two ends of the telephone line are 
shown in Figs. 1 and 2. Their operation will be described 
in three stages, according to the setting of the operations 
switch S1. 


Operation A 


The operations switch $1 shown in Fig. 1, the control 
end of the line, is set to position 1 (“dial”), thus con- 
necting the contacts of relay K1 across the input end 
of the telephone line, at terminals T2 and T3. The 
dialing contacts are normally open; they close as soon 
as the dial is moved off normal, then open and close 
successively as the dial returns from the number dialed. 
In Fig. 2, which depicts the remote end of the line, the 
relay K2 is normally released as shown. When the num- 
ber is dialed, relay K3 successively operates and releases, 
energizing and deenergizing the coil of stepping switch 
SS1 via the contacts of K3, thus moving the wipers of 
the two-bank ten-position stepping switch SS1 to the 
desired position. The contacts of one bank (B) are all 
tied together, except for the first contact, and con- 
nected to the slow-operate relay K2, which operates a 
fraction of a second after dialing has been completed. 
Its contacts prepare a new circuit arrangement de- 
scribed in operation B. 


Operation B 


The operations switch S1 (Fig. 1) is now set to posi- 
tion 2 (“operate”), connecting the telephone line ter- 
minals T2 and T3 across the momentary contact double- 
pole-double-throw key KS1, which is normally open, 
and which can connect the 48-volt de power supply 
with either polarity across T2-T3. Since relay K2 in 
Fig. 2 is now operated, and the stepping switch SS1 is 
set on bank A to connect through to the selected relay 
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Fig. 2—Telephone remote control (remote end). 


coil in the relay set (Fig. 3), the selected relay now oper- 
ates. This closes a circuit containing a motor and 
power supply, both of which are at the remote end, for 
as long as the key KS1 at the control end (Fig. 1) is 
operated. Because of the rectifiers in series with the re- 
lays, the polarity determined by the key KS1 controls 
the direction of motor drive by selecting the appropri- 
ate relay. The correct setting having been obtained 
with the selected motor, operation B is completed. 


Operation C 


It remains to clear the line so that the next motor can 
be dialed. The operations switch S1 (Fig. 1) is set to 
position 3. This shorts T2 and T3 and connects both to 
the 48-volt dc power supply. Since there is no voltage 
between T2 and T3, but only from either of them to 
ground, relay K4 (Fig. 2) operates. Its contact energizes 
the reset coil of SS1, via its off-normal contact ONG, 
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which was closed in operation A, thus driving all the 
wipers of the stepping switch SS1 back to position 1. 
At this time, the off-normal contact ONC1 opens, and 
the wipers of the stepping switch SS1 stop. At the same 
time, the circuit of K2 is opened on bank B and it re- 
leases. Everything is now set as it was before operation 
A. When the operations switch S1 (Fig. 1) is again set 
in position 1, the next motor may be dialed as described 
under operation A. 


SPEAKERS 


At both ends of the line, speakers are connected 
across terminals T4 and TS. The line is balanced (hence 
the two coils for relay K3), so that the speakers may be 
left permanently connected across the line through 
large capacitors and used at any time. 


CONCLUSION 


A circuit which uses only available and reliable com- 
ponents, and which can readily be constructed at low 
cost, has been described. It gives great flexibility in re- 
mote control with only a single telephone line by allow- 
ing any one of several motors to be dialed, then con- 
trolled and set, and, finally, released to select a new 
function. Each motor adjusts some quantity, such as 
frequency or direction of polarization, which can there- 
by be remotely controlled. 
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